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PREFACE 


Tlii? book treats, in a self-contained manner, the problem of the 
deformation of an clastic solid without making the assumption of the 
classical infinitesimal, or linear, theory of elasticity that the deforma- 
tion is .so small that the squares of the strain components are negligible. 
Matrices are used methodically, and an account of those properties 
of matrices which are used in the book is given in the first chapter. 
The importance of non-i.sotropic media is emphasized, and the fact 
that applied strc.ss can make non-isotropic a medium which is isotropic 
when free from stre.cs is pointed out. The relations between the 
third-order, a.s well as the second-order, clastic constants for various 
tyj)e.s of crystalline media arc derived. The agreement between the 
])rcdictions of theory and Bridgman’s experiments on the compressi- 
iiility of media under hydrostatic pressures up to 10*' atmospheres is 
shown. The final chapter applies the theory to the problem of the 
deformation of spherical shells and circular tubes under extreme 
extonial and internal pres.surcs. 

Wo hjvvc, so far, very little experimental knowledge of the third- 
order elastic constants or of the cfTcct of applied stress upon the second- 
order ola.stic constant.s. If the mathematical treatment given here 
servc.s to stimulate the procurement of experimental knowledge of 
thc.se phenomena we shall have attained our aim. 

Part of the material of the book was given in lectures at the Dublin 
In.'-tituto for Advanced Studios (summer, 1948), the Carnegie Institute 
of Technology (fall, 1918), and the Escola Tecnica do Exercito, Rio de 
.lanoiro (1919). We have thought it well to let the book maintain 
the elomontary form of treatment adopted in the.se lectures and have 
not attempted to provide rcfcrcnce.s to other treatments of the topics 
di.^cuKsod in thc'bpcjh.- Owing, doubtless, to technological demands, 
those topitxs are currently mqch.fJiscus.'^cd. and the mierested reader 
tm.iy find many references in the late.st issues of sUc'h abstract journals 
iiK Matha7m!ical Ihrmvs and Sekner Abstracts. 

'Hie inscription that closed the preface to my prcviotis book in the 
Wiloy Applied Mathcmntic.s S’orfoJCTTjYroductfoa ./o .d.7ip/fcd Mathc- 
inaU'cs, may well l-»o repeated here: • , 

'I'o the Glory of God, Honor of Ireland, 
and 

tiolidnrity of the Americas. 


t.'io Jc-.ri’-.i 


F. D. MriiNAGiiAX 
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VECTORS AND MATRICES 


Tlic tjipory of finite deformntions of im clastic solid is most easily 
presented and understood by the use of matrices. Very little beyond 
(he elements of the tlicorj’ of matrices is needed, but since many 
enpineerb who may be interested in the theorj' of finite deformations 
may not be familiar with matrices tins introductory chapter gives a 
self-contained and elementary account of that part of the theory of 
mat rices which is necessary for our purpose. 

1. The innlrK concept 

The concejit of a matrix arises naturally and inevitably when we con- 
sider several difTerentiablc functions of several variables. In ele- 
ment ar>' c!dculu.s wc learn that, if u = n{x) is a diffcrontiablo function 
of a single independent variable x, then the differential du of n is 
obtaincii by multiplying (he dilTercntial dx of x bj' the derivative Ux of 
» with re.spect to x: 

i/» = Uj dx. 

Thus the ilerivative Ux of )/ witli respect to x may be regarded as a 
tmgnificalion fnrinr (hat converts, by multiplication, the differential 
dx of X into the differential dti of w. Wlien wc pas,s to the considera- 
tion of lira differentiable functions (w, r) of tivo independent variables 
(x, »/) we Imve two dilYerentials (rftt, dv) that arc connected with the two 
dilTercntinls (dx, (hi) of the independent variables by the formulas 

du ^ }tx dr -r Uy dy. 

dr " Vx dx -h Vy dy. 

Here u, and Uy i-.re the partial derivatives of it = »(r, y) witli respect 
to X and y, re*,pectively, and c.- and tv arc the p.artial derivatives of 
j Hr, y) with n-spect to x and y, respectively. The collection of 
four partial derivatives (»-, tv, r,,) novr plays the role pre^douslv 
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plaj ed bi the one dem atn c Ux, and this collection of four denv ntives 
ma> be considered as a magni^cotion /ttcior that converts, by muUipli- 
cation, the pair of differentials (dx, dy) into the pair of differentials 
(dii, rft) R e regard the pair of differentials {dx, dy) as a single entitj’, 
the differential lector of the independent variables, and, similarly, v.e 
regard the pair of differentials (du, dv) as a single entity, the differen- 
tial vector of the dependent variables We denote these differential 
vectors b> bold-face type, thus dx \s the differential vector of the 
independent vanables and (dz, dy) are the coordinates of dx. In 
s> mbols 

dx “ t(di, dy) 


(the translation of this stenographic notation being dx is the vector 
whose coordinates are dx and dy) Similarlj , 

du = t (du, do) 


(translate this stenographic symbol) It is convenient to imtc the 
coordinates of dx m a column, i c one above the other, thus 


(S) 


Wo refer to this sjmbol os a 2 X 1 matrix, i e , a matrix of tno rons 
and one column The first coordinate dx of the vector dx is the 


element m the first row of the matnv 


©■ 


and the second coor- 


dinate dy of the vector dx is the element in the second row of the 
Similarly, nc write the coordinates of du os the clc- 


■C:) 


menta of a 2 X 1 motnx 


(i> 


and then the collection of four 


denvatnes (u„ u„ t,, i,) maj be regarded as constituting a mngnifica- 


(dx\ 

'W 


into the 2X1 matnx 1 


write the four derivatives as the 


iron factor that converts, bj muftiphcation the 2 X i matnx ( 
c(*;) Wc»r 

elements of a 2 X 2 matrix (to , n matnx of two rows ond two col- 
umns), and wc denote ibis matrix bv the svmbol — ^ 

(x, y) 

(u, r) 


(n. r) ^ /«, uA 

(x.y)’“V*'* ‘■‘'Z 
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Note that the upper letters of the S 3 Tnbol 

the columns, of the matrix 


tell the rows, 


whereas the lower letters tell the columns, of the matrix. Thus 
the elements of the first row are the derivatives of u with respect to 
X and y, and the elements of the second row are the derivatives of v 
with respect to x and y. On the other hand, the elements of the first 
column arc the derivatives of and v with respect to x. and the ele- 
ments of the .second column are the derivatives of v and v with respect 

y. This 2X2 matrix is known as the Jacobian matrix 


In 


(3-, y) 


(after K, G. .1. Jacobi [1801-1851], a German mathematician) of the 
pair of functions (7/, v) with respect to the pair of independent variables 
(x, ?/); we denote it bj' the sj'mbol ./. We now introduce a shorthand 
notation in which we denote anj' one-column matrix by the symbol 
for the element in its first row; thus we denote the 2X1 matrix 

simplj’ by f/x and the 2X1 matrix simply by du. In 
this notation the pair of equations 

du = Uj dx Uy dy, 
dv — 7v dx -1- Vy dy 
appears iti the nbljrcviatod form 

du — J dx. 


'I'hc 2X2 matrix ./ has two 2 X 1 matrices, or cohmn vectors as we 

hhall term them, namely, and and two 1 X 2 matrices, 

or row vectors ns wo shall cull them, namely, (j/^, and (j'x, Vy). We 
term the c\prp«sion dx -f a,, dy (which furnishc.s the first element of 
the 2 X 1 matrix du) the product of the 2 X 1 matrix dx by the finst 
row vector of J; similarly, the expression v^dx + Vydy (which fur- 
ni.shc-! the .second element of the 2 X 1 matrix dit) is the product of 
the 2 X 1 matrix dx by the .‘^cond row vector of J. In this ter- 
minology. then, the product J dx of the 2 X 1 matrix dx by the 2X2 
Jacobian matrix./ is a 2 X 1 matrix (/« whose elements are obtained by 
methodically multiplyinp; the 2 X 1 matrix dx by the row vectors of 
J, one after the other. 

Example. The formulas conncctinR plane polar coordinates (p, <}>) 
and rectatiKuhir Cartesian coordinates (x, y) arc x - p cos <;b, w - p sin 6. 

The Jarobian matrix of (x. y) with respect to (p, < 3 ) is 
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=(: 


'cos <tt —p sin 
i^stn ^ p cos 


:) 


and the product of the 2X1 matnx dp = hy J la 


dx 


_ / cos ^dp p sm if) dip 
\sm ifkdp -i- p cos ^ dp^ 


:) 


The first element of this matnx is the differential of x and the second 
element is the differential of y 

Remark "Vou should note the analogy between the relation du ~ 
J dx and the basic relation dy = y^dx of elementary differential 
calculus Realize clearl} that it is the 2X2 Jacobian matrix J that 
plajs for two differentiable functions of two independent vanables 
the role played by the denvcdiie of one differentiable function of one 
independent variable The four elements of J have indiMdually no 
fundamental ei^ificancc It is the 2X2 Jacobian inalnx of which 
they are the elements that is the important thing 
There is no difficulty in extending the theory of the preceding para 
graphs to three differentiable functions of three independent vanables 
or generallj to n differentiable functions of n independent variables 
n " 3 4 5 Thus for three mdependent variables we denote 

/di\ 

bj dx the 3X1 matnx I dy | and b> du the 3 X 1 matnx I di | 
\(i2 / \du)/ 

T1 e Jacobian matnx J = • — is the matnx 
{x y z) 

( M, Wy 

p. V, p. 1 

IT, Wy Wf/ 

of which for example the 1X3 matnx («, Uy u,) is the first row 
vector J maj be regarded ns a magnification factor that converts 
bj multiplication the 3 X 1 matnx rfz into the 3 X 1 matrix du 

f/« ^ J dx 


In this kind of multiplication of di by J each element of the 3X1 
matnx du is the product of the 3X1 matnx dx b> the corresponding 
row \ector of J for example the element dt> in the second row of (he 
3X1 matnx du is gn cn b^ the formula 


dp = t,rfz + Tyrfp + t* dz 
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^•:xnlnI)lc. The formulas connecting space polar coordinates 
(r, 0, if>) and rcclnngular Cartesian coordinates (x, y, z) are 

X = r sin 0 cos 4>, 

y ~ r sin 0 sin <f>, 

z ~ T cos y. 


'i'lie .Taeolnan matrix ./ 


of (i, y, z) wth respect to (r, 0, (f)) is 

(r, 0, <f>) 


J - 


rsin 0 cos 
sin 0 sin 6 
cos 0 


r cos 0 cos p 
r cos 0 sin 4> 
— r sin 0 


■ r sin 0 sin <^>^ 
T sin 6 cos <t> 
0 


Tlic elements of t lie 3 X 1 matrix dx arc furnished the formula 
i!x ./ f/r; thus 


dx ~ sin 0 cos lii dr + r cos 0 cos <^d0 — r sin 0 sin (p dif>, 


dy ^ sin 0 sin dr -f- r cos 0 sin dO r sin 0 cos 4> d<l>, 


dz = cos 0 dr ~ r sin 6 dO. 


It should lie clear from the preceding discussion wliat is meant by 
an 111 X n matrix -1 wlicrc m and n are any two positive integers. This 
is a coih'ction of inu numbers arranged in in rows and n columns. We 
ilenote the ilnnrnt in the rth row and .stli column by the syml>ol 
«/. Note carebjlly that the upper label tells the row and the lower 
label tells the cohimn. In the applications of matrix theory which we 
shall have to make neither in nor it will surpass 3, and we can use the 
readily vi,-ualiz<'d concepts of onc-diinensional vectors (i.c., vectors on 
a line) or of two-<limcnsional vectors (i.e.. plane vectors) or of thrcc- 
<iiinert;iona! voe(oi>' (i.e., space vectors). For e.xample, if in = 2, 
n 3, so that .1 has two rows and tlm'c column.^ .4 possesses three 


plane column vectors of which the second, for instance, is and 

\«2V 

it aPo t\vo .vpacc row vectors of which the first, for instance, is 

bti’i oj*. uj*!. If U is an n X 7 > matrix (so that B possesses a.s many 
tov.\ SIS .1 po'se.-«c>> columns) it is po'^siblc to multiply nuy column 
vector of /? i>y any row vector of el .‘.ince the column vectors of B and 


?!ie row vectors of .1 arc .all of the same dimension, namely, n. We 
orite tlie-c variou*: prcuiucis of column vectors of B by row vectors 
of A ns an »?! X p matrix (\ tb.o clement in the jlh row and /.th column 
of r iHunsj the pr.viuot <n the /.th column vector of B by the jth row 
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vector of A, ; «= 1, m, fc ** 1, , p, and we term C the 

product of the n X p matrix B bj the m X n matrix A We indicate 
this product operation as follows 

C = AB 


EXERCISES 


1. Shon that if A is polite in multiplication, i e , if AB is the same as B for every 
B, A must be square (i e , ha\e the same number of columns as rows) with all its 
diagonal elements 1 and all its other elements zero 

A ole The diagonal elements of a square matrix A arc the elements oih os* 
whose row and column numbers are equal 

Hint I( A 13 an m X n maliix and 6 an n X p matrix, AB is an m X p matrix 
For A B to be the same as fine roust have tn = n bo that A is square 
Remark The matrices that are polite m multiplication arc known as unit 
matrices Terming the common number of rows snd columns of a square matrix 
its dimension, we denote the unit matrix of dimension n by the symbol Er Thus 



2 Form the following products 

-Dit -!)• 1 )(^ -?). 


(C) 




(«) ( “ m (5 ‘f). w ^^2^, 


A ole Tlie tlilTrrcnce between the answers to (a) and (6) shows that when both 
products AB and BA can be formed (show that each product is a square matrix 
when this happens) and when they have the same dimension (show that A and B 
arc each square when this happens) the products arc not necessarily the same 
MullipUcalion of square malriees of lAe same dimension tsnot, in general, commiUatiie 
The results of (c) and (cO show that if we denote by cB, the n X n matrix whoso 
diagonal elements arc sll — e and whose non-diagonal elements arc zero then if A 
IS any p Xq matrix each of the products (cBp)A and A (eEf) is obtained by multi- 
ply mg each and e\er> element of A by e Wo denote each of these products 
simply by cA 

3 Form the two products 


Anstfer 


(a) (-4, 1,2) 


(a) (-2';), 




(-4.1,2) 
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4. Show tlml, if D is n dinRonRl matrix, DA is obtained by multiplying each row 
vector of A by the corresponding diagonal element of D whereas AD is obtained by 
multiplying e-'icli column vector of A by the corresponding diagonal element of D. 
S'oif. A din^ional matrix /> is a square matrix of which cverj' non-diagonal element 
is zero; thus d/ 0 if I: y- j. 

5, Show that if Di and /Jjarc any two diagonal matrices of the same dimension 
then DiD\ ra DiDi', in words, muHtpUcalion of diagonal matrices is commutative. 

f). Show that if A is any m X n matrix and li is any n X m matrix, so that both 
AB and BA arc square matrices, then the sum of the diagonal elements of AB is 
the same as the sum of the diagonal elements of BA. Hint. The sum of the 
diagonal elemenlH of either .AB or BA is the sum of all products obtained by multi- 
plying any element of A by the associated element of B where we understand that 
the element of B which is associated with Op’ is bg’’ (note the interchange of row 
and eokuun’, for esau^plo, the clcmetU of D which is associated with the element in 
the (iccoiul row and third column of A lies in the third row and second cohinin of 
Bj. Xotc. We term the sum of the diagonal elements of any square matrix the 
trace of the matrix, and we denote the trace of A by Tr ri. The result of the pres- 
ent exercise appears then ns follows: 

TrAB = TrBA 


(althoiigh AB and Be\ arc in general different and, indeed, in general of different 
dimen'-ions). 

7. Verify the result of Exercise 0 for the two matrices of Exercise 3 and for the 

/ 6 ’\ 

two mnlrices A « (oj, nj, oj) and B ~ lb"). What are the dimensions of 

wJ 

AB and BAI Note that AB is the product of the column vector B by the row 
V'-etor A , 

B, Show that matrix multiplication is associative; i.e., if A is any m X n matrix, 
B any n X p matrix, and C any p X q matrix then both the products {AB)C and 
A{BC) may be forme<l and they turn out to be the same. Hint. The element in 

n p 

the rth row and xth column of either product is the double sum n of hj' C,’. 

t " 1 J •* 1 

Aefr. \\r dfnotf the common product simply by AFiC. 

9. bhow that matrix multiplication is distributive with respect to addition, 
i.e., that {A 4- B)('. >» AC -f- BC. Xotc. Addition of matrices is defined only 
(or mnlriees of the same type, i.e., of the same niimber of rows and columns, and 
the rule of n'hiitinn may be pbra“wl a.s follows: Respect the row and eolumn position. 
Thus if ,l and B are e.aeh rn X „ matriet'S the clement in the pth row and 9th col- 
umn of A 4- B is (by definition) 

Of’’ 4* bgV_ 

10. bhow th.at addition of matrict"- is n commutative operation (what doc« this 
mean?). 

11. J-luiw that A -f .t •» 2,1, 

12. Show th.»!,t Tr (.4 4- B) - Tr .4 4- Tr B. 

^13, >hftw that th.e of .4 4- B is the cum of the transposes of A and B. 

• If A i- any n K n matrix ue term the matrix obt.ained by interchanging the 
rvTu f.sd eahmiv.s of .1 the t'anspore of .4 and Wr denote the tran=ix>se of .4 bv an 



s 


VECTORS AND MATRICES 


attached star A* Thus the present exercise may be written as follows 

sa A* + B* Therow vectors of ^4 * are the column vectors of 4, and the column 

^ ectors of A • arc the row vectors of A 

14 State in words and prove that (A*)* ■= A 

15 Form the cofactor matrix of each of the matrices ^ Note 

If 1 13 any square matrix of dimension n we understand by the cofactor matrix of A 
(denoted by co A) the matrix obtained by replacing eacb element of A by its 
cofactoi the cofactor of Op* being tbe product of the determinant of the (n - 1)- 
dimensional matrix obtained by erasmg the ^th row and pth column of A by 

(-! - 1 ) (4 4 ) 

16 If A IS a square matrix of dimensioo n show that A(co A)* = (det A) £« 

where det A denotes the determinant of A flTofe We shall not stop for a defini- 
tion and discussion of the elementary concept of the determinant of a square 
matrix bote, however that deterovinants may be dejined, by induction from 
n - 1 to n by means of the result of this exercise, it being understood that tbe 
determinant of a ooe-dimcosional roatnx is the (one end only) element of the 
matrix Thus to define det A where A is the two-dimensional matrix d^"® 
have, first, co A (see Exercise 16) Hence (co A)* ” 

andso A(eo A)* - ^ ® - (arf - 6c)£i Hencedet A - od - fc 

/ 2 4 5\ 

17. Calculate co A and det A for A » ( 0 -1 — 1 J 
\-3 -2 0/ 

/ 2 d -3\ 

Answer co A *• I 10 15 — 10 J det A »• 1 
\ I 2 -2/ 

18 Show that if A is a square malru of dimension n whose determinant is not 
zero the matrn B ■« (det A)“’ (co A)* has the property that AB — BA Ei, 
Note We term the matrix B defined m this way the reciprocal of A and we denote 
it by the symbol A~^ A square matrix whose determmant is not zero is termed 
non Singular 

19 Show that if AB •• Bnt where A is a square matrix of dimension n, then A 
13 non-singular and fJ = A~^ Jfrnt A is non-eingular since del (AB) — det 
“1 and det ( IB) •• det A det On multiplying the relation AB 
through on the left by A“* and using the associative properly of malnx multiplica- 
tion ne obtain B •• A“* 

20 Show that if BA — where A is a square matrix of dimension n, then A 

IS non-aingular and B A“* Note The results of Exercisea 19 and 20 may be 
briefly expressed as follows Only non^ingular matrices A can satisfy either of the 
relations AB — £» BA •• En, and when A is non-amgular the afeoeiatcd matrix 
B IS unambiguously determinate, being A“* 
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21. ShoTV" thnt if A nnd Ji arc non-pingular matrices of dimension n so also is AB 

and dctcrrame (/IB)"’. J/tnl. => En. 

/Irwircr. (/IB)"’ «=« B~’/l~’. 

.Vofc. The imi>orlant result of this exercise may be phrased ns follows: The 
reciprrieitl of Ihe product of hro non^singular jiwinccs of dimension n is the product of 
their reciprocals in the reverse order. 

22. Generalize (he result of Exercise 21 to the product of any number of non- 
singiilnr matrices of dimension n. 

2-1. Show that if .-1 is any m X n matrix and B any n X p matrix, so that the 
product AB can be formed, then the product B*/l * can be formed and B*A* =* 
(AB)*. 

21. Generalize the result of Exercise 23 to the product of any number of matrices. 
A‘e/r, The important result of this exercise may be phrased ns follows: The trans- 
pore of the prcnhict of amj number of matrices is the prcHiucl of their transposes in the 
reverse order. 

23. .‘ihow that if /I is any square matrix then Tr A* — Tr A. Hint. The 
diagonal elements of .-1* arc the same ns the diagonal elements of /I. 

2. The matrix element of arc and the element of area in tlic 
plnnc 

I.ct tlio roctanRular Cartesian coordinates (x, y) of a variable point P 
in a piano he functions of a single independent variable, or parameter, 
a. Then P traces, as o varies, a curve (which reduces to a point if 
■I nnd y arc each constant functions of or). We term the 2X1 matrix 

the malrix elemoU of arc of the curve that is traced by P. The squared 
tnutmitudo of the vector v{dx, dy) is furnished by the formula 

{ds)‘ — {dx)~ + (dy)" = (d.r)* dx 

whore the 1X2 matrix (dx)* = (dx, dy) is the transpose of the 2X1 
matrix dx. We (enn d-s, i.e., the positive square root of (dx)* dx, the 
scalar element of arc of the curve traced by P. 

Lot, now, the rcctangiilnr Cartesian coordinates (x, y) of P be 
functions of two indei)eiidcut variables, or parameter.^, a nnd /3. We 
have now two 2 X 1 matrix elements of arc rfoX and d^x where 



and wc set up the 2X2 matrix whose column vectors arc doX and 
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The absolute value of the determinant of this matrix is the element oj 
area m the plane, it being understood that this determinant is not 

zero, 1 e , that det 5 ^ 0 Since an interchange of the param 

(“i P) 

eters a and /3 changes the sign of det 7 ^-^ we may choose these 

(«» p) 

parameters m such an order that det p - > 0 and -ne suppose 

this done Denoting the element of area in the plane by dS, we have, 
then, 

dS — det ^ da 
(cc, 0) 


It frequently happens that x and y are not given directly as functions 
of the independent variables a and 0 but rather as functions of two 
other variables, uhich are functions of a and 0 For example, m the 
theory of plane deformations a point Pa whose coordinates relative to 
any convenient rectangular Cartesian reference frame are ( 0 , b) is 
deformed into a point P* whose coordinates relative to any convenient 
rectangular Cartesian reference frame (not necessarily the same as 
before) are (x, y) We have now two elements of area, an tmhal ele 
ment of area dSa defined by 

dSa “ dot ^ da d0 
(«> p) 

and a final element of area dS* defined by 

dSa = det 7 ^^ da d/3 

(« < 3 ) 

It follows immediately from the rule of multiplication of determinants 
(in accordance with which the determinant of the product of tw 0 square 
matrices is the product of the determinants of those matrices) that 

dS* = det^^dSo 

(a, 6 ) 

(prove this Ihnl 7 -^— — ) The easiest and quick- 

\ (a, 0) (o, 6 ) (a 0) / 

est way of evaluating dS, is to proceed according to the following 
formal procedure Wnte dow ndx = Xada + Xb db, dy = j/a do + yb dh 
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aixl nnilfiply d’j I>y 'Oifrc it is understood that the differentials da 
and dh obey the follov.-ing allcrnaling rules of multiplication: 

da da = 0; db da = —dadb; dbdb - 0. 

We find that dx dj/ = - xiyn)da db, and on setting da dh = dSa 

we see that f/.S^ is dx dy. 

Exainjdc 1. Consider the plane shear defined by the formulas 

X = fl d- Fh, 

2/= b 

where h is a constant. Wo have dx = da + b db, dy = db, and so 
diSj = dx dy = da db — dSa- 

'I’liiis the final clement of area dSr is the same as the initial element of 
area dSa. 

INninple 2. Tx't the independent variables (a, /3) be plane polar 
eoordinatc.s (r, 0). From the equations .t = r cos 6, y = r sin e wo 
roiul off dx = cos Odr — r sin 0 dO, dy = sin 0 dr + r cos 0 dO, and on 
evaluating dx dy (where dr and dO obey the alternating rules of multi- 
plication, dr dr - 0, dOdr — —drdO, dO dO = 0) wo obtain dxdy = 
r(co-5' 0 -h .‘iin* 0) dr dO ~ r dr dO. In w'ords, the clement of area in 
plane polar coordinates is r dr dO. 


EXERCISES 

t. Siiow ttial, ill till' iiluiic f-licar defined by jr ■= ri -f hb, y — h, vectors parallel 
to llie n-a\i“ are iinehanned wliereai veotor.s parallel to the 6-axis arc turned 
Ihfoiadi the siicte .are tan I: and arc iiicrea"ed in magnitude by the factor (1 + 
fltnL <\r •» A) <lh, dy == dh. When db = 0, v{dr, dy) = v((la, dh), and, when 
i?a ♦" 0, r(ifr. dy) « r(/; db, dh). 

2. Tleicrmine the effei’t of the plane shear of Exercise 1 on a x'cetor that is 
inrlmtxi .at an tingle rif r/-i to tlie a-.axi-:. 

n. Hliotv tlint in tlif plane tnminn definctl by 

- r cin {(? -b k) r-i rj fo-: f: — 6 sin k 

’ k con“tant, 

h '>■ r .<.111 (t* -b 1) a .‘■in k -r 6 cos k 

the rii!-.! clctui'tit of :uf a dS^ {■< the mme as the initial clement of area dSo. Show 
nho tl.n! the magnitud.- of any vector in the plane is unaffected by the plane 

Viji 


X .Matrix clcim'nts of arc and of surface area and the element of 
xtrhjine in spnee 


I.rt the rci'iniigahir Carte-ian coordinates (x 
in. '=!p,ic'' be furis'tion-- of a .single independent 


y, c) of a variable point P 
variable, or paranictcr. «. 
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Then P traces, as a vanes, a curve (which reduces to a point if x, y, 
and 2 are each constant functions of a) We term the 3X1 matnx 

© /xa dc\ 

~ I da I 

\ 2 a daj 

the matnx element of arc of the curve that is traced by P The 
squared magnitude of the vector v{dx, dy, dz) is furnished by the 
formula 

(ds)’ = (di)* + (dy)* + (d 2 )* = {dxydx 

where the 1X3 matnx (dx)* = {dx, dy, dz) is the transpose of the 
3X1 matnx dx We term ds (i e , the positive square root of (dx) * dx) 
the scalar dement of arc of the curve traced by P 
Let, now, the rectangular Cartesian coordinates (x, y, z) of P be 
functions of two independent vanables, or parameters, a and j 8 We 
havenowtwo3 X I matnx elements of arcd«x and djxwhere 

( Xa da\ / xg d$\ 

Va da ), dgx = \y$d? ), 

Za daJ d^/ 

and we set up the 3 X 2 matnx whose column vectors are dje and 
dgx 

( Ta da Xg d&\ 
y. da yg dp I 
Za da Zg dp/ 

If we regard this 3X2 matrix as obtained by erasing the first column 
of a 3 X 3 matnx the first column of the cofactor matnx of this 3X3 

matrix is f dS*' ) w here dS* = det —— da dp, dS’’ = det 7 -^-^ da dp, 

(«. « («, &) 

dS* « det da dp We term this 3X1 matnx the matnx ele- 
(a, P) 

meni of area of the surface traced out by P The vector r(dS*, dS", dS‘) 
is the vector, or cross, product d,! X d^x, and so v{dS*, dS", dS') has 
the direction of one of the two normal directions at P to the surface 
traced out by P (it being understood that d,x X d^x is not the zero 
vector) In other words, the column vector of the matrix element of 
area is normal to the surface If dS* denotes the matnx element 
of area the scalar clement of area dS, is the positive square root of 
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(dS')" dS^ \vhcre (dS^)*, the transpose of dS"^, is the 1X3 matrix 
(dS', d,% dS^). 

It frcquentlj" happens that x, y, and z are not given directly as 
functions of the independent variables a and d but are given rather as 
functions of three other variables, which are functions of a and d- For 
example, in the theorj' of space deformations a point P„ whose coordi- 
nates relative to any convenient rectangular Cartesian reference frame 
are (a, b, c) is deformed into a point Px whose coordinates relative to 
any convenient rectangular Cartesian reference frame (not necessarily 
the same jts before) are {x, y, z). We now have two matrix elements 
of area; an iniliol matrix element of area, 

I dot da dp 

(«, d) 

j del - -- -- da dp , 

. («, d) 

det^^dadd 
\ (a, d) 

and (i, final matrix clement of area. 


dS' = 


Adopting the alternating rule of multiplication of differentials (in 
accordance with wliich da da = 0, dp da = -da dp, dP dp = 0), we 
ntay write dS"' in the more compact form 

/db dc\ 
dS'^ = I dc dfl I, 

\da db/ 

and VC may write dS' in the more compact form 


^lot 



- da dp ^ 


(«, d) 

del da dp 

(«, d) 

dal da dp 

(«, d) 



dK-' 


V/y dz\ 
dz dz j, 
^dr dy/ 
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Since t/o = j/aOa + j/fcba 4* ytfia, etc , ^vc have 


det 


{y, g) 

(«, /3) ■ 


(6, c) (a, /5) (c, a) (o, 0) 


_i_ ^ f (y» . (°> fc) , 

4- det ; — rr det ; — —> etc , 
(o, b) (a, 0) 


(prove this), and so 

, fa. z) 


dS“ = det iiS“ + det (iS‘ 
(b, e) (c, o) 


+ det ^ dS' 

fa, 6) 


Jacobian matrix • 


This and the t\\o similar equations, which furnish dS^ and dS’‘ (wTite 
dowTi these equations), tell os that the 3X1 matrix dS® is the product 
of the 3X1 matrix dS“ by the 3X3 matrix coJ where J is the 
fa, V, i) 

(o b, c) 

<fS* = (coJ)dS“ 

Taking it for granted that J is non singular (i e , that det J 7^ 0), we 
may \vrite this relation m the equivalent form 
dS* « (dety)(d*)"^dS“ 


Note that det 7 = dety* and that {J*)~ 


Vdet J*/ 


CO J 

Cy. 2) , 


It 


,dat|!i^dS- + de.,^'dS‘ + 

(6, c) fa, o) 


is clear from the relations d5* 

det 7^^ dS‘, etc , that dS* may bo calculated by multiplying out 
(o, 6) 

dy dz = (j/o da 4- yip dfe 4- Ve dc)iza da 4- Z6 d& 4- Zc dc), subjecting the 
differentials da, db, dc to the following alternating rules of multiplica- 
tion do da “ 0 = db db = dede, dedb = —dbdc, dadc — — dedo, 
db do = —do db, and finally setting db dc = dS“, dc do *= dS*, do db = 
dS' 

ExatDple 1. Determine the matnx element of area of the sphere 
a; = r sm cos 0, y = r sm 0 sin 0, z = r cos 6, r constant We have 
dy = r cos sm dfl 4- r sin B cos ^ d0, dz = —r sin 6 dd, and so 
dS® = dydz — r’ sm* B cos tftd^dO — r* sin* B cos ^ dd d<t> Simi- 
larly, (show this) dS^ == dzdx ^ sm* 6 sin <#> dd dtp and dS® = 
r* sm B cos d dd d<ti, and thus the 3X1 matrix element of area is 


dS‘ 


( sm 6 cos 4>\ 
sm 0 sm ^ I 
cos 6 ) 


dd dtt> 
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Hcncfi (sliou- thif!) the scalar element of area dSs — r" sin 6 dO d<f>, and 
(he onliimn vector of the matrix clement of area has the direction of 
the outward dra^^^l radius. 

lAuinple 2. Determine the matrix elements of area of the surface 
into which the surface of a rectangular box whose faces arc parallel to 
(he various coordinate planes is deformed bj' the .shear x = a + kc, 
ij ^ li, z = e, it being understood that J: is a constant and that the 
point Palin, b, c) is deformed into the point PtI{x, y, z). 

We have dx — da + k dc, dy = db, dz - dc, and so 

(dydz\ /dbdc \ /dS^ \ 

dS~ = ( dz j = ( dc dn ) " ( )• 

\d.T dyj \da db + k dc db/ \dS' - k dSy 


Tims the matrix clement of area of a face of the box which is parallel 
to the a-eonrdinate plane (for which dS*' = 0 — dS"^) is changed by the 


shear from 



to dS" 



The scalar clement of area is 


magnified by the. factor (1 + The matrix clement of area of a 

faci* of the box which is parallel to the 5-coordinatc plane (for which 
tiS'' >" n d.S'’) is unchanged by (he shear (prove this), and, similarly, 
the matrix element of area of a face of the box which is parallel to the 
c-coordinate jilunc (for which d.S** = 0 = dS^) is unchanged by the 
shear (prove this). 

Let. now, (he rectangular Carte.sian coordinates (x, y, z) of P be 
functions of three indi'pcndent variables, or parameters, a, /3, and 7. 
We have now thrt'c 3 X 1 matrix elements of arc d„T, d^x, and d-^ 
where 


/x„dct\ /r^ddX /xydy\ 

d„x - I i/„ da j. dfX = { Ilf d|3 j, dyX = I 1/7 dy j, 
\z„ da/ \zf dS/ \zy dy/ 


and we i.et up (he 3X3 matrix whose column vectors are daX, dsx, 
and dai 

( X„ dn Xf dS Xy dy\ 
f/„ da y. dS t/y dy j. 
z„ do Zf dd Zy dy/ 


Tie- 

Xi'Ui* 


al^ulute value of the determinant of this matrix is the element of 
in space, it being undcr.-totKl that this dctcrmin.ant is not zero, 
. , Cx. V. r) 

ttini lift 0. Ve can onier the independent variables 
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so that det > 0 (show this), thus, on denoting the element of 

(a, 0, y) 

volume by dF*, we have 

dV X = det a* da df 

(«i y) 


It irequent\y happens that i, and z are not given diiecdy asinnctions 
of the independent variables a, p, and y but rather as functions of 
three other variables, which are functions of a, P, and y For example, 
m the theory of the deformation of a three-dimensional medium, a 
point Pa whose coordinates relative to any convenient rectangular 
Cartesian reference frame are (n, 6, c) is deformed into a point P, 
whose coordinates relative to any convenient rectangular Cartesian 
reference frame (not necessarily the same as before) are (i, y, z) 
We have pow two elements of volume an imhoZ element of volume 
dVct defined by 

dVx “ det da d/3 dy, 

(a> 0, y) 


and a final element of volume dVx, defined by 

dVx — det da d/3 dy 

(a, 0, y) 

(it being understood that the independent vanables a, 0, y have been 
so ordered that dFa and dF,, defined in this way, are positive) It 
follows immediately from the rule of multiplication of determinants 
that 

dF, = det dF. - (det J) dV. 

(o, b, c) 


where J is the Jacobian matrix of (jc, y, 2 ) with respect to (a, b, c) 
(prove this) An easy way to evaluate dF* is to multiply out 
dxdydz = (xada + zj dZ> + Xedc){yada + y^d& + ycdc){zada + 
Zb db + Ze dc) (it being understood that the differentials da, db, dc obey 
the alternating rules of multiphcation da da da ~ 0 = da do db — 
etc, do dc db = —da db dc, etc ) and to set do db dc = dFo If Pa and 
Px denote the imtial and final densities of the medium that is being 
deformed the principle of conservaiton of mass tells us that 
Px dVx = Pa dVa, 

and so the relation dF* — (det J) dF« may be ivntten m the equiva- 
lent form 

pa = Px det J 
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EXERCISES 

1. Fhon- that the final 3 X 1 matrix element of area dS^ is connected with the 
iniiinl 3 X 1 matrix element of urea rfS* by the formula 

(is- = - (.;*)-* d.S“. 

Ps 

2. Fliow that voliimc-j arc jireservcd under the shear a- = a + t'c, ij == b, z ^ c. 
Uinl. dr dn dz da db dc. 

.1. Fhow that volumes arc preserved under the torsion r = r cos ((? + be) = 
„ Ic ~ I, sin /.c, y «= r sin (0 + he) «= a sin kc + b cos he, t ’= c, k constant. 

t. Show tluit the element of volume in cylindrical coordinates (r, 0, z), where 
f tK T co‘< 0, y ^ r sin 0, z >=> r is r dr d 0 dr. 

3. Show that the element of volume in space polar coordinates (r, 0, 4>), where 
f -m r sin !? cos <iS, !/ ” r sin 0 sin <{>, z ^ r cos 0 is r* sin 0 dr d 0 d <}>. 


t. The cnnonicnl form of a linear operator 

'fhe disctission throughout the present section will concern a square 
tnalrix of dimension 3 but the argument is applicable to square matrices 
of tiny dimension (repeat the discussion in full detail for a square matrix 

(u, V, w) 

of dimension 2). Wo have seen that the Jacobian matrix J — - 

(x, y, s) 

of three dlffcrentiahlc functions (w, v, w) of (x, y, z) may be regarded as 
!v magnification factor that converts, by multiplication, the 3X1 
matrix dx into the 3X1 matrix du: 

du — J dx. 


Himilarly, any 3X3 matrix M may be regarded as a magnification 
factor that converts, by multiplication, a 3 X 1 matrix a into a 3 X 1 
matrix 6: 

h — Ma. 

li.s now consider the .space vectors a and b whose coordinates 
relative to any convenient rectangular Cartesian reference frame are 
ihc elements of o and 5, respectively: 

a = r(o^ a", a^), b = Ir, 6®). 

l.f't the coordinates of a and b relative to a second rectangular Car- 
terian reference frame (which may he obtained from the first coordinate 
reference frame by a rigid rotation) be furnished by the 3 X 1 matrices 
id and h\ re.=:nectivcly. (A'ofe. A prime is often used to indicate the 
tranepoH' of a matrix but we have adopted a star to indicate the tran'j- 
pfKf* jvfjfj lejjves the prime free for the present purpose; a' is 7 iol 
the tran^pav'C of o; in fact, n' is, like o, a 3 X 1 malri.x, whereas the 
s,.'ui«po-v< fi* of o a 1 X 3 matrix.) T.et the coordinates of unit 
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vectors along the axes of the new reference frame relative to the original 
reference frame be furnished by the three 3X1 matrices Ui, uz, and 
143, and denote bj R the 3X3 matnx whose column vectors are «i, 
uz, and Uj The fact that the three unit vectors along the axes of 
the nev. reference frame are not only of unit magnitude but are also 
mutually perpendicular is conveniently expressed by the matrix 
equation 

R*R = Ei 

(prove the validity of this equation) The fact that det R* = det R 
assures us (why’) that det R = ±1, and the fact that the second 
reference frame may be obtained from the first by a ngid rotation 
restricts us to the relation det R = 1 (why’) (Hint det R is a 
continuous function of the elements of R, hence it can assume only 
one or the other of the values ± 1 Since it is originally, when R = 
Ei, +1, It must be permanently +1 ) Thus R is a 3 X 3 matrix 
rvhose determinant is 4-1 and uhose reciprocal is its transpose We 
term any such 3X3 matrix a three dimensional rotation matrix If 

c — I c* I is any 3X1 matnx we have Rc = c^wi + c^ut + c^us 

vv 

This equation tells us that the coordinates m the onginal reference 
frame of the vector, hose coordinates in the new reference frame are 
furmshed by the elements of c, are provided by the elements of Rc 
On setting c =■ R“o and using the fact that RR* = Rs, ''e have the 
following result 

The coordinates m the onginal reference frame of the vector, whose 
coordinates m the new reference frame are furmshed by the element of 
R*a, are provided by the elements of a In symbols 

of = R*a 

or, equivalently (why?), 

a = Ra! 

Similarly, b = Rb', and so the relation b =« Mo yields Rb’ =* MRo' 
or, equivalently, 

V = M'a’ 

where M' = R*MR We say that the elements of M furnish the 
coordinates in the onginal reference frame of a linear vector operator 
M and that the coordinates of this same linear vector operator in the 
new reference frame are furnished by the elements of M' = R*MR 
(just as the coordinates of the vector a in the onginal reference frame 
are furnished by the elements of the 3X1 matnx a nhercas the coor- 
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and b = d( 6^, b®, b^) are said to be orthogonal (or perpendicular) when 
b*a = 0 (Show that the orthogonality of two complex vectors is a 
commutative relation, le , that b*o = 0 if, and only if, a*b = 0 
Hint b*a = (o*b) * ) A 3 X 3 complex matnx each of whose column 
vectors is of unit magnitude, these unit vectors being mutually per- 
pendicular, IS termed a vntlary matnx, and we denote a unitary matnx 
by the symbol U 


EXERCISES 


1 Show that a 3 X 3 matrix U la unitary if, and only if, V*U = Et Note 
The result of this exercise tells us that a 3 X 3 matnx is unitary if and only if, its 
reciprocal is its star, i c , the Go&jiigate of its transpose 

2 Show that 8 3X3 matnx U is unitary if, and only if, UU* = Es 

3 Show that eiery rotation matrix R is umtary 

4. Show that if V is unitary det V ts & complex number of unit modulus i e , 
that det U * 0 real Htnl det U* is the complex conjugate of det V 

5 Show that if each of the column vectors of a unitary matrix is multiplied by a 
complex number of unit modulus (not necessarily the same for the various column 
vectoia) the resultmg matnx is still uiutary 

6 Show how to construct a unitary matrix whose first column vector is an 
arbitrary given unit vector Htnt The second column vector ut is any unit 
vector perpendicular to the given unit vector ui, and the third column vector is a 
unit vector whose coordinates are proportional to the conjugates of the coordinates 
of (Oi X uj) 

7. Show how to construct a ummodutar uiutary matnx, i e , a uiutary matrix 
whose determinanl u unity, whose first column vector is an arbitrary given unit 
vector fftnl Combine Exercises 6 and 5 

8 Show that the general two-dimensional unitary matrix is of the form 



where $ is real and fia -(• ^6 ” 1 

9 Show that the general ummodular two-diraeoaional unitary matrix 
form 


it 1 ) 


where Oa + 6b = I 


of the 


10 Show that the product of two nmtaiy matrices is unitary Htrii 


{ViUsVUiU, = U,*Ui*UiUi - U2*Ut = Et 

11. Verify that the product of two matnces of the type where Sa + 

B 1, IS again of this type 

Ummodular umtarj’ matnces U play m the complex field the role 
played by rotation matnces R m the real field Thus, if a is any 3X1 
complex matnx whose elements furnish the coordinates of a three- 
dimensional vector a, the coordinates of this same vector jn a second 
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complex rectangular Cartesian reference frame are furnished by the 
matrix a' = U*a where the column vectors of U furnish the coordi- 
nates, relative to the first coordinate reference frame, of unit vectors 
along the axes of the second coordinate reference frame. If M is a 
linear vector operator that converts, by multiplication, the complex 
vector a into a complex vector b, 

b = Ma, 

wo have the two equations 

h = Ma, b' = M'a’ 

where the elements of the two 3X3 matrices M, M' furnish the 
coordinates of the linear vector operator M in the first and second com- 
plex rectangular Cartesian reference frames, rcspectivch^ The con- 
nection between the two different presentations of the same linear 
vector operator M is furnished by either of the two equivalent formulas 

M - U*M'U, M' = UMU* 


(prove the equivalence of those two formulas). 

Wc arc now prepared to attack (he problem of determining a particu- 
lar rectangular Cartesian reference frame (in general, complex) in 
which the presentation il/' of a linear vector operator M (whose pre- 
sentation in a given rectangular Cartesian reference frame, usually 
real, i.s a given 3X3 matrix M) is particularh' simple. Let X be a 
rani of the cubic equation det (ilf — XE 3 ) = 0. Then there exists a 
non-?,oro 3X1 matrix a lliat .satisfies the equation Ma = Xa (why?). 
If X happens to be real and if the elements of M are real the elements of 
a may be taken to be real, but if X is a non-real complex number the 
elements of a will not all be real when the elements of Af arc real (for 
if the elements of a were real the elements of Ma would be real and we 
could not have Ma = Xo since X is, by hypothesis, non-real). In any 
event let uj bo a 3 X 1 matrix, whose column vector ii\ is of unit 
magnitude, which is proportional to a (i.e., «i = f:a whore k is some 
complex number). Con.slnict a 3 X 3 unimodular unitar}- matrix 
f- ( wluKc first column vector is «i. Then, on denoting bv ei the 3X1 
/ 1 \ 


matrix 


0 

A 


I. we have the two relation.'^ 


Mill *= X«i, tti -- UiCi 


(sluMv tliis), and these imply the relation MUiCi = XD'iCi or, equiva- 
lently. 


Xf,. 
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In other ^\ords, the first column of Ui*MXJi is ( 0 1 Thus the 

\o/ 

comple\ rectangular Cartesian reference frame determined by Ui js 
such, that the presentation, in it <if the linear vector operator M has 
zeros m the second and third rons of the first column We may nrite 
this presentation of M m the form 


r-') 


where is a tno-dimensional matrix and the stars denote numbers 
(complex) about iihose lalues ne can say nothing further Applying 
to N the argument just presented concerning M, ive construct a tno- 
dimensional ummodular umtary matrix V which is such that 


The matrix 172 * 


/I 0 0\ 

— I 0 ^ ) IS a three>dimensic 

\0 / 

umtary matrix (prove this) which is such that 

i it follows that 

U2*Ui*MViUz = (o *] 

\0 0 y/ 


three>dimensional ummodular 


(prove this), and it follows that 


U2*Ui*MViUz = 


The matrix U 1 U 2 is a three-dimensional ummodular unitary matrix 
that we denote simply by U We have, then, the following funda 
mentallj important result 

If IS any three-dimensional square matrix, there exists a three 
dimensional ummodular unitary matrix U which is such that U*MU 
has zeros below the diagonal We term the matnx M' = U*MU, 
which has zeros below the diagonal, the canonical form of the linear 
vector operator whose coordinates, in the original reference frame, are 
SwTmshed hy the elements of If 
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If .R is ssicii ihn( .R* is the Fame as *R we term it symmetric if its 
('lcnu‘nt« are real and llcnnilinn if its elements are complex. It is 
clear that df' = U*MV satisfie.s the same relation (.1/')* = dP 
(prove llii.y) and so M' is diiigonal with real diagonal elements. If 
tlie elements of M are real it follow.s (pnn'o this) that the elements of 
r may he taken to he real (.so that U i.s a rotation matrix). We have, 
then, tlie following im]mrtant special case of our principal result; 

If .1/ is atu' threo-dimen.sional .symmetric matri.x there exists a rota- 
tion matrix It .such that It*MIt i.s diagonal with real diagonal elements: 

fmi 0 0 \ 

R*.VR = ( 0 m-2 0 • 

\0 0 mj 

The three real numbers (mi, uin, ws) are known as the characteristic 
numbers of the linear vector operator M whose coordinates in the 
original reference frame tire furnislied by the elements of the real S3'm- 
metric matrix il/. 
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real complex number x*Mx is never negative (no matter what is the 
complex vector x) we say that the Imear vector operator M, whose 
coordinates are furmshed by the elements of ilf, is non-negative It is 
clear that M is non negative if, and only if, its charactenstic numbers 
are non-negative (Prove this lltnl x'*M'x' = -f 

niix'^'x^' + mzx^ x ^' ) We non construct the 3X3 matrix 



where denotes the unambiguously determinate non negative 
square root of mi, etc , and it is clear that the square of this matrix is 
the canonical form M' of M We denote it, therefore, by 
and it is clear that {M')^ commutes with M', i e , that = 

Af'(Ar')^*, since both M' and (3/')^ diagonal matnces On setting 

M” = 

(so that the elements of 3/^ are the coordinates in our ongmal refer- 
ence frame of the linear vector operator whose coordinates, m the 
reference frame m which the coordinates of M are furmshed by the 
elements of 3f', are furmshed by the elements of {M )^), it is clear 
that 18 a non negative linear vector operator (why?) and that the 
square of is M (prove this) It is also clear that commutes 
with M, 1 e , that ** MM^ (prove this) 

If A 18 any 3X8 matnx it is clear that the 3X3 matrix M defined 
by M — A* A furmshes the coordinates of a non negative linear vector 
operator In fact, M* = M (prove this) and x*Mx = x*A*Ax = 
(Ax)* Ax IS never negative, being the squared magmtude of the vector 
whose coordinates are furmshed by the elements of the 3X1 matnx 
Ax If the matnx A is non singular, M is not only non-negative but 
also positive, i e , x*Mx is positive save when x is the zero 3X1 
matrix (m which event x*Mx is evidently zero) In this case M 
IS non singular (why?), and this implies that 3/^ is non singular 
(why?) Denoting the reciprocal of 3f^ by 3/“^, we can easily see 
that the matrix AM~^ is unitary In fact, (AM~^)*(AM~^) = 
M-^A*AM-^ = = (3f“^‘)*3/ = 3f-^3/ = Et On 

denoting the unitary matrix A3/“^* by U, we have A = UM^ We 
obtain, m this way, the foUowmg fundamental result 

Every non singular 3X3 matnx A whose elements are complex 
numbers may be written in the form 


A = UM^ 
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-ivhcrc Uisti unitan' 3X3 matrix and iV'‘ is a positive (and, hence, 
Ilermitian) 3X3 matrix. 

If the dements of A arc real. M = A* A may bo transformed into its 
canonical form by a rotation matrix (why?), and so the elements of 
.If"'" are real (why?). Hence the elements of U = Ail/"'* are real. 
If def A > 0 it follow.? (.since det > 0) that dot C/ > 0 and so 
f; is a rotation matrix. AVc have, then, the following particular 
in.stancc of our prineip.al result: 

If A is a non-singular 3 X 3 matrix, with real elements, whose 
determinant is positive, A may be written in the form 

A = 

where li i.s a rotation matrix and M = A*A; is a symmetric 
matrix whose characteristic numbens 7ni^‘, arc all positive 

(aij, vit, f?i 3 being flic positive characteristic numbers of M = A*--l). 
Thi.s re.sult is of fundamental importance in the theory of elasticit 3 \ 

6, Conelusion 

If one has miustcred fairly well the contents of thi.s introductorj' 
ehapler he is in po.s.session of all the matrix theory' required in this 
treatment of the finite deformatioms of an elastic solid. For an intro- 
ductory ftud.y at a level .suited for senior high school or beginning 
college 8tudcnt.s the reader is referred to the first five chapters of 
Amhjlk Gromc/ry,* where analytic geometry is treated from the verj’^ 
heginning from (he point of view of vectors and matrices. More 
advanced matrix Iheorj- is presented in the first two chaptera of 
Inimludion io Aji/di'cd MalhcviaUcx,- .and a treatment of those parts 
of the Ihcorj' of matrices which are of importance in group thcorj' is 
given in Theory of Group RcprcsnUnlionsA Try to become Raniiliar 
with matrices and to use thorn regularly as a helpful tool. Do not 
fail into the error of regarding them as a complicated device invented 
hy mathematicians to make the theory of elasticity harder than it 
.neJitaliy i.s. 
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2. Write down the most general 3X3 matrix that commutes with the diagonal 
matrix 


where di di 
Answer A = 


/di 0 0\ 
Z) - { 0 di 0 ) 

\ 0 0 dj 


we can say nothing 


C 0° 3 


where the stars indicate elements about whose values 


{x, y, 2 ) 

3 Calculate the matrix Af = where AulheJacobianmatrix- — ; — of the 

(a, b, c) 


shear x ^ a + ke, y ** b g ^ e (k constant) 

4 Hepeat Exercise 3 for the torsion x = r cos (fl + fcc) « o cos ic — 6 am kc, 
y = r sin (fl + kc) = o sin fee + 6 cos Ac, 2 = e (k constant) 

5 Determine the axis and angle of the rotation whose matrix is 



Answer The line from (0, 0, 0) to (1, 1, 1), ^ 



2 


THE SPECIFICATION OF STRAIN 


1. Ilir slrnin nintrix 


Wc consider tAvo states of n three-dimensional deformable medium: 

(1) The imlial or itnsiraincd state. Wc denote the rectangular 
Cartc-sian coordinates of a typical particle of the medium, when in 
the unstrained state, noth rc.spcct to any convenient reference frame, 
by (a, h, c). 

(2) The ftna} or s^lraitia} or deformed slate. We denote the rectangu- 
lar Cartesian coordinates of the particle, whose coordinates when the 
medium was in the initial unstrained state were (a, b, c), by (x, y, z). 
The reference frame with respect to which (x, y, z) arc measured need 
not be the same (although it is usually convenient to have it the same) 
as the refemnee frame with rc.spcct to which (o, b, c) are measured. 

If wo consider a collection of particles occupying, when the medium is 
unstrained, a certain volume a, b, and c will each be functions of 


throe independent variables (a, d, y) which are such that del 


(a, b, c) 


(a, i3, y) 

0. Wo suppose the independent variables, or parameters, «, 0, and 

> .‘■0 order'd that dot > 0, and wo denote bv dVa the product 

(ff, d. y) 


dot 


(o. b, el 

(rr, B, -j) 


dn df? di. 


dV'a i.*? the imlial element of vohime. (i.c., the 


dement of volume of the medium when it is in its unstrained state), 
and we may, by adopting (he alteniating rule of multiplication (what is 
lhi^?) of the dilTenmtials da, dB, dy. write it in the form da db dc: 


dl'n = dadbdc ^ dot 


(fl. f>. c) 

(«, d. 7) 


da dB dy. 


} he coordinates (x, y, :) of the final po.sition of the particle whose 
initial coorvliisates are (a, b, c) arc a.^^sumed to be differentiable func- 
tion*? of (a, h, r), ami it h .assumed further th.at the determinant of the 

27 
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Jacobian matrix ^ f ' \ is positive Then the final element of volume 
(a, b, c) 

18 given by the formula 

dV, - dx dy dz = J) 


(a, b, e) 


where J denotes the Jacobian matrix 


jx, y, z) 


If Pa and Pi denote, 


(a, h, c) 

respectively, the initial and final mass densities, the principle of con- 
servation of mass assures us that 

P*dF. = PadFa 


and so the compression ratio — is furnished by the formula 


— » ^ 
p« det J 

Let us now consider a collection of particles that be initially on a 
curve Ca Then a, b, and c are functions of a single independent 
variable or parameter a, and the matnx element of arc of Cg is the 
3X1 matnx 



The scalar element of arc dsa of C® is the positive square root of (do) *da 
(dsa)* = (da)* da 

The particles of the deformable medium that onginally lay on the 
curve Ca he, after the deformation, on a curve (7*, and the matnx 
element of arc of C* is the 3X1 matnx 

/da;\ 

dx = i di/ j = J da 
\dz/ 

The scalar element of arc ds* of (7* is the positive square root of 
(dx)* dx 

(dsx)® = (dir)*d* « (da)*J*Jda 


If dsj is the same as dSa for every curve C« (i e , for an arbitrary 3X1 
matnx do) all curve lengths are unchanged by the “deformation” so 
that, in particular straight lines ate “deformed” into straight lines 
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fv.-hy?). Tliifl moans that x, y, and r are linear functions of a, b, and 
c (uhy?) and ho J is a constant matrix. Since this constant matrix 
sritisrscs the relation 

i(Ia)*J*J(!a = id a)* da 

for arbitrary choices of the 3 X 1 matrix dn we must have J*J = 
E% (prov<‘ this). Since det .7 > 0 it follows that J is a rotation matrix. 
^'hu‘? the transfonnalion from the initial coordinates (a, b, c) of an^' 
p.'irticle of the medium to the final coordinates (x, ij, z) of the same 
particle is merely that involved in the tran.-^formation of coordinates 
from one rectangular Cartesian reference frame to another. In other 
words, the change from the initial to the final position of the medium is 
a mere rujid displaccmad and no deformation or change of shape of 
the medium (or of any jiortion of the medium) is involved. When 
./ i.s not a rotation matrix we use the difference between J*J and Es as 
a mca.cure of the deformation of the medium in the neighborhood of 
the point (o, b, c). We .«-ct 


i(./M - /i3) = r, 

and wc s.ay that the 3 X 3 symmetric matrix y measures the deforma- 
tion of the elastic medium at the point (a, b, c). This matrix n is 
hnown as the strain matrix, and it satisfies (ho relation 

(d.v)' ~ (f/s„)* = 2ida)*r]da 

(prove thi.s). The u.sc of the factor i in the definition of y or, equiva- 
lently, of the faetor 2 in the relation just written, is due to the fact 
that this relation involves the difTcrenco of the squares of dsi and ds„. 
If the strain or deformation i.s .small enough so that wc can replace 
dt, hy '2ds„, (dsj)' — (r/s„)' hceomes 2{ds^ — dsa)idsa) (why 

cminoi wc replace ih^ — ds„ by zero?), and we. have the approximate 
relation 



The column vector of the 3 X 1 matrix ~ is a unit vector; if this 
^ dSa 

3 X 1 nmtrix is denoted by u, the quadratic form »*jjh in the elements 


of u measures, approximately, the retntivr rnaimificntions 

dSn 

m the neighborliood of the point (o, fi, r). It is clear from the relation 
(rfif)' - (aVt* ‘2{da)‘r,da that the quadratic form u*riu measures 

. 1 » idrfr 

fN.'U’tly ^ (prove this). 
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The most direct way of calculating for any given deformation is 
first to set up .1 = f and then to calculate the symmetno matrix 

Oj C/ 

J*J The diagonal elements of J*J are the squared magnitudes of 
the column sectors of J, and each non-diagonal element of J*J is a 
scalar product of ti\o of these column vectors After J*J has been 
evaluated, subtract Ea from it, i e , diminish each of the diagonal ele- 
ments of J*J bj unity one half the resulting matrix is the desired 
strain matnx 

Example Calculate tj for the shear x = a + kc, y — b, s = c, 
where k is & constant Here 


and so 


/I 0 A /I 0 k \ 

J*|o 1 o1j*J»Io 1 0 I 

\o 0 1/ Vi 0 1 -h IV 



If k 18 80 small that its square may be neglected we say that the shear is 
infintlenmal thus the strain matnx for an infinitesimal shear parallel 
to the a axis m the <zc-plane equals 



EXERCISES 

1 Calculate the strain matnx v for the uoiforiD compression or dilatation z » 
a-f-Ad i/»*b+A6 *“c-l- ic k a constant Note This is a compression if 
k <0 and a dilatation if I > 0 

Anawer ij = (A +71*)^* 

2 Uliat IS the strain matrix for an loliaiteaimal uniform compression or dilata 
tion? 

Anatcer ij » kFt 

3 Calculate the strain matrix ij for the torsion z — r cos (0 + Ac) — a cos ke — 
6 Bln Ac V — r sin (9 + Ac) ■•a Bin ke + b cos kc z — e A a constant What 
does e reduce to when the torsion is mliiiitc'iinalf 


( ° 

0 

-kh \ 

/ 0 

0 

-kb\ 

Answer 7 1 0 

0 

Ao ) i 

r{ 0 

0 

ka I 

V-Ai 

Aa 

kHu’ + 6')/ 

\-lb 

ka 

0 / 


4 Calculate the strain matnx 7 for the simple tension and cross-scetional con 
traction of a rod defined by z “ a — oka y «■ 6 — okb * c + Ac where A and a 
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nre con'i^nt*. Note. The constant is kno\rn tis Poisson's ratio (after S. D. 
Poi‘ 5 on jl781-lS'10], a French mathematician). 

/-^f : 0 0 \ 

Autv:er. >j •” { ® — tr/; + - 5 ^*/:* 0 )• 

\ 0 0 k+ik"-/ 

.*>. What docs thi' strain matrix of Exercise 4 reduce to when the strain is infini- 
tfiimnl, i.c., wlicn tlie square of k may be neglected? 

/-<r/; 0 0 \ 

itnifiw,'. { 0 —ok 0 j- 

\ 0 Ok/ 

2. ^^lc liclinvior of the strain matrix under transformations of 
the initial and final coordinate reference frames 

A rotation of the initial coordinate reference frame, i.e., of the rectangu- 
lar Cartesian reference frame with respect to which the coordinates of a 
typical particle of the clastic medium are (a, h, c), is described by the 
formula 

o ~» a' = R*a 


where the coordinatc.s with respect to the old reference frame of unit 
vectors along the axes of the new reference frame are furnished by the 
column vectors of /?. We understand by the notation o— > o' = R*a 
that the new initial coordinates of the particle, whose original initial 
coordinale.s were furnished by the elements of the 3 X 1 matrix a, 
are furnished by the elements of the 3X1 matrix R*a. It follows 

(o, h, c) 


from the relation o' = R*a that n — Ra' and, hence, that 


H. Since 


(t, I/, c) (x, y, z) (a, b, c) 


(o', 6', c') 


wc have the following result: 


(o',6',c') (a, 6, c) (o', c') 

Under the tran.'-formation o— * 0 ' = R*a of the initial rectangular 
Carte.rian reference frame, .1 —> J' — JR, the final coordinate reference 
frame being unchanged. 


EXERCISES 

1. Show th.nt undor the tran^formiifion x— ♦ x' <= li'x of the fimil rcctanguler 
CSrie^ion fMtno J —* J’ *J. 

2. Shoo thnt under the simultaneous fransformntion.s a~> a’ «=< R*n, i— ♦ x' = 
h’*xot both the initial niul the final n'Ctangular C.artcsinn reference frames J — * 
J' « R'JtL 


Since M w V and since the transformation J J’ = JR implies 
tlie tran.«formntion J* -+ J*' = R*J* (why?), it follows that M -* M' 
H*MR. Since the strain matrix tj = i(.U - E^), this fact yields 
the following important result: 
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Under the transformation a—*a' = R*a of the initial rectangular 
Cartesian reference frame ij — » V *“ 12*1)5 
In other words, the strain matrix ij is sensitive to the initial rectangu- 
lar Cartesian reference frame There exists (at each point of the med 
mm) at least one particular initial reference frame in which ij appears 
in canonical ( = diagonal) form We term the axes of any such par- 
ticular initial reference frame principal axes of strain and «e write the 
strain matrix i? when referred to a reference frame whose axes are 
principal axes of strain, as follows 



Note carefully that, unless n is a constant matnx, the principal axes 
of strain will vary from point to point of the medium Note further 
that if 5 IS a rotation matrix that transforms n into its canonical 

Ai 0 0\ Ai 0 0\ 

form, so that R*vR “ I 0 i)j 0 J we have »)5 =* 5 j 0 ve 0 1 

\0 0 1)3/ \0 0 i]s/ 

Hence the column vectors of R are characteristic vectors of ?) (why’) 
The easiest way to determine R is to determine the characteristic 
vectors of i) It frequently happens that one of these is relatively 
easy to determine (the determination of the others being comparatively 
difficult) The best thing to do in this event is to construct a pre 
limmary rotation matrix Ri whose first column vector has the direction 
of the easily determined characteristic vector of »j Then Ri*tiRi will 
be of the form 



and our problem becomes the comparatively simple one of a plane 
strain whose matrix is indicated by * * 

In contrast with its behavior with respect to a change of the initial 
rectangular Cartesian reference frame the strain matnx ij is insensitive 
to a change of the final rectangular Cartesian reference frame In 
fact, under the transformation x— »x' = R*x, J—* J' == R*J, J*—* 
J*' = J*R, M -* M' = J*RR*J = M Hence = v (why?) 

Thus the elements of the strain matnx ij are invariants under trans 
formation of the final rectangular Cartesian reference frame, to 
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omplinHZO t1n“ point wc denote these elements by the symbols 

^iiai Vn'fi ClC.. 


f 

Van 

Vcb 

Vac\ 

Via 

Vhh 

ntf 1 

\Vta 

Vcb 

VeeJ 


Note flml thf first of tlie I wo subscripts tells the row and the second the 
(‘ohimn; linis 55?.^ is the element in the second row and third column of 17. 

llnnoTh L We use nah, etc., instead of 170", -n", etc., since, n is 
not the presentation of a linear vector operator but i-ather tlie matrix 
of rocflicients of a cpiadratic form «*J7?r where w is the 3 X 1 matrix 

element of arc ~ of Ca whose column vector is of unit magnitude. 

The fact that the .strain matrix 7? is symmetric is c.xpressed hy the 
relations 

"^ch " Vbc* Van ” Vcfit Vhft Vah* 

RrmnrI: j?. The inscmitivcncss of 77 to a change in the final rectangu- 
lar Earte.sian reference frame may be oxprc.sscd by the statement that 
the .strain matrix i.s unaffected if the given deformation \s followed by a 
rigid rotation of the medium; thus both ,7 and /?*./ yield the same 77. 
'I'he eciiailivcncss of 77 to a change in the initial rectangular Carte.sian 
reference frame may be exprcs.ccd by the statement that the strain 
matrix i.s affected if the given deformation is preceded bj” a rigid rota- 
tion of the medium. Tlius .7 and .777 yield dilTerent strain matrices. 


EXERCISE 

.t. Show that tlip .strain matrix 13 tr.ansfonns under the simultaneous trnnsformn- 
linji n~- o' h‘*rt, x-“* 1“ Il'jr of the initial and final rcctanpulnr Cnrtesi.an 

refi'ft !tre franu". accordinp to the formula ij — ♦ ij' *= Ilinl. 13 is insensitive 

to the trati'forniation x' «• /f*x. 


Exaniple 1: Determine the principal axes of the plane shear x = 
« A /•/», 7/ h, k consl'int. Here .7 = Q j^and.S0 77 = ^2^- 

Setting R r- ^5.0 {],(, initial reference frame 

S'? tibtfiincd by rotating the old initial reference frame through the 
angle n). wc have 


2rJ^ ^ * \ 

\k cos a A k' sin a —k sin « A k~ cos «/ 

^ fk *-511 2 n A k" sin* o k cos 2 o: A i-k~ sin 2 a \ 
\l rri» 2a A kk* .‘•in 2tT ~k sin 2 « A fr cos‘ a / 


2R*v,R 
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In order that be diagonal we must have cot 2a = “oOr, 


2 


equivalently, sm 2a = — — — (o lying m the interval J 2' 

if A; > 0) On eliminating a, we find that ^ 


(4 + kY r 


4(4 + r>)«’ + 1 “ kd 

tesimal theory of plane shear (m nhich I* is neglected) the angle a 


independent of k To this degree of approximation iji = 


K -k 

■ 2”'' “ 2 

The next approximation to «ii and ijt (in uhich we neglect A* and 

The 

-k 


k k} 

higher powers of k but not A.*) m “ - + — , » 


formula cot 2a : 


■ which furnishes the exact angle of rotation is 


so simple that there is no point m approxtmaUng a, if 11 e imte a -* 
7 + /3 we have tan 2d “ r and so the exact angle of rotation may be 


1 


Example 2 Determine the principal axes of the torsion x ■ 
r cos {d + kc) y — r sm (0 + Ac), z ^ c, k constant Here x - 
a cos kc — b 8in kc y = a sin Ac + 6 cos kc s = c, and so 


/cos kc 

— sm Ac 

-ky\ 

/ 0 

0 

I sm Ac 

cosA^ 

kx I 

ij = i( 0 

0 

\ 0 

0 

1 / 

\-kb 

ka 


~kb\ 
ka I 

FrV 


Since det jj = 0, one of the characteristic numbers of rj is 0, and it is 
clear that tj(a, b, 0) is an associated charactenstic vector Setting 
up the preliminary rotation matrix 
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3o 


/O 0 0 

2Ri*nIi\ = ( 0 0 hr 

KO kr 


and we have only to deal with the plane 


fhoar (in the i/'c-plane) whose strain matrix is -§• /-v)' 

if we start with an initial reference frame whose first and second axes 
have the directions of plane polar coordinate lines at the point (a, h) 
in the (a, 6)-plane the principal axes of the torsion arc obtained by 
rotating these axes around the radial direction through the angle 

The canonical form of the torsion strain matrix 

kr 


T 1 kr 

- + - tan j- 


is furnished by the formulas = 0, i ?2 


,.3,3 

-1(4 + k”r‘y^’ 


(4 -f 4 


kr 


(4 -f Jrr") 




+ 


k‘r~ 


k^r^ 


4(4 -f Irr-) 


a- 


If we are 


dealing with a cylinder of length I which is twisted through a quarter 

r r 


turn at c = / we have kl = - and so kr - , , 

4 4 I 


Thus the approxima- 


tion furnished by the infinitesimal theory will not be valid at the sur- 
face of the cylinder (for a torsion of this magnitude) unless the ratio of 
the radius of the cylinder to its length is so small that its square may 
be neglected. 

If we use the some reference frame for the initial and final positions 
of the medium we may introduce the displacement rector v{U, F, TF) 
who,‘<e eoordinatc.s are defined bv the formulas 


I’ *= X — a, 


y 


IF 


It follows that J Fa -f /C where /v is the Jacobian matrix 

nml so .V « ,/V Fa + K* + K + /v‘/v'. Hence 


(U,V, IF) 

(n, b, c) 


V = iPf - F,i) = i(K* + K) -f- iK*K. 

In the infinitrsmal theorj' of strain, K is regarded as an infinitesimal 
matrix, i.e., as a matrix who.«c elements are so small that their squares 
• and products are negligible. To this degree of approximation, then, 

/ Ua 4(F„ 4- U,) A(f', + 1FJ\ 

n - UK* 4- K) - UVa + U,) Ft, A(ir, 4- Fe) j. 

\UKc + F'c) i(ire. 4- F.) ]]\ J 

Fo obtain the e-xnct oxprc.ssions for the elements of r\ wo must add 
\K*h to UK* 4- ?v); remembering that the elemont.s of K.*K are 
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the squared magnitudes and the scalar products of the column vectors 
of K, we obtain 

n.. = V, + + V.' + IIV), 

nu = n + HUi' -h + ir.=), 

1 ., - ir. + + V/- + ir.=), 

li, = 4(iri + r,) + i(u,i/, + ViV, + inirj, 

-I.. - i(t/. + ir.) + + r,v. + TF.irj, 

5.1 - i(r. + t/i) + + w.Wt) 

3. The invariants of tlic strain matrix 

We have seen that the strain matrix ij is sensitive to a rotation of the 
initial rectangular Cartesian reference frame, under the transformation 
o-» a' » R*a, 1} -» ij' = R*vJi It follows that, no matter what the 
number X is, V — \Bi = R*iv — \Ez)R (prove this) Hence det (ij' 
— X£?j) e det (ij - XEj) (why?) If we denote by /j, 7j, and 1 %, 
respective!) , the trace of ij, the trace of co v and the determinant of ij, 
w e have 

det (17 — XEj) “ /j — AX + /jX — X* 

(prove this) If Ii, /j', /j' are, respectively, Tr e', Tr co ij', and det n' 
it follows that, no matter what the numlicr X is, 

/s' - /i'X + /I'X* - X’ = A - /sX + /,X® - X> 

Hence (why?) 

W = h W = h, h’ = /i 

Thus, although 17 is sensitive to a rotation of the initial rectangular 
Cartesian reference frame, the three functions /i, A, and Jj of the 
elements of 17 (or, as we shall simply say, of 17) are insensitive to any 
such rotation We term them the three invariants of rj They are 
furnished by the formulas Ii = i?** + 1766 + »7ce> in words, Ii is the 
sum of the three diagonal elements of 77, 

T libel I |i7ee I7e« i r7oa I7ab 

-12 “ * *»* ' 

|17eb 17ee) |17ae I^Tba llbbl 

in words, /2 is the sura of the three two-rowed diagonal minors of 17, 


la = det 17, 
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in v/ords, /a is the determinant of v- On ^Y^iting in its canonical 
(e diagonal) form, avc see that 

jfj c= Tjj 4 - 132 + nzi h ~ nans + vsnt + nivz, h — niV2V!i- 

If ?3 is nil infinitesimal matrix h, h, nnd h are, respcctivel}', infini- 
icsimals of (at least) the first, second, and third orders. 

It i.« ciisy to express the comprc.«sion ratio — m terms of the invari- 

Pa 

ants /i, /«, and /s of ij. In fact, p„ = Pi (det J) and dot (Es + 217) = 
de( M 5= det ./*,/ = (det J) - and so 

r= del 4- 2i}) = 1 + 2/j + 4/5 + 8/3- 

Hence 


— = (1+2/14- 4/2 + 81 3) 

Pa 

The approximation to (his exact formula which is furnished by the 
infinitesimal theory is 

^ = (1 -h 2/,)-^* = 1 - Ji. 

Pa 


To this degree of apiiroximation, /i = — is the relative decrease 

Pa 

in density or, cquividently (in view of the principle of conservation of 
ma-vs), the relative increase in (local) volume. 

EXJilRCISES 

J. ?Iu)iv that for llic miiform dibtation (or comprcsiion) s «= (1 + /;)a, y ■= 
(1 (1 + constant, 17 «=■ (/; + and c.alculato /i, /j, and /j. 

dnetw. /j « W: + U‘), h " 3 (/: + ilr)-, h ■= (t- + lk-)\ 

2. Verify that, for tlio uniform dilat.ation (or compression) of Exerci.se 1, /; = 
|/i*, fj uV/i’, 

3. Kiioxc Ui.at if the relations of Exercise 2 arc satisfied the strain matrix is a 
j'-'.sl'ir matrix, i.e., a muUisile of Ej. ll\n(. The cubic equation that determines 
the charactcriitic numbers of is (X - ihV ° 0. 

-t, (Siiim- that two deformations having the same strain matrix have their 
in.atrirra J j and J • conneeted by the relation J- => RJi where R i.s a rota- 
tt-m matn.t. //jnf. J, » J. « and m J: - RJi where R « 

■’fjfts*. 

a, Hiew th.at the compit-.s'uon ratio for the uniform dil.at.ation (or compression) 
<■1 i.i.s'r,"*,' 1 furnish»'fi by the formula — «» (1 -f h)~^. 
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6 Show that for the uniform dilatation (or compression) of Exercise 1 the strain 

1 M"” 

matrix 17 is the product of Ei bj 2 \~y ~ ^ Eiit det Af « (I +2r)* 

where ij = cEj 

7 Calculate the invananta of the shear i««o+A:c y = ii 2 = c and deter 
mine the relations between them 

Answer /i «=• ■ji' Is — •Jfc* /* = 0 /i »» — ff/i /| =• 0 

8 Repeat Exercise 7 for the torsion x = r cos (9 + ic) = a cos fcc — 6 sin ilf 
y ••rsnffl+ic) = a 8 nie + fc cos ke t <=> e and verify that this is at any 
point (a b e) a shear whose coefTcient is ir = k{a* + 6*)^^ 

Insiffr Ii = Is = — 4^*^* /a = 0 

4 Tlie compatibility relations between the elements of the 
strain matrix and their derivatives 

Suggestion The mathematical computations involved in thib 
section are quite formidable for a reader who is not well trained m 
tensor analjsis One may safely omit this section m a first reading 
at some later time when the reader has leisure and wants to under 
stand what is behind the compaltbtlUy relations of the mfmitesimal 
theory of elasticity he may return and work through the computations 
of this section 

The basic relation {dx)* dx ida)*M da tells us that the quadratic 
form with coefficient matrix ilf m the elements of the 3 X 1 matrix da 
appears when written as a quadratic form m the elements of the 
3 X 1 matrix dx as a quadratic form with the coefficient matrix £j 
whose elements are constants In the terminology of metric differen 
tial geometry this constancy of the elements of the coefficient matrix 
of the quadratic form (dx)* dx is expressed as follows 

The curvature tensor of the three-dimensional metnc space the 
coordinates of w hose metneal tensor m the (o b c) system of coordi 
nates are furnished by the elements of ilf is the zero tensor 

This fact furnishes six relations involving the elements of M and 
their firsU and second order denvativcs with respect to (o b c) or, 
equn alently since If = Fj -f 217 six relations involving the elements 
of 17 and their first and second-order derivatives with respect to 
(a b c) These arc the (exact) compatibility relations between the 
elements of 17 and their first and second derivatives with respect to 
(o be) If 17 is an infinitesimal matrix and we neglect in the (exact) 
compatibility relations all products of an element of 17 or of anj 
denvatue of an element of n by an element of 17 or any derivative 
of an element of 17 we obtain the compatibility relations of the infini 
tcsimal theorj of elasticity 

A knowledge of tensor analysis wnll not be assumed here and the 
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diKcuhsiun will fsimply fvirnisli that portion (and no more) of tensor 
anfilypis which is necessarj' for the purpose. Following the notation of 
tt'n.‘^or analysis we use (o^, o*, «“) instead of (a, b, c) and (x^, x‘, x^) 
instead of (x, y, r). The double occurrence (once below and once 
above) of a Greek label in a term indicates summation, with respect 
to that label, over the range 1, 2, 3. For example, if x*" is a tj-pical 
clement of a 3 X 1 matrix x and if y. is a typical element of a 1 X 3 
matrix y, the matrix product yx is denoted by the symbol yex" (or, 
equivalently, y^x®, y-yX'', etc.). In this notation the fundamental 
relation (dx)* dx = {da)*M da may be written as 

7Ti„aa'’,ra®., = I, if s = r, 

= 0, if s r. 

Here nip, is a typical clement (the clement in the pth row and gth 
column) of M and the comma in the sj’mbols a“.r, denotes difTcren- 
tiution; a’’,, is the derivative of dP with respect to x*". The occurrence 
of the two Greek letters a and fi in the symbol indicates a 

double summation (one summation on a and one on ^) so that 
is the sum of nine terms (write down five of these terms). We now 
dilTcrcntiate the relation = constant with respect to any 

one of the x’s, say x'. In doing this dilTercntiation we remember that 
the ('lcmont.s of M arc given as functions of the o’s so that the deriva- 
tive of aipj with rc.spect tox' is where ?np,,r denote.s the deriva- 

tive of nip, with respect to n’’. We obtain 

wliere for example, denotes the second derivative of with 
respect to x’’ and z‘ (so that — o^,<r [why?]). We solve this 
.\v.stc-m of equations for the second derivatives as follows: First 
iiitorchange r and I and rewrite the triple summation in 

the equivalent form T7i,,i.„a".ra‘’.,oLf (^'hy arc these two triple summa- 
tion.t the same?). We obtain 

-f „ = 0. 

Interchanging a and I in our original relation and rewTiting similarly 
our triple summation, we obtain 

On ^nbtracting our original relation from the .sum of the two relations 
obtained in this way and denoting by Tp,;, the combination d. fmpr., 
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nitr p ~ wip, ,) of the first denvati\es of the elements of M, we obtain 

Tb# ,a“ ra*.»a’ I + mBja“ ,,o® , = 0 

(m denvmg .vhich we have used the fact that the two double sums 
niaflO® fO* „ and mefa” r»a* , are the same (why?) (Hint The matrix 
HI IS symmetric ) (Note that = Fp,,) This relation can be 
written in the equivalent form 

(Fb^ ,o" ,a® t + r»)o^ t ~ 0 

On multiplying this relation by z* p (the denvati\ e of x‘ with respect 
to o’*) and summing with respect to t, we obtain, since ^ p is zero 
unless j = p, m which event it is unity (whj ?), 

r«j pO,” , + mp^o* „ = 0 

Finally, on multiplying this relation by and summing with 
respect to p, where is the element in the jth row and pth column of 
A/-*, we obtain 

+ - 0 

We now differentiate the relation F.^ pO* , 0 * , + mpta" r* “ 0 with 
respect to interchange s and < m the result, and subtract the two 
equations obtained m this way Since o'’ ,t — o*’ u and » 
a* , 1 , we cfctam 

(Fb/j p) ^0“ ro' ,0'' t - (Fbs p) »a* ro" tOT , + F.jj^c" . 

— Fba pO" r,o^ I + trtptf ya“ rfO’’^ t fnps ri®"' t = 0 

We simplify this relation by first wntmg the triple summation 
(r,sp),a*‘^‘’,a’', 

m the equivalent form (r„, p) so* ,o' .o"^ j (why are these two tnple 
summations the same?) and then eliminating the second derivatives 
of the o’s with respect to the x's by means of the relation „ = 
— nv'''r „5 r®'* t On the appropnate rearrangement of the summa 
tion labels every term has a“ ,o* bo’ t as a factor, and then if we multi- 
ply through by a:’’ qX* ,x ‘ » and sum on r, «, and t there remain only the 
terms for which a = j, ^ j and y — k We obtain m this way the 
following relation 

(Xqj p) k ~ (Fgl: p) ^ ~ pTqk 0 "P *®**Fgj jF,* p 

— jfeFg; I + m*‘mp, }Vqk 0 — 0 

From the definition of the symbols Fp , , it readilj follows that mpg , = 
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r,, , -f r^r;p, on iiFing this to eliminate the first derivatives with 
respect to the a’s of the elements of the matrix 3/, we obtain 

Per./?- ~ (P«/:r)A “ (Pitip)./ + — m'T piicrgysf ~ 0- 

Tlw^e are the (exact) compatibility relations connecting the elements 
of M and their first and second derivatives with respect to the o’s. 
The terms involving the second derivatives of the elements of M are 
furnished by ~ 0\i-,p).) and arc ~ 

d- niqkjp). It follows readily that the cxtri'alurc tensor Rqp.jk is alter- 
nating in p and q tind that, it is symmetric in the pairs qp and jk. 
Tiiere are, then, tmly six distinct compatibility relations, namely, 

/?:n;23 = 0, R31;3J = 0, Ifl 2;12 = 0, 

/f23;31 ^ 0, /?31;12 ~ 0. /?12;23 ” 0. 

When the strain matrix is treated as an infinitesimal matrix, the 
symbols 1%.,., are infinitesimals of at least the first order and so the 
term'' involving the products of P’s in the compatibility relations are 
of at least the second order. Neglecting second and higher order 
infinitesimals tlio compatibility relations reduce to 

^^^pxc.tq T = 0 

or, equivalently, since M differs from 2t] In' a constant matrix, 

’tjc.'jt 'b ’tot.JF “ Vjt.pk ~b nph.iq- 

d'liere are three rclalion.s of the type 

’723.11 + ’711, 23 = ni2.31 + ’731,12 
and three of the type 

’722.33 + ItSS.CC — 2’723,23. 

In our original notation, wliere we used (a, b. c) instead of (o’^, a% n^), 
the six comp.'ifihility ndations (of the infinite^sirn.al theory of elasticity}^ 
are 
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EXERCISES 

1. Verify the (infinitesimal theory) compatibility relations by setting 
Vat *» Ua, Vbc “ + Vt), CtC 

2 Show that the (infinitesimal theory) compatibility relations are satisfied by 
every homogeneouv strain, i e , by every strain matrix whose elements are constants 
Is this statement valid for the exact compatibility relations? 

Ansiier Yes 

3 Show that the (infinitesimal theory) compatibility relations are satisfied by 
any strain matrix whose elements are linear functions of (a, 6, c) Is this state 
meat true for the exact compatibility relations? 

/Inswcr No 
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THE CONNECTION BETWEEN STRESS 

AND STRAIN 


L 'J'lic mntrlx 

Wo I'oii'^icicr ;iny tlncc-fiimonsional portion of our deformable 
medium (when in tlie final or deformed state) and we denote the sur- 
face of Vg by .S-. Since our deformable medium is supposed to be in 
equilibrium, when in the deformed state, Vi is in equilibrium under 
tfic action of various forces. These forces are of two tj'pes: 

(1) Mas.s or iiody forces (such as the weight of the material in 
Kr)- Ps being the density at a fyjiiea! point f*j.:(a:, y, z) of Vx, the 
innss of the element of volume dVx i^ px dVx and we denote the force 
|>er unit iims"; acting on this element of mass by F; the coordinates of F 
an* furnidicd by tlic elements of a 3 X 1 matrix 
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the elements of the 3X1 matnx 



and the force acting on Vx across the element of area dS® has its coor 
dmates furnished by the elements of the 3 X 1 matrix / dS, where 
dS* = {(dS*)*dS®)^* IS the scalar clement of area of Sx We term 
the vector whose coordinates are furnished by the elements of the 
3X1 matrix / the stress on the matrix element of area dS® Thus 
each coordinate of the stress \ector has the dimensions of a force 
divided by an area Stress is force per unit area 

We now make the basic assumption that these two systems of forces 
which maintain that portion of the deformable medium which occupies 
Vx in equilibrium would maintain it m equilibrium if it were rtgtdified 
This implies that the two systems of forces (mass forces and surface 
forces) which act on the deformable material occupying V* must 
satisfy, when taken together, the conditions imposed upon a system of 
forces which maintains in equilibrium a ngid body These conditions 
are most conveniently expressed as follows 

The virtual uork of all the forces of the system in any virtual ngtd dis 
placement is zero 

A virtual rigid displacement is a particular case of the more general 
concept of a virtual deformation In order to define the latter we must 
imagine that the coordinates (x, y, z) of P, are functions not only of the 
coordinates (o, b, c) of Pa but also of on accessory independent vanable, 
or parameter, which we shall denote bj $ The four variables a, fe, c 
and 6 are independent, and if / is any differentiable function of these 
four vanables w e wTite its differentia! m the form df + 5/ where 

df = fa da + fb dh + U dc, fi/ = /» dB 

Thus df IS calculated under the assumption that $ is held constant 
whereas 5/ is calculated under the assumption that o, b, and c are held 
constant A virtual deformation is defined when the 3X1 matnx 



IS given as a function of x, y, and «, i e , when is furnished at each 
point Px of the deformed position of the medium When the 3X1 
matnx fix is such that S(dsi)* — 0 we say that the virtual deformation 
IS a t irlual rigid displacement 
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Since the vnrinblca a, hi c and 0 arc independent the order in which 
fccond-fjrder derivatives with respect to tliesc variables are calculated 
in immaterial tihese second-order derivatives being assumed con- 
tinuous). 'i'liti.s, for example, .r (where x,as denotes tlie 

derivative with respect to 0 of x„, and T,ta denotes the derivative ndth 
respect to a of a-?) and, after wc muUipl}' through b}- dO, this relation 
appears in the form 5(x„) — (ot)„. If M is any matrix, of which the 
element in the ^tli row and 7)th column is wc denote by the 
matrix of which the element in the tfth row and pth column is 5m j,'’. If 
wc denote by (6x)r the 3 X 3 matrix 


/(5x), 

(ox)„ 

(ox).’ 

(5x), = ( (5)/). 

(5y)j, 


\(5r). 

(5r), 



(Sr, 5 1/, 62) 

(^, ?/, s) 


the relation S(x,,) = (ox)o and the similar relations obtained by replac- 
ing a by b or c and x by y or 2 yield 


5.7 = {5xhJ 


^ 1 * '71 

where ./ i.s the Jacobian matrix IProve this. Ilmt. (5x)„ = 

(a, b, c) 

(5i‘)/Xn -f (oxlyf/n + (5x),-„.] It is clear that 5/* = (5J)* (prove 
this), and St) 


gj* = J‘((5x),]*. 


riitieo Jf « J*J j( fullow.s- that 


5,1 f ^ + .7*5.7 

- + (5x),!J. 

Me denote hy D the symmetric matrix 


7) Jj|(5.r).}* -J- (gx)^.} 

/ (5x), i!(5y)..-f (5x),j 

I ^1(5;/).- -r (5x)„} ( 6 ,,)^ 

M -?• (5:)r i i ( ( 02 ). -f (5y), 1 

and then wc have the relation 


i{(5x)r + (02)., !\ 
^•!(52)y -f (5y),j j, 
(5r), / 


o.ir ^ 2J*DJ, 

or. eqtsiv.'slently, since .If Tij -b 2 tj 


J*DJ. 
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Since (dsx)^ — (da)*J\I da, we have [remember that c is independent 
of 0 so that 5 (da) =* 0] 

S(ds^y = (da)*MIda = 2{da)*J*DJ da 
= 2(dx)*D dr 

For the virtual deformation to be a virtual ngid displacement we must 
have 5(ds,)* = 0 for arbitrary choice of the 3X1 matrix di Hence 
ue have the following cntcnon for a virtual rigid displacement 
A virtual deformation Sx ts a virtual ngid displacement if, and only if, 
the symmetric matrix D = ^{[(ix),]* + (ix)*) is the zero matrix 
It follows that the system of forces, both mass and surface, acting on 
Vx must be such that their virtual work, m any virtual deformation for 
which D IS the zero matrix, is zero The simplest virtual rigid dis- 
placements are the virtual translations for which the elements of the 
3X1 matrix ix are constant functions of (x, y, z) Setting, m par- 
ticular Sx - 1, 5y »» 0, Sz =5 0, we find that the virtual nork of the 

mass forces is p^FdF* uhereas the virtual work of the surface 

Jy. 

forces J* / Hence 


Ly''-+L 


fdSx 


(which merelj expresses the fact that the resolved parts in the direction 
of the positive x axis of all the forces, both moss and surface, acting on 
Vx add up to zero) On letting the volume of F, — ► 0 we see that 

lim ( » 0) J* f dSx = 0, and on applying this result to a flat cylm 

dncal volume whose height tends to zero we see that / changes sign 
when the direction of dS* is reversed Appl>ing our result to the 
tetrahedron whose faces are the coordinate planes and a plane whose 
coefficient vector is v(u^ u^, u®), we find that 

where /* is the x component of the stress on an element of area whose 
normal has the direction of the positive x-axis, etc Wnting down the 
corresponding equations for g and h, we see that the 3X1 raatnx / 
whose elements furmsh the coordinates of the stress across the matnx 
element of areadS®whichhasthedirectionoftheumtvectora(«\M®,w®) 
is given by the formula 

/ « T« 
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the other hand, the stress matrix has nothing to do with the initial 
position of the medium, it is determined by the final position of the 
medium Thus it does not make sense to ask Is the stress matnx 
sensitive to a rotation of the initial rectangular Cartesian reference 
frame’ It is logical to inquire, however Is the stress matnx sensitive 
to a rotation of the final rectangular Cartesian reference frame? 


EXERCISES 

3 Shon that if the stress on aa element of area u always normal to the element 
the stress matnx is scalar le a multiple of E$ Note When T => —pEz 
p > 0 the stress on an element of area i9 always normal to the element and is a 
pressure rather than a tension Hus is the situation in hydrostatics and we term 
for this reason ascalar stress a ftydroslafic stress A hydrostatic stress is & pressure 
when r 13 a negative multiple of Ez and a tension when T is a positive multiple 
of Ez 

4 Show that m a hydrostatic stress the magnitude of the streaa on an element of 
area is independent of the direction of theelement of ares 

5 Show that the resolved part in the direction of the element of area dS* of 
the stress on dS* is u*Tu where u is the 3X1 matnx whose elements are the eoor 
dinates of the unit vector which furnishes the direction of dS* 


2 The conditions of equilibrium and the virtual work in any 

virtual deformation 

The virtual work of the mass forces acting on any portion V* of the 
deformable medium in any virtual deformation is furnished by the 

lolume integral J' pt(Sx)*F dVi and the virtual work of the surface 
forces acting on F* in this virtual deformation is furnished by the 
surface integral J* {Sx}*f dS^ Since f = Tu and since (d5i)ii = 
dS* the virtual work of the surface forces may be written in the form 
J (Si) *T dS' If, now, $ = (t, % f) is any 1X3 matnx the surface 

integral I fdS^ = J (Ids* + vdS" + fdS*) = (Idyds + 

r)dzdx + ^dxdy) may be written, by virtue of Greens theorem, as 
the volume integral I (£• + 1 , + f.) = Jj, Wiv, f *) dv^ where 

v,e understand bj the divergence (in the x-coordmate system) of any 

3 X m matnx the 1 X m matnx obtained by taking the divergence of 
each of its column vectors Thus 
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EXERCISE 

1. Fhow thnt when the ttre-"! Ik a hyrlrostatic pre«'iurc the equations of cqui- 
!i!)rium are -(Rrn'lr p)* -{- pif' = OamJ that the equations of motion arc -feradj 
/•)* ” tr'a ''here Rrad- p is the 1X3 matrix (})z, Vp, Pi)- 


On iisifip t!io fact tlmi 7'* = T, Ic., that tlic stress matrix is sjmi- 
mefrir. in ihv vsprtiji'-ion /or t))e virfuaJ work, wc obtain (show this) 
the fol!ov.'iiie futularacnfal result: 

'I'll' lirlual rrorl; of all (he Jorccfi acting on any portion U. of our 
ilffiirmnhh vudium in any virt^ial deformation is obtained by integrating 
the (rare of (he prodnrt TD over TQ: ^ 


( bw 

Virtual vork = f Tr{TD) dV^. 

.) r. 


r 

i. 


r>^ 


IVe have seen that og ^ .J*DJ and so TD = Since 

7‘r Alt « 7V It A it follows (show this) that 

'"Si 

7V(7T7) =. Tr{d~'^T{J*rHy). 


lienee oiir expression for the, virtual work of all the forces acting on any 
porti(»n I . of the deforniuhlc medium in any virtual deformation may 
he written in the ecpii valent form 
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/6‘\ 

and 6 = ( 6* I are two real 3X1 matrices, the 

w 

scalar product of the column vector of 6 by the column vector of a is 
the 1X1 matrix a*h This scalar product may be written as Tr a*b 
or, equivalently, as Tr b*o If A and B are two real 3X3 matrices, 
either of the two equal numbers Tr A*B, TtB*A may be termed, by 
analogy, the scalar product of the two 3 X 3 matrices A and B Thus 
the virtual work, m any virtual deformation, of all the forces acting on 
any portion F* of our deformable medium is found by integrating the 
scalar product of and Sjj over F, It is clear that this 

scalar product is (as the name implies) msensitive to a rotation of the 
final rectangular Cartesian reference frame In fact, under the rotation 
S’-* aj' = R*s, J' « R*J, T-* T’ = R*TR and so 

J-^RR*TRR*{J*)-^ = thus is insensi 

tivQ to a rotation of the final rectangular Cartesian reference frame and 
we know that ij, and hence iij, is insensitive to any such rotation 
In the infinitesimal theory of elasticity the (approximate) result con 
cernmg virtual work may be phrased as follows The virtual work of 
all the forces acting on any portion F* of the deformable medium m 
any virtual displacement is found by integrating the scalar product of 
T and Jtj over F* Since T is sensitive to a rotation of the final rec 
tangular Cartesian reference frame, whereas h is not, it is clear that 
the scalar product Tr(T5j?) of T and $ij depends on this final reference 
frame, this makes it evident that the statement of the infinitesimal 
theory of elasticity concerning virtual work cannot be an exact state 
ment, for the virtual work cannot depend on the accidental choice of 
the final rectangular Cartesian reference frame 

3. The virtual work when the stress is hydrostatic 
The expression for the virtual work of all the forces acting on any 
portion F* of the deformable medium, m any virtual deformation, 
takes a particularly simple form when the stress is hydrostatic, i e , 
when T = — pRj In this case = —pJ~\J*)~^ = 

— p(il/)”^ and so the virtual work is found by integrating the scalar 
product of — p(il/)"^ and St} over F, Since M = Es 2i? 51/ = 
25 t 7 and so 

Virtual work = — J J p Tr{M~^SM) dV^ 

If, now, A 13 any 3X3 matrix, 5(det A) la the sum of three determi- 
nants (each obtained by differentiating one of the column vectors of A) 
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502* 

03^1 



003^ 

ifioi* 


fl3*i + 1 


oaz" 

fl 3 “j * 1 * pi* 

Oo* 

503- 

Isoj® 

02^ 

03 '! 1 

Oi^ 

ba-A 

03^! 


0:^ 

503® 


/prove this), and this result mn}' be written in the equivalent form 
5(det A) = Trico AYhA. 

If A is non-singular we obtain, on dividing through b}' det A, 
7V(/1-*5.4) = (del A)"^5(det A) 

and HO the virtual work maj-, when the stress is hydrostatic, be written 
in the form 

Virtual work = --I- f pfdet Af)“^fi(dct M) dV^. 

~ J V. 

Since dot M = (det ./)" wo have 

\'^irtual work = — j* jjfdct ./)~'5(dcl ./) dV^, 
dV * 1 

and since dot J = 7 —' fio that 5(dct^) = o{dVx), wc obtain, 

(i\ (I ^ » a 


(itially, 


Virtual work 


■X, 


7/6(dFx). 


If J) is a constant function of (r, ij, z), i.c., if the stress is a constant 
hydrostatic. pres.-:ure (or tension), this equation simplifies further. On 

writing V'- -- f d]\ ^ f dFa, we obtain oF^ = I 

Jy. Jy,dV„ Jy. dVa 

j oidVf). (Note that, in the integration over V„, the limits of 
J 1% 

integration arc independent (jf the accessory variable, or parameter, 
(h where.n.s in the integration over I', the limits of integration vary', in 
general, with 0.) Hence, when the sires.s is a constant hydrostatic 
pres.-im', or tension, we have 

Virtual work = — p j 5(rfFr) = — poF-. 

J y. 


4. The clnMic energy and the relation lJet^s•ecn .«;trcps and strain 

We CA ume that the work (hme by all the force.s acting on any portion 
Ti of our deformable medium in any deformation i.*; .stored up in as 
r.vuftV er.i'rfry or er.enjif of deformation (none of it being dissipated in 
hwu'; euf'.'-K)’ o! deformation may well depend not only on the 
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actual state of strain of the medium but also on the previous historj 
of the medium, m other words, the energj of deformation may be 
distnbuted throughout Vg with a mass density if' 

Energy of deformation = J.. Px'f' dVx 

where is a function not onlj of ij but also of the previous values of 
Tj (from 0 to ij) Such hysteresis effects are matters of common expcn 
ence but we shall ignore them here We make the fundamentiil 
assumption that the mass density ^ of the energy of deformation (i e , 
the energ} of deformation per unit mass) is a function of the strain 
matrix y so that if/ is unambiguously determinate whenij is given 

if/ ~ if/(ri) 

Remark 1 It is more conicnient to deal with the energy-per-unil 
mass rather than the energj -per unit-volume since the mass of nns 
clement of volume dVt of the medium remains constant m anj defor- 
mation whereas dVg changes, in general 

Remark 8 Since >7 is a symmetric matrix, ^( 17 ) is unambiguously 
determinate when the six following elements of i? are given — t 7 a«, 
»766» i7e«> viei Ve»t ’Jaft — ^^^d WO Can Write 

if/ “ |/'(l7aa> l74kj Veti *7f«> ’laj) 

However, the gam obtained in reducing the number of varmblea of 
which ^ is a function from nine to six is overbalanced by the loss of 
symmetry involved in this reduction We shall agree to write any 
fiinition of the symmetric matrix v sytnrmtncalhj wlicrc we under- 
stand that /(» 7 ) IS written symmetrically if /(A) = /(A*) where A is 
any 3X3 matnx If /(j 7 ) is any function of 17 it may be w ntten sy m- 

mctncally by setting /( 17 ) = / ^ pC’j)* obvious that 

g(A) = { 7 (A*) where A is any 3X3 matrix (prove this) The mam 
advantage of wnting a function /( 17 ) of ij symmetrically is that we do 
not lo‘-e symmetry when vve take the gradient of / with respect to 17 
Just as the gradient of any function of (z, y, s) is the I X 3 matnx 
obtained bv taking the three partial dcnvatives of the function with 
respect to x, y and z, so the gradient of any function /(A) of a 3 X 3 
matnx 1 isthe3 X 3 matrix obtained by taking the partial denvativcs 
of /(A) with rc«pect to the vanous elements of A Thus the clement 
m the ptii row and 9 th column of the gradient of /(A) with rc‘*pcct to A 
13 the partial dcmalivc of /(A) with respect to a,® In calculating 
these partial denv ativ cs we pay no attention to the particular structure 
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of A ; v.e regard the nine dements of A as independent whether they 

actually are independent or not. Thus, in particular, if .4* = A wc 

regfird'a^'" as constant wlicn wo arc differentiating with respect to 

H « even though we know that a/ = a/. Wc denote the gradient of 

5f . d/ /dfy , . 

f(A) by the .symbol ”• It is clear that — = { — j and so we have 

the following important re.sult: 

///(A) u! n /unciion, tmUen symmclricaUy , of the symmetric vmlrix 
Of 

/I then —7 is, I the .4, a symmetric matrix. 
dA 

Evainplc 1. The first invariant /j = i)ar. + yih + yec Is a function, 

dl 

written .cvrnmolricallv, of the symmetric matrix y. ~ is the sj'-m- 
' ‘ ay 

metric matrix Ej: 


Oy 


Es. 


Esaniple 2, The second invariant /j = iyi,iycc ~ "niryd) + (yeeVaa 
— ’trrt’tnr) + {vn»yf.b “ yoiybn) is a function, written .symmcl rically, 

<7/" 

of tlic symmetric matrix y, is the .symmetric matrix IiEs — yt 


din 

dy 


/lE, 


Exnmjde 3. The third invariant h ~ dot j? is a function, (.sup- 

Of 

po-^odly written symmetrically), of the symmetric matrix y. — - i.s the 

dy 

;.yjnmolric matrix co y: 

dis 
dy 


CO y. 


The prmeij)!(> of cnn.-crvnlion of energy may be cxprc.<scd as follow,s: 
the virtu'il work of all the forco.‘« acting on any volume 1'^ of 


our 
Pz'P tlVx. 


deformable juodiurn, in any virtual deformation. = 5 j 

. r ' , . r 

t)n writing J in the equivalent form j p,t'pdV„, we hav 

^ ^ fr |j. PjOvt'dFs., and so tlie virtual work 

li.<% the integral over V, of "r^eyf) i. the same- as the 
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integral over T* of Pxi'p Since this must hold for an arbitrary volume 
Vx of our deformable medium, ue have (why?) 


Since f IS (by hypothesis) a function (written symmetncally) of ij, we 


have (by the very definition of a differential) = Tr 
this), and so 




(bhov 


Try-'ryj-'j,) - p.rr 

Since this relation must hold for an arbitrary (symmetric) matrix 
5i}, we have (why’) 


or, equivalently, 






On setting <t> - p«^, so that ^ is the energy of deformation per unit 
tmltal volume (why?), we obtain, finally, 



This is the fundamental relation of elasticity theory connecting stress 
and strain If tf> is known as a function of 17 and */ J is given, we can 
calculate the stress matrix T from this relation (ij being determined by 
means of the formula 17 = — Ez) and the compression ratio 

being determined by means of the formula — = ;) In the 

pa uet J 

infinitesimal theory of elasticity the difference between J and Ez is 

neglected, (so that the difference between the compression ratio — and 

Pa 

1 IS neglected) and so T = — In words, 
drf 

To the degree of approximation afforded by the infinitesimal theory 
of elasticity, stress is the gradient of energy density mth respect to 
strain, that is, the stress matrix is the gradient, with respect to the 
strain matnx, of the energy of deformation per umt initial volume 
Be sure to understand that this is not an exact statement, it is only 
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an itpproximnf ion to Uie exact result, %vJiicliis furnished by the formula 

\pa/ (’V 

Since ./*/ ~ d/ we have J - wiiere R is a rotation matrix; 
hence the fonntila furnishing T may be written as follows: 


7’ 


- RM-^ — AP^R*. 

Pn HV 


On comparing this with the relation 

JJ* = RA/R*, 

we obtain tlic following result: 

fn iJial particular final rectangular Cartesian reference frame in which 
the cxiordinalcs of .hi* are furnished by (he elements of Af the coordinates 


of T arc furnished by the elements of 


M 


dr] 


If, in particular, the initial rectangular Cartesian reference frame is 
so cho'-’cn that A! is in diagonal form, we see that: 

In the final rectangular Cartesian reference frame that is furnished by 
(he principal axes of .JJ* the coordinates of T are furnished by the ele- 
ments of 




Hit, 
nil r— 


(mjTth>)'- 

. d<t) 
drjuh 

{niimf}'- 


(t(t> 

dtp,., 



’ drp,], 


(m'lnisf^ 

(mjw,)’’* 

d<tt 

(m3?;!;)'' 

t d<i> 

d6 

?;i3 


djjrn 

dr]fh 

di]cc 
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of unit magnitude and directed a^vay from the medium) Since 
dSxU = dS" we have 

TdS* « dSxf 

Since T — J~J* and 

dS- = (det J) (J*) ■ rfS* = (J-)-' ds- 

(See Chapter 1 Section 8) this relation may be ^\ntten in the form 
dA 

J = dSxf 
on 

e denote by Ta the 3X3 matrix J — and obser\ e that Ta dS“ — 
dn 

T dS*, m \\ords, the result of operating on the Jinal matrix element of 
area dS‘ uitk the stress matrix T ts the same as the result of operating on 
the initial matrix element of area rfS® tnlh the matrix Ta Since, m a 
deformation problem, it is the initial matrix element of area dS* which 
is given (the determination of the final matrix element of area dS* 
being part of the problem) it is mth the matrix Ta and not with the 
stress matrix T, that we must operate when we Mash to satisfj the 
boundary conditions These conditions are 

Ta dS" = dSJ Ta^J — 

dn 

If, for example, part of the surface of our deformable medium is free 
from applied force f is the zero matnx over this part of S* and so 

Ta d)S“ = 0 Since J is non singular it follows that — diS“ == 0 and 

dn 

dtft 

so dS“ is a charactenstic vector of — corresponding to the character 
dij 

istic number zero Note In some problems the stress matnx T 
(rather than the stress/) is furnished over S* Thus the medium may 
be subjected to a constant hydrostatic pressure Then the boundary 
conditions simply state that T must comcide, over Sx, with the given 
stress matnx For example, an important problem that we shall later 
examme thoroughly is the followmg Determine the deformation of a 
circular cylindncal tube under the action of internal and external 
pressures p, and Pe Assuming the displacement of each particle of 
the cylinder to be radial and a function of the distance of the particle 
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from the- of tho cylinder, (he deformed tube is again a circular 
rylindric.'d tube. The vector TadS° must have (he direction of the 
radius and, if r. and r, arc, respectively, the iiniial internal and external 
radii of (he tube, (he magnitude of Ta dS'^ must be pidSx when r = 

n and P( dSg when r « r,. From (he relation dS‘ = dS°, 


we obtain 


dSx 
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ISOTROPIC ELASTIC MEDIA 


1 . The energy of deformation of on isotropic elastic medium 

We have seen that the strain matrix 17 is sensitive to a rotation of the 
initial rectangular Cartesian reference frame Under the rotation 
a — * a' =* R*a, n — » = K*ijR This change of the strain matnx does 
not mean that the rotation of the axes affects the strain or deformation , 
the strain is the same as before but its coordinates (1 e , the elements of 
the strain matnx) with respect to the new axes are, in general, different 
from Its coordinates with respect to the old axes The ene^ of 
deformation, per unit initial volume, must be insensitive to the rota- 
tion of the initial rectangular Cartesian reference 

frame This fact implies, since ij is sensitive, in general, to such a 
rotation, that the form of the function which furnishes the eneigy 
per umt initial volume, must be, m general, sensitive to a rotation of 
the initial rectangular Cartesian reference frame In the new refer- 
ence frame the energy of deformation, per unit initial volume, will be 
furnished by a neio function of the new strain matnx and the 
fact that the energy of deformation, per unit initial volume, is a scalar 
(1 e , that it IS insensitive to any rotation a a' = R*a of the imtial 
rectangular Cartesian reference frame) assures us that 

<#’(’!) = = 4>'(R*vR) 

If it should happen that the form of the function 0 is insensitive to a 
given rotation a~* a' — R*a of the initial rectangular Cartesian refer- 
ence frame, we say that the medium is ehashcally insensitive to the 
rotation R Thus the medium as elastically insensitive to a given 
rotation R if, and only if, <i>'{R*tiR) = il>(R*riR) for an arbitrary sym- 
metne matnx ij In view of the relation ^(»j) = tf>'(R*T]R) we can 
rephrase this definition of elashc symmetry (with respect to a given 
rotation of the initial rectangular Cartesian reference frame) as follows 
The deformable medium is elastically tnsensitive to the rotation a 
60 
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a' = ira of Ihc initinl rectangular Cartesian reference frame if, and only 

if, ' 

thin relation being an identity in thr (symmetric) matrix g. 

Iiiifrad of rcf^arding flie cleinonLs of g' ~ R*gR as furnishing the 
roordinntes of tlic strain with respect to tlie new axes we may consider 
t)ic neic strain who«e coordinates with respect to the original axes arc 
the same as the coordinates of the old strain with respect to the new 
O-xes (so tiiat tiic coordinates of tljc new strain with respect to (he 
original axes are furnished hy the elements of g' — R*gR). Then the 
U'lation <t>ig) = <}>iR*‘gR) may Ik* interpreted as indicating that the 
< nergies of deformation, per unit volume, of these livo (in general, 
<lifTeren() strains are the same. The new strain is that which would 
he produced if the original strain were preceded by the rigid rotation 
a ■— a' " R*a of the medium. Note that (his is an actual rotation of 
the mcfiiiim and not mcrclj* a rotation of the initial rectangular 
Carte,^ian reference frame. 

A defonnahlo medium is said to be elastically isotropic (or simply 
isotropic) if it is elastically insen.sitivc to every rotation of the initial 
reeuuigolnr Cartesian reference frame. Thus, 

A deformable medium is isotropic if, and only if, 
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30 30 

isotropic medium “ “ 77- -Ej + ~ (/1E3 — »?) + ~ co i;, hence 

01J oil oij 01 3 

when the initial rectangular Cartesian reference frame is so chosen that 

n (and hence V and 1 /*^) are diagonal so also is — It follows that 

3 jj 

— commutes with so that 1/^*— If** = If — Hence, 

3>j 3ij 3ij 

For (in isotropic medium the coordinates of the stress tensor T arc 
furnished in the final rectangular Cartesian reference frame in uhich tie 
coordiniJ/es of JJ* are fumtsked by the elements of M = J*J, by the 

/pj\ 50 

elements of the matrix I ~ T" 


In the final rectangular Cartesian reference frame that is fur 
nished hy the principal a\es of JJ* T is diagonal (it being understood 
that the axes of the initial rectangular Cartesnn reference frame are 
the principal axes of M = J*J) 



30 30 . 30 30 I 1 

where = “ 4 - (vi + ns) ry + vns-ry etc Thus the pnncipal 
3 »ji 3 fi dl» 3/3 

axes of T are the principal axes of JJ* In an arbitrary final rectangu 
!ar Cartesian reference frame ue have 



where J “ RU** Note The pnncipal axes of ij arc the principal 
axes of J*J whereas the pnncipal axes of T are for an isotropic 
medium the principal axes of JJ* Thus the principal axes of T arc 
not even for an isotropic medium in general the same ns the pnncipal 
ixcs of ij If howerer \\c follow the deformation bj the ngid rota 
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^ "*” shall not carry the development of 

<f>, as a power senes m Ji, 7j, I3, farther (so that we agree to neglect, 
in the development of mfinitesimala of higher order than the third 
in the elements of the strain matrix 17) We have, then, 

^ = ^0 + ^>1 + + 4 > 2 t 


<j>i = all, 


d<t>i 

drf 


= Bi, 


02 = ' "t ^1* “ = yJiEa + 2fiyj\ 

o df] 

= (f/x* — 2 mJ^E 3 + 2Tn/ii7 + n co n, 


“ = aEz + (X/x^s + 2>(17) + (Z/i* — 2m/2)£3 + 2nilii} + n co 17 
aij 

Hence T - y(aEj + (XI, B, + 2m) + (lli‘ - 2mI,)E, + 

2«7itj + n CO ij}y* When the strain is zero, 17 is the zero matri\ 
and J IS a rotation matrix On denoting by To the inihol stress 

matrix, 1 e , the stress matrix when 17 = 0, it follows (since — « 1 

Pa 

when »7 = 0 [why?]) that 

To ^ aE} 


In other words, the initial stress ts either a hydrostatic pressure or a 
hydrostatic tension The only assumption about the deformable 
medium which we have made is that the medium is isotropic and this 
has been enough to force the conclusion that the initial stress must be a 
scalar stress, 1 e , a hj drostatic pressure or tension Hence, 

No elastic medium that is tm/tally tn a stale of stress u-hich ts not 
scalar, te,of the nature of a hydrostatic pressure or tension, can be 
isotropic 

Assuming that the medium is isotropic, we denote the initial hjdro- 
static pressure by po, then a = —po and the stress matrix T, when the 
strain matrix is 17, is furnished by the formula 

r = 7| -poBj + (XI,E, + 2m) + (!7i= - 2mh)E, 

+ 2m7i»7 + n CO 1717* 
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formula iiivnlvp.s the nix consiniits po, A'l ^ ^ (tlic significance 

of the l''nn eonsfant being that tlicy are independent of the strain 
matrix v)- Of tiiese six constants the first is tiic initial pressure and 
so its dependence on the initial pressure is obvious. It is natural to 
expect, then, that the remaining five, namely, X, g, I, vi, and n, depend 
also on po. These five constants are known as the clastic constants of 
ihc isotropic medium. In the infinitesimal theorj', which neglects, in 
the formula for the .stress matrix, inf)nitc.sim.al.s of the second order in 
the ele.ments of the three cla.stic constant.s (/, in, n) do not appear, 
and T i.s furnished by the formula 

7- = h) J\-poE, + (X/i^.i + 2pv)\J* 

Xpo/ 


wherf; the compre.«sion ratio — = I ~ I The two elastic constants 

Pa 

X and fi which ni)pe;ir in thi.s formula are knovoi as the clastic constants 
of Lami> (after Gabriel Lamti [1795-1870], a French mathematician). 

Warninrj. Do not be mi.slcd In' the tci'm constant into overlooking 
the fact that X and p depend on the initial hydrostatic pressure. 

Writing J ^ we have, to the degree of ni)proximation con- 

lemiilatcd by the infinitesimal theory, 

irrn - (i - /,)(F?, + 2>,)i -po/^s + (x/,/^.-, + 2^*?)! 

-- —ptiEs + (X -i- po)J lEs + 2{p — po)i?. 
Wlien the initial .strc.«s' is zero, Ihi.s- e<iuation reduces to 
U*TR ~ \qIiEz + 2pov 

where (Xo, po) an* the values of (X, p) that correspond to po = 0, i.e., to 
?eni initi.'il j)rc.<sure. On comparing this with the previous relation, 
whieh may be written in the form 

d* PdEs)R = (X -r pojliEs d* 2{p — po)ij, 

it i;. clear that the numbers Xn and pc, that appear in the formula Rir- 
nidsing the .stress, when the medium is free from .stre.ss in the initial 
state front whieh the .strain is measured, are the values of X -f- Po 
.and p ~ p,-,. n\ss>octivcly, when po is zero. If. for example, we neglect 
th.e tiependence of X and p on po the effect of a non-zero initial hydro- 
••tatie pressure is, in so far as the increment T d- poEz in the strc.s.s 
m.atrix !- coneomed, to inerr fisr ,\ by po and to decrease p by po. 

in the Mxond-order a]»proximation. in which we neglect, in the 
tormuht for ?, tlurd-ortier (but keep .second-order) infinitesimals in 
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the elements of the strain matrix ij, the formula lorR^TR is 


R*TR - (Ej + 2 i 7){ -poEj + (X/iEj + 2,i^) 

Pa 

+ (III* 2t7iri)Cs + 2ra/\»j + n co ij| 

nlicre now ^ - (1 + 27. + 4/,)-“ = 1 - /i + (| /,* - 21 ^ On 

collecting terms, we obtain R*TR = — poEj + (X + Po)IiEi + 
2 (ii - po)’J + 2(X + po — M + »t)/i + I(I — X — ^o)h^ — 2 (ni — 
Po)/2lC3 + 4/n?* + n CO 1? When the medium is initially free from 
stress this reduces to R*TR « (X0/1E3 + 2/ioij) + 2(Xo — ao + tno)Ii 
+ }(lo ^0)^1 ~ 2moIi]E3 + 4/101}* + no co ij "here {lo mo no) are 
the values of {I, m, n) that correspond to po = 0 

In comparisons of theory inth experiment tho compression ratio 

— 13 often furnished by experiment and it is unnecessary to use the 

Pa 

approximation (1 — 7i in the linear theory and 1 — /i + (f/i* — 
2 / 2 ) m the second order approximation) If "c do not use these 
approximations v\c obtain the folloiviog results 

lATtear Theory 

R’TR = (B, + 2,)l-p,E, + (X7,Bj + 2m)! 

" (f) i 


When Po = 0 this reduces to 


R*TR 


(XofiEj + 2/101?) 


Second order Approxtmalion 
R’TR = (£. + 2,)| -p, 775 + (UiB, + 2m) 

+ (I/i^ “ 2 mIt)E 3 + 2mlir] + n co ijJ 
“ t ~Po2^3 + hliEj + 2{n — po)i? 

+ - 2 m/j)E 3 4- 2(m + \)Iiv + n co ij + 

When po = 0 this reduces to 

R*TR =5 1^0^ 1^3 2/4(1^ + (Io7i* — 2171(3/2)^3 

+ 2(mo + Xo)/ii7 + no co IJ "H ^/loij*) 



SCALAR STRAIN OF AN ISOTOPIC Ml'.DIUM 


07 


3. The ?lrc«s in nn inotropic medinm wlien in a stale of scalar 
strain 

We shall write our scalar strain in the form 


7j = — cF?. 


(the reason for the negative sign being that most of the experiments 
with wiiieh we wish to cheek the tlieorx* deal with comprci^sion rather 
tlian dihitnlion. If we wrote i; = cEz the coefTicient c would be 
()()ritivc in dilatation and negative in compression; when we write 
t) ^ ~rKz the cochieient r is positive in compression and negative in 

dilatation'). Since the medium is i.solropic, — is a scalar matrix (i.c., 

dt] 

a multiple of Fj) when n is a scalar matrix (prove this) and it follows 
(why?) that R*TR is a .‘•calar matrix. lienee (why?) T is a scalar 
matrix (being the .same .as R*TR), and on .setting 


T 


■pEz 


(so that p is the pressure in the medium) we obt.ain the following rela- 
tions between c and p in the linear theory and the second-order approxi- 
mation, respectively; 


J.iiunr Theory. 
written ns follow.s 


Since (1 — 2c) 



the formula for p may be 


p ^ 



IPif "b (3^ + 2p}cj. 


^\hcn we u'-e the approximation 1 + 3c to ~ this equation becomes 

Pa 
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m other words, the compressibility is the derivative with respect to p 
of the loganthm of ~ Terming the reciprocal of the compressibility 


the tncompresstbthty, r\e have the following result 

To the degree of approximation fumuhed by the infinitesimal (or 
linear) theory of elasticity the (local) incompressibility is given by the 
formula 

Local tncompressiMlily = X + -f + -^o 


When the initial pressure is zero the local incompressibility is Xq + ^fto 
Note carefully that the (local) incompressibility is sensitive to the 
value of the initial pressure If we neglect the dependence of \ and /i 
on po (so that rve replace X and p by Xo and no, respectively), the (local) 
incompressibility is the linear function Xo + %no + Po If we 

make the less drastic assumption that X and p are linear functions of po, 
the (local) incompressibility is a linear function Xo + -fpa + ^Po of 
Po (the coefficient k of po being a constant that must be determined by 
experiment) Assuming that the strain and po are constant over the 
homogeneous medium, ve have, then, 

(dp r)p-p, Xo + fpo + kpo 

Since Po IS arbitrary, this equation is equivalent to the relation 

A| 1 « I 

dp V Xo + fpo + kp 

On integrating this relation, we obtain the following equation of state 
for the isotropic deformable medium 

I ^ _ (ZsX 

\F/ 

or, equivalently, 


where Vo is the value of V which corresponds to p = 0 The value 
of k which corresponds to the (drastic) assumption that X and p are 
independent of po is •^, and the correspondmg equation of state is 

p - (3X + 2rt |(^)“ - ij 
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The Sccond-ordrr Approrimatton. Here 

p = Ip, (3X -p 2n)c - m -f 7;)c=}, 

and (his reduces, wiien po = 0, to 

p ~ {(3Xo + 2/io)c — (9fo + 7io)e'j* 

Sinee = (1 - 2c)~^‘ = 1 -f c -p ^ c* -p • • • , these results 

\Pa/ 2 

ma}’ be written in the equivalent forms 

p = Pa + (3X -p 2;i -p 7 Jo)c -p (3X 2/i — 01 — n -{- f-po)c‘, 

p « (3Xo d* 2/io)c -1- (3Xo + 2fio — 01 o ~ 7io)c". 

Limitinn ourselves to the case in which po and c are constant over the 

medium, we have ~ ~ = 1 d- — and so (1 — 2cf' — 1 -p — 

Pz In y 0 Fo 

2 


ir 1 . 

Hence c ~ ~ — h 


•11 / . J.2AE 1/An2 

, OT>,v„lo„,l,-, , - 2r _ (l + --) - If + 

— 2 

— ) -p • • • and so the formulas 
0/ 

ows: 

Ap tr- _ d- ^ p -p “ 7)0^ 


for p may ho written as follows: 

I \ ^' 
Fn 


+ 5(x + |--2(-?»+^p.)(|L^ 


P = - (x. + I ^ + ^X. + ? - 2 (. - ^ ((-!■)■ 

Fiop Y^i*- denoted, for the moment, by it follows that 


dp 

d 


P A , 2 1 \ d-p 2 ,2 4 

s “'P+i;''+F) rf?= + + 


On dinerentiatinc the first of the,'-e relation*- with resjTect to $ and sub- 
'■tstutint: the rv>ult in the second, we obtain 
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and if this value is di\ ided bj the expression for ^ it follo^vs that 




+ “Im) ~ 21 — ip 
^ + f + ip 


For values of p -which are negligible m comparison wth Xo + iuo a 
good approximation to the fraction on the right is 

2( 9fo+nol 


3Xo *h 2/xo) 


Thus the A, that appears in the equation of state 
(Xo + Imp) .) 


'm 


of the linear theory is (approximately) I — 


2(Ofo + np) 


The fact 


3(3Xo + 2po) 

that Xo + Iwo maj be determined by experiment (as the value of 
^^^hen p “ 0) furnislies a method of determining by experiment 

the combination 9lo + no ^ 

Remark 1 We have termed the equation 

(Xo + W f/FoV ,1 




the equation of state of the linear theory Actually, the equation of 
state of the linear theory is the equation 

p = {3Xo + 2/(o) e, 

which we obtain by setting po = 0 m the equation 


P = 


Ipo + (3X + 2n)e\ 


If we use the approximation — = 1 + 3e this linear approximation 

Pa 

appears in the form 

p = (3X|) + 2tto)e = ^Xo + - ^ ~ 
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cquniion is in tli'^agrecmcni v.'ith experimental results unless ~ 

j ‘3 infinitc'imal. and wc shall not stop to discuss it. The equation that 
we prefer to call the equation of state of the linear thcon' is rcallt’ the 
equation of state of .an iutrgratcd linear theory. We use the linear 

AE 

fheoiy to cxqtre.ss A?) as a multiple of — not merely at p = 0 but also 

at an arbitrary value of p; the multiple X d- la + an vnhwwn 
function of p, and we approximate it by a linear function of p (the 
eo'-fTieient of p in this linear function being determined by experiment). 
The final relation between p sind V’ is obtained by integrating the rela- 
tion between dp and dV obtained in this way. 

Itrmnrk 2. The formula p - \po -f (3X + 2a^)c — (91 + 

n)c*} of the .second-order .approximation may be checked again.st the 

relation 7 ^ — rr^ wlicrc'h = Fn'^ is the energy of defonnat ion of the 

entire medium. In fact, on selling /j = ~3c, In = 3c‘, h = — c® in 
the formula 

X -4- 9a 1 -i- ‘>m 

« ~pc/i d 1\- ~ 2filn 4. /j3 _ 2tnIiIo d- Tj/a, 


we obtain 
and .so 


<t> ~ 3pne d* |f3X 2/i)e‘ — (9/ d- «)c" 
«!• - V„ \Zpor d- ?(3X -b 2a)r- - (91 -b 


/FA" 


I 1 wo have e 

J 

2 1 V''./ 1 

iW> 


Hence - — . - 

jr 

(A iP. + ( 


“ S \v.) ir- - .n-, ^ (i^) If- + 0^ + 2.)« - 

!W + n)f-l (f ) ip« + (•'» + ‘-'(.V - (iM + Ti)f=l. 

1. Cotnpnrivon of theory with experiment in the en«e of extreme 
IijKlroHfatic pres‘.urr 

The cxjH'riment.s of Bridgman on the eomproKsibilities of variovm media 
up to the extreme prts-.ure of 10^ utmo-pheres ( ^ 10^ grams per S(|urire 
cuUimetcr) are available to le-t the various relations connecting 
pO'..!ire v.ith voUmu’ whidi an* funii-hed by the various theories we 
have convIderiHl (tie- .simple linear themy. the integrated linear 
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theory and the second order approximation) Bridgman* gives the 
V 

relative volume — at pressures up to 10® atmospheres (1 atmosphere 


being I kilogram per square centimeter) The values for the metal 
sodium are as follo'ws (the unit of pressure being taken as 10^ atmos 
pberes = 10^ grams per square centimeter) 


V = 


25 

3 

4 

5 

6 

7 

8 

9 

10 

Z.l 

Vo ' 

^ Vo) 

- 0 789 

0770 

0 737 

0 70S 

0 683 

0 661 

0 641 

0 623 

0 606 

“ 7o 

V 

" ^ Vo 

= 0211 

0230 

0263 

0 292 

0 317 

0 339 

0 359 

0 377 

0 394 


It la immediately clear that the formula furnished by the (simple) 
linear theory is in disagreement ivith the experimental results This 
formula is 


P 




and if we use the experimentally determined value of ^ at p « 2 5 

Vo 

to determine the value of the (single) constant Xo + § 1*0 appeanng in 
the formula we find that 


Xo + ho “ 11 86 

The follmving table gives the comparison between the experimentally 

V V 

determined values of and the values of — calculated from the 
formula 


p •= 


JL 

Vo 


1 +^ = 1 , 

Po Xo + Imo 


25 3 4 5 6 


- 0 0844p 
7 8 9 


10 


Proc Am Acad Aria Sn Vbl 76 pp 55-70 (1948) 
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\Mien p ^ 4 


the value of 


_F 

Fo 


given by the formula is too small by 


about 10 per cenf.. The experimental values of — have an accuracy 

f 0 

of 2 or 3 per cent, and so the simple linear thcon.' is not valid up to a 
pressure of *1.10’ atmospheres (les.s than twice the pressure used in 
determining the constant, ^tl + of the formula). Y e shall not con- 
sider further, therefore, the formula furnished by the simple linear 
theory- and shall pa-ss to a study of the formula provided by the inte- 
grated linear theory. 

There arc two con.stants, namely, k and Xo -f |go in the formula 


_ (^0 + t Mo) 1^— y — ]| 


of the integrated linear theory. If pi and pi arc two pressures to wliich 
corrcsipond the volumes Fi and 1 ^ 2 , respectively, we have 



and so pz ~ ~ 0>2 - P\) - 0. The value of k 

may lie determined from this equation (by the method of trial and 

Y 

error) by using the experimental values of ~ and — : k having been 

1' 1 V 2 

determined, the value of Xo + g/‘o may be obtained from the equation 


, , 2 kpz 

Ao + r ao - TTTTT" 


In orrler to .see how the integratwl linear theory works both backwards 
.and fonvjmls, i.e., when we extrapolate to lower a.s well .ns to higher 
values of p, we select the values pi = 5, p. ^ 0 in order to determine 
the two con*-tants, h and Xf, -r gpe. of the formula. Our calculation 
will s'Tve to dcl'cmiine the value of the integrated linear iheorv a.s a 
mean'*- of pri-'.iidion; ilms we m.ny imagine that Bridgman had suc- 

Y 

c*>s!fxi in me.nsuring ~r only up to p G (i.e., up to a pre,=.-ure of 
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CO, 000 atmospheres) and we wish to predict the values of — up to 

Vo 

p = 10 (i c , up to a pressure of 10® atmospheres) If our predicted 
values check well with the experimentally determined values we may 
use these \nluos to predict the (as yet experimentally undetermined) 
V 

valucsof — uptop = 20 (le, up to a pressure of 200, 000 atmospheres) 
The equatioh that serves to determine K is 



Y y 

where — = 0 708 — = 0 683 On denoting the left-hand side of 
1 0 Vo 

this equation by f(k) we find that (check these calculations) /(3) = 
0 214, /(4) - -0 09G5 /(3 8) = -OOOC, /(3 70) « 0 001 In view 
of the errors of experiment the value 3 70 for k is close enough Using 
this value of k we, have 



22 74 ^ 
3 242 ' 


7 014 


The formula (of the integrated linear theory) that furnishes — is, 

r 0 

accordingly, 




+ 1 851) 


The following table gives the companson between tho values of 
computed from this formula and the experimentally measured value 



p = 26 3 4 

6 6 

7 8 9 10 

V 1 

— (measured) = 0 780 0 770 0 737 
Va \ 

0 703 0 683; 

0 661 0 641 0 623 0 606 

V 

— (calculated) =. 0 798 0 776 0 73S 

0 703 0 883 

0 662 0 643 0 637 0 613 


The greatest divergence between the calculated and the measured 
values IS at p = 2 5, where the calculated value is in excess of the 
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vnluc by lots (Ivan 1.5 per cent. Since the experimental 
valtiLv do not claim an accuracy of more than 2 or 3 per cent -vve may 
fay that the values furnished by the two-constant formula of the 
integrated linear theory check the e.xperimcntally measured values 
over the entire ranfvjo of pre.^sure-s from p = 25000 atmospheres to 
fi ■■ 100,000 atmo.spheres. 

This check between theorj* and experiment gives us confidence to 
V 

predict the values of “ that will be measured when the experimental 
Fo 

icclmiciuc is sunicietifly developed that mea.surcment.s of ~ up to 

p 20 (i.e., up to a pre.^sure of 200,000 a(mosphcrc.s) can be made. 

We ti'-e the experiment alh’ moa.sured values of at p = 0 and at 

p r- 10 (o rodctcrminc (he values of I: and ,\o + |pd. The equation 
thill serves to dclonnine h is now 


10 



« 0 


Pi P" 

where ~ -023, — - .GOG. Wc find that (perform the calculation 
' 0 » 0 

or verify the result) 1: - 2.80, Xq -f = S.9GS. Note. The differ- 
ence lietween the.u’ values for /: and Xo + ajuo and those obtained 
previously (when pi wa-s 5 and wa.s G) merely indicates the differ- 
ence between the linear .approximation to X -j/li in the interval 
0 ^ p *7 10 and the linear approximation to X -f in the interval 

^ . P 

5 t; p ^(3, The formula determining-— is, now, 

Po 

- (0.310)(p -1- 3.135) 



and the comparison between the experimentally determined value.? and 
the eomiuited value.s over the range p ~ 2.5, 3, • • • , 10, together 
nith the jiredietcd values over the range 7> — 11, • • • ,20, i.sfurni.shed 
by the table on jv. 70. 


Thu- til*' theory predicts (hat a pressure of 200,000 atmosplieres will 
(■omptC"' sodium to 0.-197 of its original volume. At a pressure of 

P 

SOO.WK) atmospheres the formula yield.s- == 0.100, and if wc have 

* 

roafidoace rnuuch in the formula to use it at p = 100 xve find that 

P 

the c<wn-sponiling value of is ().29.->. In other word.-:, a prc-ssurc 

• f 
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of 10* atmospheres will compress sodium to less than three^tenths its 


original volume 


The agreement between the values of computed 
Vo 


from the formula and the expenmentally measured values, is withm 3 
per cent over the entire range p = 2 5, 3, , 10 of the experiments 

The greatest discrepancy is at the pressure 2 5, where the calculated 


value of — IS 0 815 and the measured value 0 789 This measured 
'0 

value IS somewhat uncertain as it was obtained by an adjustment 
between two senes of experiments with different apparatus Over 
the range p = 4 to p — 10, for which direct readings from one appara- 
tus were obtained, the agreement between the values given by the 
formulas and the experimentally measured values is \vithin 2 per cent 
(the discrepancies for all the values save the first being less than 1 per 
cent) 


P - 35 


4 

5 

6 

7 

8 

9 

10 

y 









— (measured) 0 789 
Vo 

0 770 

0 737 

0 708 

0 688 

OCCI 

0641 

0623 

0606 

■— (calculated) 0 815 

0 782 

0 760 

0 715 

0688 

0 6C4 

0 642 

0 623 

0 606 









p- 11 

12 

13 

M 

16 16 

IT 

18 

to 

20 


y 

— (calculated) 0 591 0 577 0 564 0 552 0 541 0 531 0 522 0 513 0 605 0 497 


Remark A noteworthy feature of the formula connectmg p and 
V 

— which IS furnished by the integrated linear theory is the unsym- 
Vo 

metneal way m which it treats pressure (p > 0) and hydrostatic ten- 
sion (p < 0) Thus it IS clear from the formula 

that, as P — » 0, p — ♦ w and, as P — » , p — » — — ° 

other words, the medium cannot, according to the integrated linear 
theory, support, without rupture, a hydrostatic tension of amount 

» although it can support an arbitrarily large hydrostatic 
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|)re‘i-itirc. From the formula in which tljc constants (for sodium) were 

determined by using the experimentally measured values of — at 

p, 5 and p: ^ 0 it would appear that the medium cannot support a 
hydrostatic tension of 1.851 X 10^ atmospheres, and when the con- 
stants are calculated from the experimentally measured values of 
V 

“ 7 - at Pi ~ fl and pz ~ 10 thi.s mplurivg hydrostatic tension is raised 
I'o 

to 3.135 X 10‘ atmospheres. 

We now turn to the fonnula furnished by the second-order approxi- 
mivlion : 

j) == !(3Xo + 2ao)e — (OIo + uolc'} 

/ y \ 

~ ^ ~ “(OIo ri* ^*o). 

/ F\’‘ 

Here 1 — 2 c = ( ’p' j • sodium may be presented 

n.s follow.H: 

j. ~ 2.5 3 -1 5 C 7 8 9 10 


£ 

i”p 


0.780 0.770 0.737 0.70S 0.08.3 O.CGl O.Gtl 0.023 O.GOG 


0,0731 0.0700 0.0020 0.1028 0.1122 0.120G 0.1283 0.1353 0.H20 


GO'*-- 


310 2.750 3.013 4.-1.5G 5.2,81 G.09S 0.898 7.G87 8.4G2 


re have 


Denoting, for a moment, by g the product of p by ( 'iT'j • ''c 

the following two equation.^ to determine the two con.stants « and 0 of 
our forsmila; 

ocj ri- ^er qi, 
nrz -}- ~ 7;, 

t’j and e* being the values of <• which correspond to any two values pi 
and ;h, r<‘>]K‘otively, of p. The fommla-s that .‘■erve to determine a and 
|:1 are. acconlinglv, 


- I'i) 


^ - 


7;fl - 7,0; ^ 
ryzici - fi) 
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isothopic elastic media 


In order to test the \ alue of the second-order approximation for the 
purposes of prediction we a^am suppose that we are m possession of 
experimental measurements up to p = 6 only and ne use the formula to 
V 

predict the values of — corresponding to p = 7, 8, 9, and 10 Since a 
yo 

and /3 are very sensitive to small vanations in ci and ez if ez — Ci is 
small, we use the values pi = 25 and pj = 6 m determining a and 
/3 We find (check this calculation) that 

a = 2 619, /3 = 390 3 


The results of the comparison between theory and expenment is 
furnished bj the follow mg table 


p(obscrvc<l) = 

25 

3 4 5 

1 0 

7 8 9 10 

V « p I ~ ) (observed) » 

2 310 

2 750 3 C13 4 450 

5 284 

1 6 098 6 89S 7 687 8 462 

^(calculated) 

2 300 

2 730 3 605 4 457 

5 282 

6 080 6 859 7 608 8 362 

p(calculatcd) ~ 

2 409 

208130805001 

5 DOS 

6 980 7 056 8 608 0 882 


No(e For reasons of convenience we have calculated the values of 
p that would correspond (when we use the formula of the second 


order approximation) to the cxpenmentally measured values of tt 

KO 

rather than, as with the formula of the integrated linear theory, the 
V 

values of — that w ould correspond to the experimentally controlled 
'0 

values of p The greatest divergence between the calculated value of 
p and the experimentally controlled value of p occurs when the 
experimentally controlled value is 10 and the calculated value is less 
than this by less than 1 2 per cent Since the experiments do not claim 
an accuracy of more tiian 2 per cent, we may say that the formula of 
the second order approximation yields correctly the connection 
V 

between p and — up to p = 10, i e , up to a pressure almost twice 
r 0 

as great as the larger of the two pressures that were used in determining 
the constants of the formula 

If we determine the constants a and 0 by means of the expen- 
mental measurements at pi = 2 5, ps = 10, instead of by means of the 
experimental measurements at pi = 2 5, pz = 6, we obtain the 
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follov.'ing values; 

a - 1.902, P = 400.2. 

'Plio conspariFon heUveen the c.vpcrimcntully measured values and the 
values obtained from the formul.a 

p^-Vy - 1.902c + 40G.2€= 


is furnished by the following table: 


jiOitnf rvrd) 

*2..') 

3 4 .5 0 7 8 9 

10 

{vY 

ij ^ f’ 1 T" ) (obsprcod) 

2.310 

2.750 3.013 -1.4.50 5.284 0.098 0.S9S 7.087 

8.402 

q (('alciilalial) 

2.310 

2.740 3.015 4.480 5.327 0.13S 0.931 7.094 

n 

j) (raleiil.nted) 

2..'’)00 

2.990 4.002 5.037 0.019 7.010 S.03S 9.008 



'riius the e.ah'ulated valuc.s of p differ from the obsen'od values of p 
by le.'^s than 1 per cent (i.c., by lc.«.s than the c.vpcrimcntal margin of 
error) over the entire range from p = I io p — 10. 

To find the value of p which would correspond to a given value, saj', 

r , /V\'^ 

0.197 of ~r we proceed as follows: From ( — ) “ 0.G275 we obtain 

' 0 ^ \I'o/ 

r - 0.b‘>ri2 and so p - 14.44, p ~ 18.2.3 (as compared with the 

value ji 20 furnished by the extended linear theory). 

Xof/'. Tiie valuc.s 2.019 and 390.3 for « and /5. respectively, deter- 
mined from the experimental ob.*-ervation.s at pi ~ 2..o and pj = 6 
differ from the values 1.902 and 400.2, re.spcc lively, determined from 
the mca.surements at- pi - 2.5 and p; = 10. This merely indicates 
that the assumption that f/ is a quadratic function of r over the range 
0 ^ p 5; 10 is only an approximation to the truth; the pnrahnlir. 
cpprorwiaiioii over the range 2.5 ^ p 1$ G diffei-s slightly from the 
p-im/folic npproxinutffnn over the range 2.5 ^ p ^ 10. The fact that 
*3 i** '■n large compared with o indicates the irnitortancc of considering 
the third-order ifastie eon.sl.ant'.-. Whereas Xo -I- Ipo is around 0.7, 9/ 
t: is negative and around —400. 
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In order to test the \ alue of the second order approximation for the 
purposes of prediction we again suppose that ^\e are m possession of 
experimental measurements up top = 6 only and n e use the formula to 
V 

predict the values of — corresponding to p =* 7, 8, 9, and 10 Since a 
To 

and /3 are verj sensitive to small variations m ei and c? d c? — ei is 
small, we use the values pi = 2 5 and p 2 = G m determining a and 
i8 We find (check this calculation) that 


a = 2 619, ^ = 396 3 


The results of the companson between theory and expenment is 
furnished by the following table 


p(observed) •» 

25 

3 4 5 

e 

7 8 9 10 

fvY 

9 - P 1 ^ 1 (observed) «• 

2310 

2 750 3 613 4 450 

5 284 

6 098 6 898 7 687 8 462 

^(calculnted) > 

2 309 

2 739 3 595 4 457 

5 282 

6 080 C 859 7 608 8 862 

^(calculated} •• 

2 499 

298139805 001 

5 908 

6 080 7 056 8 DOS 9 882 


Note For reasons of convenience uc have calculated the values of 
p that would correspond (when we use the formula of the second 

V 

order approximation) to the experimentally measured values of 
rather than, as with the formula of the integrated linear theory, the 

V 

values of — that would correspond to the expenmentally controlled 
►0 

values of p The greatest divergence between the calculated value of 
p and the expenmentally controlled value of p occurs when the 
expenmentally controlled value is 10 and the calculated value is less 
than this by less than 1 2 per cent Since the experiments do not claim 
an accuracy of more than 2 per cent we may say that the formula of 
the second order approximation 3ueld3 correctly the connection 
V 

between p and — up to p = 10, i e , up to a pressure almost twnce 
'0 

as great as the larger of the two pressures that were used m determimng 
the constants of the formula 

If we determine the constants a and by means of the expen 
mental measurements at pi = 2 S pj = 10, instead of by means of the 
experimental measurements at pi = 2 5, p 2 = 6, we obtain the 
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following values: 

tt = 1.902, /3 = 40G.2. 

'Ilic comp.'iiri'-on between the experimentally measured values and the 
vulims obtained from the foitnula 


- L()02c + 406.2e- 
r fnrni-'hed by the following table; 


gfi.li'enetl) 

~ 2 .r> 

3 4 5 0 7 8 9 

10 

/cV^ 




,j„j, (oh-ervodi 

2.:uo 

2.7.50 3.013 4.150 5.2Si 0.098 0.898 7.087 

8.462 

'1 ff'ile'ilnUai) 

2.310 

2.746 3.015 4.4S9 5.327 0.138 0.931 7.691 

S.4G2 

;» fr'iK'ciatoi) 

2.:)00 

2.00r» 4.002 5 037 6.019 7.0 16 8.038 9.008 

10.00 


Thu« the ealoulated valuer? of ?> differ from the observed values of p 
i>y le.'''; than 1 per vent (i.c., by less than the experimental margin of 
error) over the entire range from 7) = 1 to p = 10. 

O'o faui the value of p which would rorrc-spond to a given value, say 

r /r\'^ 

0.497 of ™ we proceed as follows: From f — j = 0.G275 we obtain 
e • . (1,}SG2 :u)(l so p p ^ 18.23 (ns compared with the 


vuhie ;i - 20 furnished by the extended linear theory). 

XoU . The v.'duc-s 2.1519 and 39(5.3 for « and j5, respcctivelv, deter- 
mined from the exjK'ritnental observations at pj = 2 5 and'p, = 6 
differ from the values 1.902 and 400.2, respectively, determined from 
1 he measuremmUs at p, .= 2.5 and p, == 10. This merely indicates 
(hat the a-muption that q is a quadratic function of c over the ranee 
(I A p 10 is only an ajiproximation to the truth; the parabolic 
cppri.M.,n oyer the range 2.5 ^ p ^ G differs slightly from the 
appuKurmlmi over the range 2.5 .< p ^ 10. The fact th.t 
an. M. large compareil with a indicates the importance of considering 
tn- Uisnl-ordcr ohastii- eon.stanls Whereas X J- " : siaering 

, • , ojns. *' ucreas Ao -f- -j/io (s around 0 7 97 

— f. {• negative imd around ~.]00 ^ 
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1 . Definition of a non-isotropic elastic medium 
Certain materials (for example, wood) that are important m con 
struction are not isotropic i e , they are sensitive to at least some rota 
tions of the medium (before application of the strain) We have seen 
also that even if a medium is isotropic v. hen it is initially free from 
stress or when the initial stress is a hydrostatic pressure or tension, it 
cannot be isotropic when the initial stress is any stress other than such 
a hydrostatic pressure or tension It is, therefore, highly important, 
from the scientific point of view, to obtain the relation connecting 
stress and strain for non-xsolroptc media Unfortunately, houever, 
this relation is much more complicated than for isotropic media, even 
m the mfimtesimal theory many more elastic constants than the two, 
\ and a that suffice for isotropic media are needed In the second 
order approximation the number of elastic constants becomes, usually, 
so large that the theory is too complicated for practical applications, 
but it IS possible to construct an tnlegrated linear theory that is useful m 
practice Many deformable media have a crystalline structure so that, 
although they are not isotropic, i e , el^tically insensitive to all pre- 
lotations of the medium thej are elastically insensitive to certain such 
prerotations Thus the relation 

ip(R* 7 jR) = 0 ( 9 ) 

Mill be nou valid for certain gi\cn rotation matrices R but not for 
all rotation matrices ^Vhen 0 is analyzed into a sum 

0 “ 00 + 01 + 02 + 0s + 

of terms of different degrees in the elements of »j (0j being a homo- 
geneous function of degree 3 in the elements of j = 0, 1, 2, ), 

the relation 0(fl*J7H) = 0(17) implies the senes of relations 

0 j(R*i 7 R) = 0,(17), 5 = 0 , I, 2 

(Pro\e this Hint The relation 4>{R*iiR) = 0(17) is an identity in 
the elements of It remains true, then, if 17 is replaced by A.7 where I 
80 
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i? an arliitmrv' i)arametcr; when this substitution is made <}>j(ri) is 
replaced by The relation >i>j(R*vR) = imposes certain 

conditions on the coofiicients of the homogeneous function, of 

flegrcej, of the elements of rj (save in the ease j = 0 in wliich we are not 
interested (why?)) and we proceed to find the nature of these condi- 
tions. for various given rotation matrices R, Avhen j = 1 or j = 2 or 
^ '™ 3. 

EXERCISES 


1. .Siiow tlin? if n (.Irformnble incflitim is ctaslicnlly insensitive to n piven rotation 
R it is clantirally insensitive to every integral power, positive or negative, of R. 

iteplare 17 liy A’ and by R^R* in the relation 4>{R*rjR) >= 

2. Kliow fljnt if the nu'diiiin is elaslicnlly inscnsitK'c to a given rotation R it is 
rla"(ieally ins^n'-ilive. to the reflexion —R. Xole. A reflexion is an orthogonal 
Iraftsforination of rcetangnlnr Cartesian coordinates the determinant of whose 
matrix is —1 (instead of +1 ns for a rotation). 

3. .‘^how that if R is a rotation through t around the a-nxis, — /? is a reflexion in 
the fi-platie, 

•i. Sliow that if the medium is elastically insensitive to two given rotations 
Aland Ajit is ehsstieally insensitive to their product AiAj. 

.*>. Flmw that if the medium is elnsfienlly insensitive to a rotation through x 
around the n-axis and to a rotation through *• around the h-axis it i.s elastically 
iUN-tiKilive to a rotation through x around the r-axis. Hint. .Any onn of these 
rol.ation' is the protlucl of the other two. 


(>, .‘'how that if the medium is elnslieally insensitive to a rotation through 


2 


around the o-iixi'- and to a rotation through ^ .around the h-axis it is elastically 


insrnsitive ti, a rotation through ^ around the r-.a.\i.s. //in/. 


Denoting these 


rot.’stiiui*. Ity /I’l, /,*;, nnd Ri, respertively, R3 «= /?iA;Ai*. .Vofe. Certain eiibie 
rrysiflls j>ov>.<--s the elastic symmetry of this exercise. 

4. that if a medium po<se.sses the eln.stic symmetry’ described in Kxerci“e 


^ . .... 2r 

0 It 1" el.’Wti'-.ally insensitive to a rotation through ” around each of the four 
iliaroniib of the eiilie whose vertiees ore the point.s (il, il, il). Hint, The 


produet of a rotation through - around the fi-axis by a rotation through ^ around 

*>_ 

tlie r-asi- !■. a rotation through " around one of the diagonals of the cube in 


2 . lilt* form of ii>i(r}) for variou«i rotations R 

n iinc-ar jttnetton of tlic dements of r,. AVlicn written sjtw- 
mttricnlly (what does thi« mean?) it appears in the form 


- rric.v) 
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%\here C\ is a symmetric 3X3 matnx The elements of C\ are the 
values, vhen ij = 0, of R*TR for 

R*TR = ~ 

and 'when ij = 0, — = 1, M = Ez, and-;^ =» C\ Thus, 

Pa 5*7 

The coefictents of the linear function of the elements of ij, when 
this linear function ta mitten symmetrically, are the coordtnaJes of flie 
initial stress in the reference frame tn which JJ* — J*J On denoting 
the initial stress matrix by To ue hate 

R*ToR = Cl, R = 

We now denote the variables Vbii Vee, Vie, Vea, Vab) by (in, vs 
V3> Vi> Vi, Vi) and write in the form 

^l(’?) ” Cim + Czm + CSVZ + C4?74 + C5176 + C6I76 

Thus the symmetric 3X3 matnx Ci is 

( Cl ici 
ici Cl ici I 
ics ^4 cz/ 

Under the rotation a-^a' = R*a, R*riR, and since oa*-* 

a'a'* = R*aa*R it follows that the elements of ij transform, under the 
rotation a— » a' = R*a, like the elements of oo* This fact iaiph<s 
(show this) that (»;i, ijj, 175, 174 175, 17,) transform like (0*, h^, c^, Iw, co, ct) 
Let us consider the case in which R is a rotation through an angle 8 
around the a-axis, since the a axis is arbitrary, this takes care of every 
rotation There is no lack of generality in taking 0 < ® ^ since R 
may be replaced by J2* and we shall do so We have 

/I 0 0 \ 

R = I 0 cos e — sm d 1 
\o sm $ cos 9/ 

so that a' = a, b' = (cos (»)6 + (sm fl)c, c' = -(sm e)b + (cos 8)c 
Hence 17— ♦ 17' where 

iji' - i?i, vi' = (cos* 0)172 + 2(co3 9 sm 0)174 + (em* 6)vt, 

173' = (sm* 0)172 - 2(cos 0 sm 0)174 + (cos* 0)173, 

174' =* -(cos 0 sm 0)172 + (cos* 0 — sm* 0)174 + (cos 0 sm 0)173, 

175' = (cos 0)175 — (sm 0)i7e, ly,' a (gm e)rji + (cos 0)i76 
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Sinro t;;' -f ra - vz + vz it is couvenicut to introduce, as replacing 
tlx: variables (iji, • * - , jjo), the variables ($i, $4, fs, fc) "’here 

ii “ VI, = i'ivz d- m), fa = Hnz ~ vz), 

f4 = Vi, f.'> ~ Va, fc = ’76- 

Then Tji " it, vz =•' fa + h, Vz - sZ - fa. Vi = f4. Vi, ~ fs, Vo ~ fr. 

and, under tin; rotation « o' = 7?*o, ^ where 

fi' s!f f;^ ~ fsi fa^ “ 20)^3 + (sin 20)^.i, 

U' - “(sin 20) h d- (cos 20) U h' = (cos 0)^5 - (sin 0)U, 

tp' = (sin t?)f5 d- (cos 0){c- 

'I'lieso relations suggest the introduction of the (complex) variables 
(Cl. f:. Ca, i' 4 , Cs. io) defined by the cciuations 

Cl *-"■ fit h ~ fs, fa — ddfs — i^i), 

ii ~ i-(f?. d- i^i), fr. '= ddfr. — ffr,). fc = -iCfr, d- (fr.)- 


Of (h(‘ six variables (fi, • • ' , fe). fi nod fn arc real, and fj and f 4 arc 
conjugate eomjdcx numbers ns arc also fp and Cc. It is clear that 
fi "■ fi, fi " f:, fa = fa d- f4t f4 = f(fa — ft), fs = fs d* fc, fo = 
j{fi - i'() and it follows (show this) that, under the rotation a —* a' = 
!i*n, f "> C' where 

fi' "= fi, f:' = f:, fa' = c=‘'f 3 . f4' = c-=*'f4, 

fr,' = c-’Ys. fr,' = c"ff.. 

Tlie relations connecting the original (real) variables (iji ■ • • 77c) with 
iIh- new (complex) variables (fi • • • j'c) are 

fi ft; s'! = fi - Jji: 

'' h d- fa f; d- fa d- fat fi “ fa ~ dd’J? d- Vz)', 


iz - f j - f; - fs - fa == 2 (fa - if4) ^ -^ivz- Vz) - ^ >74; 


I 


tt ' 


Ufa - f4); 


f 4 - o 'i' 's *) = “ (»?•• — 573) d- - ’74: 


9. a'*- • 2 

fs '= i>(h — ifij) - i-(vo — I'vo)] 

4 *'• vS *= ((is “ je)t fe = ■^(fi d* Jfc) -- d(77:. -b ivc)- 

y»e nenote tlio coefneients of when this expression is written as a 
tnndion of the variables (fj. ■ ■ • , j-^). hy (,/j, . . . , , 1 ^) 


tr fj = fs -f fp; 
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= Citji -1- C 2’72 + carj3 Citii -f cgijs + c$T]6 — dif 1 + dafs + djfs + 
^4^4 + dsfs + defe di and d2 are real, and ds and ^4 are conjugate 
complex numbers as are also d$ and dg The necessary and sufficient 
condition that the medium should be elastically insensitive to the rota- 
tion R, as far as is concerned, is that the relation 

diti 4 - + defs' = difi + + dfife 

should be an identity in the six vanables (fi, , fg) Upon sub- 
stituting for (fi'j , fsO their expressions in terms of (fi, , 
fe)) obtain the following relations involving (di, , dg) and 
the angle 6 of the rotation R 

e=‘*d3 = dj, c-'Vs = dg 

Note In writing these equations we have put down only one of each 
pair of conjugate complex equations, thus the equation e^'^dj = dg 
implies the equation c~***d« *= d4, and the equation c“‘*dg » dg implies 
the equation e'*dt “ dg 

If fl is not the quotient of 2 t by 2 , i e , if P ir, neither of the num- 
bers e®’*, «”** IS unity and we have d$ — 0 dg = 0 The medium is 
then elastically insensitive, as far as is concerned, to a rotation 
through any angle around the o axis (why^) and is a linear 

combination of fi = iji and « -^^(112 -I- ijj) We have, then, the 
following result 

If a deformable medium t$ elasttcalhj insermiiie to o rotation through 
an angle 6 9^ v around the a'Cxts it t$, as far as is concerned, elosh- 
cally insensitive to a rotation through any angle around the a-axis and 
^1(57) 18 a linear combination of iji and 1J2 -f- ns 

0 i(»?) = Pi’ll + C2(n2 + ns) 

Thus the medium cannot be elastically insensitive to such a rotation 
unless the initial stress matrix is, in the reference frame in which d/* = 
J*J, diagonal 

EXERCISES 

1 What la the necessary and sufliaent condition that the mednim be elastically 
insensitive as far as ^1 is concerned, to a rotation through any angle around the 
b-axia? around the c-axis? Hint When the (n b c)-exes are subjected to a 
cycbc permutation the two sets (123) and (456) are subjected to the same cyclic 
permutation 

2 Show that if the medium IS elastically insensitive as far as ^1 13 concerned to 

a rotation through any angle around the a axis and to a rotation through any 
angle around the b axis it is as far as ^1 is concerned, isotropic Hint ^i(s) is of 
the form c(t|i -j- ns + ’Jj) “ fVofe Tlie phrase ' as far as ^1 is concerned’ 

may be omitted since every rotation may be factored into the product of three 
rotations, two ofwhichare about thea axis the third being about the b axis 
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H 6 ~ r, c''* “ 1 and da is urbitran* (in'^tcad of beinp; zero). ^1(77) 
i', !KiVi’ H linear combination of fi, f;, and or, equivalently, of 
»j|. V:, Vh Vi- 

^liv) ~ CiVi + CiVi + ^3773 + 0^77.1. 

The initial atrc.«s' inatri.v (in the reference frame in which JJ* = J*.T) 
ift of the form 

/ci 0 0 \ 

I 0 Ci h^t ]• 

\0 ic} r.3 / 

The nccca«ary and snfiieient condition tluit the medium should be 
cla^ticjilly insensitive, as far as <l)i i.s concerned, to a rotation through - 
firnnnd (he o-axis or, etiuivalcntly (why?), to a reflexion in the a-plane, 
is that <^m( 77) “ <^iVi d* + 03773 + c^tji. 

EXICUCISES 

3. WJitu is till' M(‘rf't,>3ry njnl HafTicienl condition tlmt tl>c nicdiuni l)e eliislicnlly 
ifiscn- iti\ (•,})•. for ns iscnnri'rnofl, ton reflexion in tlie />-plnne? in the e-i)!nne7 

4. Sho'v ilfat if the iiunliuiii is el.astieiiHy insensitive, as far ns is concerned, 
to ft rotntioa ihroiiKli an jinRle (7 / r around tlie «-n\i.s it is elnstic.ally insensitive, 
ft»fnrfls v‘'i i>" coneerntnl, to .<i reflovion in each of tlie coordinate planc-s. 

3. The form of <^i:(7;) for various rotations R 

2'f’j(’)) F a homogcneou.s quadratic function of the elements of 77, and 
so it may be written in (he form 


(cjj> 7 i* d' ’ ‘ ‘ d' Cf.f, 77 c‘) -f- 2(012771170 cr.oVsVa)- 

Tlie mefiium is ela.stic.ally insensitive, ns far as <^>o is concerned, to a 
given rotjilion R if, and only if, the relation 
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independent, m this e\ent ne say that the medium possesses fewer 
than twenty-MDne second order clastic constants An isotropic medium 
possesses onlj two second-order elastic constants, the degree of com- 
plexitj introduced into the theorj of elasticity by lack of isotropj maj 
be grasped when we compare this with the twent>-one second-order 
elastic constants of a medium that is completely non isotropic (as 
far as ^2 Js concerned) 

We now consider the case m which R is a rotation through an angle 
0 around the a-a\ts, and ive introduce the complex variables (f 1 , , 

fs) of the preceding section We denote the coefficients of 2<f>2iri), 
when written as a function of (fi, , f$), by dpq 

2 ^ 2 ( 11 ) = diifi® + 4- dcafa* + 2(di2fif2 + dsersfe) 

Here du and djs are real, and dsz and <£44 are conjugate complex num- 
bcFb as arc also da and dee, di 2 , d 34 , and djg arc real, du and du arc 
conjugate complex numbers, etc The necessary and sufficient condi- 
tion that the medium be elastically insensitive, as far as is con- 
cerned, to the rotation R is that the relation 

+ 4- + 2(di2fi'i’2' + +ds8f'6'i'80 

= 4* 4* deef«* 4- 2(dj2fif2 -f +dB8f8{’e) 

should be an identity id the sit (complex) lanabJos (fi, , t«) 
This furnishes the following relations mvoUing the coeHicients dp, and 
the angle 9 of the rotation R 

e*'*dii ■“ dz3t e*‘*dii = du, c“’*d2s == d26, 

c~**'d8s = dss, = du, c'^das = djs, 

e*'*d23 = d2j, e*‘*d36 = dj# 

^ote In writing these equations we haxc put down only one of each 
pair of conjugate complex equations, thus the equation c^'^djj = djj 
implies the equation e“'‘'*d 44 «= d 44 etc 

If 0 is not the quotient of 2» b> 2, 3, or 4 none of the multipliers 
€*'*, <4*“', c**^ e®*' is unit> and so all the (complex) numbers djs, dss, 
dit, dis, d^s, d 2 s, djsr djs are tero In this event the medium is elastic- 
allj inscnsiti\e to aery rotation around the a axis and it possesses at 
most 21 — 2(8) = 5 second-order clastic constants, ^ 2 (>j) being a 
linear combination of fj*, fif*, fjfs, and fsfe or, eqiu\alently, of 
’ll*! (’J 2 4- > 33 )*, ’Ji(’7! + vi), im — irj)* 4 - and j/s* 4* Since 
»7i(»7J 4- ’js) * /i»ji - vi*t hvi may be substituted for ruirjt -{• Vi) 
this linear combination, and since ijs 4- »jj = /i — nu /i’ may be sub- 
stituted for (tj 2 4- 7 j)* Since ( 1 J 2 — nj)* + ® (173 4- »jj)* “ 
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~ »!<*). niTfS ~ nV tmy be Exibstitutcd for (r,; — tjs)* + 
iiiui, firifdiy, since /; ~ (r;:?72 — Vi“) + ^i('?2 + >/s) ~ ('!s~ 4 * ^c‘)i h 
may bf? suiwtituted for + ve"- We have, then, the following result : 
// a deformable medium is elaHicalhj insamtirc {as far as 4>2 w coa- 

2r* s- 

errnrd) to a rolntwn through an angle other than tt, -— •> or - around Urn 

O u 

ti-arh, it ir. rla.sticalhj inrenshuT {as far as <i>2 is concerned) to every 
rotation around the a-aj’is. Such a medium possesses not more than 
fie (inrtead of hrenty-onr) second-order clastic constants, and <f>2{i]) is a 
lir,rnr comlaiuition of the foUotriny five functions: 

Ir- fivi, vins — V4‘. 
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*?i> V 2 , VZ} >14 (involving ten coefHcients) and a linear combination of ij 5 *, ij*®, and 

v&vt 

7. What IS the form of <Js(i}) if the medium is elastically insensitive to a rotation 
through T around the 6 aias (or, equiv^ently, to a reflexion m the 6 plane)? to a 
rotation through n around the c axis (or, equivalent!} , to a reflexion in the c plane) 7 
8 Show that if the medium is elasticallj insensitive to a rotation through r 
around the a axis and to a rotation throi^h *- around the E> axis it possesses not 
more than nine second-order elastic constants, ^t(n) being the sum of a general 
polynomial of the second degree in the three variables tji, ijj, jjj (involving six 
coefficients) and a linear combination of ij4*, ijs* and ijj® 

When ® tte multiplier e^** is unity and so das is arbitrary, 

instead of being zero The medium possesses (at most) 5 + 2 = 7 
second order elastic constants, and ^2(1?) *s a linear combination of the 
five functions /i*, Iz, iji®, Ivtu ViVi ~ nnd, m addition, the two 
functions (ijj — rts)vi + ivi ~~ Vi)v 6 ■* (which are, respec 

tively, the real and imaginary parts of Sfjf g) 

EXERCISE 

9. Show that if a medium which possesses the elastic symmetry described m 
Exercise 8 is elastically insensitive to a rotation through ~ around the line jommg 

tho points (0, 0, 0) and (1. 1, 1) it possesses only three second-order elastic constants 
and that for such a medium ^{(i)) is a linear combination of 7t and (iniiz + 
+ ViVa) iftnt 4i2(jj) is invariant under the transformation ij-* ij' “ 

where 



thus ^ 2 ( 11 ) is invariant under the transformation »jj' = >12 172 ' = »ls, ’is^ *“ flli 
Vi' = Vi, Vi “ vt Vi = Vi Note Certam cubic crystals possess the elastic 
symmetry of this exercise 

When ® ^ unity and so dsj la arbitrary, instead of being 

zero The medium possesses (at most) 5 + 2 = 7 second order elastic 
constants, and ^2(7) is a linear combination of the five functions /i*, 
Ji, ’ll*! Ii’li, ’J2»?s ~ 174* and, m addition, the two functions (172 “ ’73)* ~ 
4 i 74 ®, (»72 — ’73)74 (which are, respectively, convenient multiples of the 
real and imaginary parts of fj*) It is easy to see that .,2'73 anfi 74* are 



they may be substituted for >72173 — 74* and (72 — 73)* — 474* Thus 
we have the follow mg result 



rnv: voibi of mv fof vahious hotatioxs n so 

The mcdhm is dmUcalhj imcnsttivc (as Jar as <^>2 is concerned) to a 
rotation throiiffh ~ arounel the a-axis if, and onhj if, 4i:in) « linear corn- 
hinaiion of the following seven functions: 

h', 1 2 , iji’f fi^u ViOSi Vi', (v: ~ V3)vi- 
EXEUCISKS 
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<^3 = iUcui’71® + ) + 3(C1121J1*92 + ) 

+ 6(ci23’7l’j2r?3 + )1 

There are fifty-six coeflacients Cp,,, p ^ q ^ r, m all six of the type 
dll, thirty of the type dis, twenty of the type dsa Thus if is 
•sensitive to every three-dimensional rotation R the medium possesses! 
(as far as <^>3 is concerned) fifty-six elastic constants, ve term thc'ie 
third-order elastic constants A medium that is completely non iso- 
tropic as far as ^2 and ^3 are concerned possesses, then, 21 -f 56 = 
77 elastic constants of which twenty-one are second-order constants 
(1 e , coefficients of ^ 2 ) and fifty-six are third-order constants (1 e , 
coefficients of ^ 3 ) The six coefficients of (or first-order elastic 
constants) are the coordinates (m the reference frame m which JJ* - 
J*J) 0 / the initial stress We proceed to investigate what reduction is 
caused in the number of the third-order elastic constants vhen the 
medium is elastically insensitive to one or more given rotations 
Introducing the complex variables (fi, , fe) of the preceding 
sections, Me denote the coefficients of 3 ^ 3 ( 1 ?), vhen this is Mnttcn as a 
function of the vanables (fi, , fe)» t>y dp,,, p ^ ^ = 1, 

. 6 

303(»?) = duifi® + + 3(dii2fiV* + ) 

+ 6(dl“8i‘lf2r3 + ) 

Of the fifty-six coefficients dp,, eight are real (Mhy^) — dm, d 222 dus? 
di 22 , di 34 , di8j, diZ 4 , d 2 j 8 — and the rest are pairs of conjugate complex 
numbers (why’) 


(ds33> 

dui). 

(ds&6, 

(dus, dm), 

(dns. dug), 

(d223i 

di24)t 

(d225i d226). 

(d|S3, di44), 

(d233, d244), 

(dsaii 

dzu), 

(d335 d446). 

(dssti d44s)i 

(d|65, digg). 

(d2i6, 

dsee), 

(dssi, d4ee)> 

(d45S, djss), 

(dgSfir dfiss). 

(di23. 

di2i), 

(di25. due). 

(dl3Si dj46)i 

(diss, dus). 

(d235( 

d246)p 

(d236. dzti). 

(d34St d346), 

{dzit d456) 


The necessary and sufficient condition that the medium be elasticallj 
insensitive, as far as 03 is concerned, to the rotation R through an angle 
0 around the a-axis is that the following relation should be an identity 
m the SIX complex vanables (fi, , fj) 

dmfi'* + + 3(dii2fi'’f2' -f- ) + Cidizsti'i'i'i's ' ) 

= dmfl^ + + 3(d„2fl*r2 + ) + C(di23f,f2f3 + > 
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EXERCISES 

1 Verify that each of the ten functions wntten above is insensitive to a rotation, 

through an arbitrary angle 6 about the a axis Htnt vi = fii *?2 + ns “ 2 f 2 
(»I2 ~ 173)* + 45J4* “ 16f3f4 1)5* + fl#* = 4fife (l|2 — 1}3)(l|5* — IJg*) — 47J41J5*|J 
reai part of and (vi — >n)9sne -h ir<(7s^ — ire^) — iniagiaarj' part of 

-Sfsfs* 

2 Show that 17273 — vi* is insensitive to a rotatioa through an arbitrary angle S 
around the a axis Hint 4(17273 — 74 *) = (72 + 73)* — I (72 — 73)* + 4 i74*l 

3 Show that 7278* + 737 &* + 2747878 is insensitive to a rotation through an 
arbitrary angle 6 around the a axis Hint 2(7275* + 7J7S* + 274787 j) = 
(72 + 78 )( 75 * + 75 *) — 1(72 — 73 ){ 75 * — 7 «*) ~ 4747578) 

4 Show that It = 71(7273 — 74*) + (7174* + 7375 * + 2747578) — (72 + 73 )( 75 * 
+ 76 *) 


It follotts from the result of Exercise 2 that ijj’is “ * 74 ^ may be used 
as a substitute for (72 — ijs)* + 4 ij 4 * and from the results of Exercises 
3 and 4 that Is may be used as a substitute for (72 — > 73 )(» 76 * ~ vs^) ~ 
4i7475r6 /i(i72i73 — 174 *) may be used as a substitute for (ij 2 + vsJivm 

— »?4*) (why’), and + » 7 $*) may be used as a substitute for 

(»72 + 178)(l75* + ’70*) (why?) Since / 1/2 - /l(l72’73 - ’74*) - i'l(’78* 
+ 176®) 4- /li7l(’72 + 173 ) = /l(’72’7i - ’74*) “ /l(l?5* + 178*) + »7I*(*72 + 
’ 73 ) + i 7 i(i ?2 + 173 )*. it follows that / 1/2 may be used as a substitute for 
’7i(’72 4 * 173 )* aiid. since /iiji* - 71 * 4 - »7i*(’72 + 173)1 that may be 
used as a substitute for i7i*(i}2 4- 173) Finally, since 7i® =171*4* 
3’7i*(i?2 + >73) + 3 i?i(i?2 + 173)* + (173 4- i73)*i /i* may be used as a sub- 
stitute for (172 4* 173)* Hence we have the following result 

A deformable medium is elastically tnsenstiive, as far as ^3(77) is con- 
cerned, to an arbitrary rotation around the a axis if, and only if, <^3(17) *5 
a linear combination of the following len functions of (iji, , 178) 

fl*» Illi , lit »7l*) 7li7l*. ’7l(’72’73 174*), /i(i72i73 — I74*)j 

171(178* + ’76*)j 71(175* + ’7S*)» 178176(172 ~ 173) + 174(175* “ ’76*) 


EXERCISE 

5 What IS the form of ^i(ij) if the medium is elastically insensitive to every 
rotation about the h-axis7 about the c axisT 


When the angle 6 of the rotation around the o-axjs to which the 
medium is elastically insensitive is the quotient of 23r by 2, 3, 4, 5, or 6, 
some of the multipliers e®**, , e~*** are umty and so some of the 

twenty three coefficients dssit » dose are arbitrary complex num- 
bers (instead of having to be zero) The medium may possess, then, 
more than the ten third order elastic constants to which it would be 


limited if 6 were different from t, ^ or ^ The results for these 
o 2*^ o 3 

vanous values are as follows (check thrae statements) 
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elastic constants, and <^3(17) is a linear combination of the ten functions 
already ^\rltten down and of 10 new functions that may be taken to 
be the real and imaginary parts of fa®, fs®, f if sfe, and f 2f3f6 (or 

any convenient ten linearly independent combinations of these) 

l\Tien ® = 2’ multipliers e*'* and e~*'' each are unity, and so the 

three complex numbers dus, d233, and d^ss are arbitrary (instead of 
being zero) The medium possesses (at most) 10 + 2(3) = 16 third- 
order elastic constants, and ^3(1;) is a linear combination of the ten 
functions of (iji, , 179) that we have already written do^ra and of 
SIX new functions that may be taken to be real and imaginary parts of 
f if 3®, f jf 3®. f 4f 5® However, the form of 0s(»?) >s most easily denved as 
follows Since the medium is elastically insensitive to a rotation 
through IT around the a axis (why?), 3^3(17) is of the form 

“ Chi’?!* + Cjsjj??* + 

+ 3(Cu8’Jl®»7a + + C344173»74*) + 6(Cl23»7l’72’j3 + 

+ C234’7J’73’J4) + 3(ci}6l7i + C2SS>72 + C365i?3 + C456’7«)’76® 

+ 6(0186171 + C2itVi + C386173 + C459l?4}’75’78 

4- 3(Ci66i7l + C266’J2 + C368T73 -f- C49674)>J8® 
Under the rotation through ^ around the a axis, namely, 
a' — a, b' — c, c' “ —b, 

17— * 17' where 171' = tju 171' = 173, vz' “ 172, v* *= “»74, W ~ ~V 6 , 
Vi' ~ Vs, and since the relation ^3(17') = ^5(17) must be an identity m 
the SIX vanables (iji, , 179) we obtain the following relations 

Cll 4 = CiS 6 = C234 = C444 — 0 , C112 = Cii3, C122 = Ci38, 

Cl24 = Ci34, Oi55 = Ci68, Cjjj = €333, C22S = C28S, 

C224 = C334, C244 = C344, C 2 S 5 = C366, C 25 G = — C 356 ) 

C266 ~ CS55. C455 = 6488 

Thus 3 <t>ziv) IS of the form 

3^3(17) = Cni»7l* + C222(’72® + ^3*) + 3{cii2J7i®(»?2 + 7)3) 

+ Cl22’7l(l72® + 173®) 4* 0223172^3(92 4- 173) + 0224174(172® “ 93®) 

4" Cl44l7l?74® 4* C244J74® (172 4- 93)} 4* 6JC12317iI7293 
4* 0:249194(92 ~ 93) 4- 02569696(92 “ 93) 4- 0456949596} 

4- 3 {ci 589 i( 95® 4* 96®) 4- 0235(9295* 4* 9396®) 

4- 0286(9296® 4- 9395*) + 0455774(175® — 96 )) 
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1. The modiRed stress matrix 

We have seen that vhen a deformable medium is m a state of equi- 
librium in a strained position the stress raatnx T is connected with the 
applied force per unit mass matrix F by means of the relation 
(div, T)* + " 0 

and with the applied force-per-umUarea matrix / by means of the 
relation 

rdS* « dSJ 

(dS‘ being the matrix element of area, and dSx the scalar element of 
area, of the surface of the medium when deformed) Since the surface 
of the medium when deformed is not given (what is given being the 
surface of the medium nhen undeformed), we introduce the 3 X 3 
matrix Ta, which is such that 

T,dS“ = TdS* 


We term this 3X3 matrix Ta the modified stress matrix Since T = 


\pa/ dn 


J* and since dS* = 


) (d*) ^ dS“ (see Exercise 1, p 17), 


Ta=J — 

dti 


Note Although the stress matrix T is syrametnc the modified stress 
matnx Ta is not, m general, symmetnc Ta ivill be symmetnc if J 
IS s 5 Tnmetnc and if the medium is isotropic, for, then, J commutes 

with 1 } = — Es) and, hence, with — 

dtj 

The relation that connects the modified stress matnx Ta with the 
applied force-per-umt-mass matrix (replacing the relation (div* T)* 
pxF = 0) IS readily found In fact, the virtual work of all the forces 
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nctiriR on any portion of the medium in any virtual deformation is 
furnished ijy the formula 

Virtual work = j dl'^ + I ( 6 .r)*/d>S’x 

- J, P„(5x)*FdF„ 4- (5x)*r« dS“ 

= f !p„(ot)*F + div„ (7’„*5x)! dF„ 

Jy. 

M 

vhere, if p "= I 7 1 1“ any 3 X 1 matrix whose dements arc dilTer- 

entifihlo funetions of {a, h, r), we understand by div^ p the sum pa + 
q,, 4- r,. Similarly, if P is any 3 X w matrix we understand by divo P 
the I X m matrix whose elements are obtained by applying the 
operation div„ to the various columns of P. 

This virtual work must be zero for the virtual translations (for which 
6t is a eonstant 3 X 1 matrix) for everj' portion V„ of the medium, and 
so 

(div. T„*y 4- PJ' = 0. 

On substituting for 7'„ its value ./ so that T„* = ( — ) J*. we 

f)n \dri/ 

obtain the basic P(}uations of equilibrium of the medium in a form 
Miitalile for employing the initial coordinates (o, b, c), rather than the 
final coordinates (j, j/, :), of a typical particle of the medium as inde- 
pendent variables; 




4 - P.F - 0 . 


When the medium is isotropic (so that ~ is symmetric^ and the 

\ J 

applied mn.s.s force is either absent or negligible, these equations 
simplify to 


or, efjuivalently, 


cUv„^.7- 

U!J 


-0, 


diVrt '■ f^- 


J. In the intiratesirnai theory of olasticitv, ~ is- replaceri bv 

■ 5rj ■ * 
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X/i + 2fir] (the initial stress being zero), and the product of any ele- 
d<f> 

meat of — by the derivative of any element of J* irith respect to c, b 

drj 

or c IS neglected The equation div^ ^ reduces, then, to 

^dtVa lA* = 0 or, equi\alently (since J* is non-singuIar), to 
div« (X/i + 2/iij) = 0 


If we assume the medium to be homogeneous, so that X and are 
constant functions of (a, b, c), this equation yields 

X(/i)a + 2/1 div„ ij = 0 

Note S In the deri\ ation of the basic equation 

(dlV„ Ta*y + PaF = 0 

we have assumed that the final reference frame is the same for every 
point of the medium (so that the virtual translations are charactenzed 
by the constancy of the 3X1 matrix ix) It sometimes happens (for 
example, when we wish to use space polar or cylindrical coordinates) 
that the directions of the axes of the final reference frame vary from 
point to point of the medium Let Rj be the rotation matrix that 
transforms any convenient fixed reference frame to the (moving) final 
reference frame If is the 3X1 matrix whose elements furnish 
the coordinates m the fixed reference frame of the vector whose 
coordinates m the moving final reference frame are furmshed by the 
elements of the 3X1 matrix 5z, we have Sx = and the virtual 

translations are charactenzed by the constancy of the 3X1 matnx if 
The basic equation is now (show this) 

{div„ + p.Ri^F = 0 

When the applied mass or body forces are zero this reduces to 
div„ (r,*Rt*) = 0 

When IV e are using orthogonal curvilinear coordinates it is convenient 
to choose the fixed Cartesian reference frame to be that determined by 
the orthogonal curvilinear coordinates at the imtial point (c, 6, c) 
Then Ri is the rotation matrix that rotates this reference frame into 
the reference frame determined fay the orthogonal curvilinear coordi- 
nates at the final point (i, y, z) In performing the operation diVo on 
the 3X3 matrix Ta*Ri*, we must use the form of the divergence 
operator which is appropriate to the orthogonal curvilinear coordinates 
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*.vft arc using (sec Introdudion io Applied Malhcviaiics, pp. 111-112’). 
Wlien (*, y, z) nrc linear functions of (a, b, c) the elements of J — 

=• arc con.stants; hence the elements of M = J*J and of >? = 
(a, h, c) 

^{^f — Ez) nrr; con.stants. Convensel}-, if the elements of -q arc con- 
stants (i.c., if the. .strain matrix is a constant matrix), the variable.s 
(x, y, z) are linear functions of the variable.s (n, b, c). In fact, the 
element .s of arc con.stant and so Exo* and SxaXs are con.stants 
(E denoting summation with respect to x, y, and z). On differentiat- 
ing the firdt of these with rc.spect to a, b, and c, we obtain SxaXaa = 
0, SxaXnt, = 0, Ex„xaf = 0. On differentiating the second with 
rc.'spect to a and using the relation SxaXd. == 0, we obtain SxfcX„o = 0. 
Similarly, 2xrXon - 0, and the three relatioms Sx„x„a = 0. 2x^.T„o = 
0, SxcXoo = 0 as.surc u.s, since the matrix ,/* is non-singular, that 
= 0, Paa ~ 0, Za„ ~ 0. Similarly, Xf.i = 0, = 0, zu, = 0 and 

Xff - 0, Pfc = 0, Zee ~ 0. IVe have already .seen that ^XaTab = 0, 
and, on interchanging a and h, this yields Extx-„f, = 0. On differentiat- 
ing the con.stant SxcX„ with re.spect to b, we obtain SxcXat, = — SxnXtc 
and, on dilTerentiating the constant Sx^Xe with respect to a, wc obtain 

SxoXj, with 

SXeX„h = 0. 

The three relations SxoXai, = 0, — 0, 2xcX„i, = 0 a.s.sure us, 

since the matrix J* i.s non-.singular, that x„j, = 0, yah = 0, z„i — 0. 
Similarly, xj,, = 0, yi,^ ^ 0, ztr = 0, and x^,, = 0, t/m = 0, = 0. 

Hence .all the second derivatives of (x, y, z) with respect to (a, b, c) arc 
rero and so the variables (x, y, z) arc linear funclion.s of (a, b, c). Wc 
say that a strain is /lemoyr/moi/.'; when the strain matrix q i.s constant. 
We have, then, the following result; 

A ftrain is hmnoyrnrnus u-hr.n, and ottly irhcn, the Jacobian matrix 

J « ^±'bjr) j-e fj coustnnl matrix. When this f.s so the vanable.'} (x, u, z) 
(a, b, c) y jr / 

arc /sorar/fnidtoa,'? of the variables {a, b, c). 




Sx{Xfa. On differentiating the constant 

J V , 


respect to c, we obtain 2^XaXu + ^xipr^a — 0, and so 


\\ hen the strain i.s homogeneous 




J* Is a constant matrix and 


so we have the following resjilt: 

When the ftrain is hotnoyenrous the applied force per 7 init mass is 
zero (m J.cd mhen the applied force per unit mass is not zero the strain 
cennot fs* honngrneous). 

'Hie following .*<ection-S di.scu.s,s two important e.xamples of homo- 
geneuu-'' strain. 


* By rrs.swj' J). Mtimnrh.sn, John tViJ,-y A 
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2. Simple shear 

In simple shear the deformation is described by the equations 
X = a + kb, y = h, z = c, 


thus each particle is displaced, parallel to the x a\is, an amount 
proportional to its distance from the 6-plane We have seen (Exam- 
ple 1, p 33) that the principal axes of strain are obtained by rotating 


the (a, 6) -axes around the c-axis through an angle a 


1 _ 
■ j + jtan 


li 

2 


The characteristic numbers of M are (14-^ sin 2a 4- sin* a, 
1 — A. am 2a -h A* cos® a, 1), and if we use the relation fe = —2 cot 2a 
these turn out to be (tan* a, cot* a, I) Hence the characteristic 
numbers of are (tan a, cot a, I), and so 


( cos a — sin a 0\ /tan a 0 0\ / cos a sm a 0\ 

sin a cos a 0 )[ 0 cot a 0 ) I — sm a cos a 0 I 

0 0 I/\ 0 0 l/\ 0 0 1/ 


( sir 
- cos 

0 


sm 2a 
cos 2a 


— cos 2a 

2 cosec 2a — sm 2a 

0 


:) 


’ (4 + l*)- 


(2 k 0 \ 

r»(i 2 + k^ 0 

\0 0 (4 + );*)7 


Hence 4/”” =■ (4 + S’) 


/2 + 0 -k 0 \ 

-« { -/t 2 0 ) 

\ 0 0 (4 + 


and BO the 


rotation matrix R = JM ^ is given by the formula 


( 2k ® \ A*" 

-12 0 )= c03 2o Bin 2a 0| 

0 0 (i + kY/ \ 0 0 1 / 


Hence f? is a rotation through ^ — 2a = —ton ’ - around the c axis 

Tlius the coordinates of the stress tensor T, m the reference frame 
obtained bj rotating the onginal (x, y, z)-axca around the z-axis 

through the angle — tan'’ are the elements of the matnx 
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The coordinates of T in the oripina! (x, y, r)-referencc frame are the 
elements of the matrix 



Since del .1 ~ Pj ~ Pa (so that no compression is involved or, 

(‘{inivalently, voliiincs are prcscr\'cd in a simple shear). The modified 

. . , d<i> 

stress matnx is = ./ • 

dr, 

We first consider the situation in which the medium is isotropic, 
and we take the initial stress to he a hydrostatic pressure p. Then 

tr- — 7 j/j + — j _ Ofili -r — — 2mliln-i- nl;^, 

2d o 

-}- {Kl \Ez -}* 2p.T)) -f- {IJ i~ ~ 2ml^hz + 2ml\ri n co y 
dn 

where 


/O h 0\ /O 0 0 \ 

^ ( /: /r 0 , CO = 0 0 0 )• 

' \0 0 0 / \0 0 -\}rl 

Since we have stojiped the development of at terms of the third 
onler in the. elements of y, we arc entitled only to keep terms not 

involvinp h to a hipher power than the second in the expression J — 

dy 

Fince /i - Z; —V*'’ have, to this order of approximation. 


dy 



1 


(X 4* 


Ez + 2fty y CO y, 


Ta - J 


£6 

dy 


+ 


-pEz 

f^iX 4* "T vi)k' (ft — p)h 


0 


0 

y (X 4* 2fi 4* T>i)h~ 0 

0 i{X 4- m - hi)h-, 


The force on a matrix element of area dSr of the boundary of the 

/ r 0 ON 

iiefomsed nu'dium is given by TodS'\ Since -(—/.• 1 0 

\ 0 0 E 

'•e have »■' ttS’, v- dS*' — f: dS^, dS" ~ dS'^, An element of 
till A 'A Inch lay initially in the o-planc (i.e., for vhich dS'^ ~ 0, dS^ = 0) 
O' tleiornuai into an element of area 0) who^e tnapnitude is 
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the product of the magnitude of dS“ bj (1 + The force, per 

unit initial area, on an element of area^hich laj initially in, or parallel 
to, the a plane is found, then, by multiplying the first column of Ta by 
(1 -j. = 1 — Thus the coordinates of the force (per unit 

initial area) that acts upon an element of area ^\hlch was initially m, 
or parallel to, the a-planc are (— p + ^(X + 2/; + m + p)fc*, nk, 0) 
Similarh , the coordinates of the force (per unit initial area) that acts 
upon an element of area nhich was onginally in, or parallel to, the 
6-plane are ((a» — p)k, — p + -^(X + 2fi m)A.®, 0), and the coordi 
nates of the force (per unit initial area) that acts upon an element of 
area n hich n as ongmally in or parallel to, the c-plane are (0, 0, ~p + 
^(X + m — ^n)l}) If p = 0 the first-order part of this force per 
unit area is furnished b> the matrix 

/O M 0\ 

A.(m 0 0) 

\0 0 0 / 

This IS the result of the linear or infinitesimal theory of elasticity The 
second-order correction that must be added to this is furnished by the 
matrix 

( X -|- 2;i -}- m 0 0 

0 X + 2/1 + m 0 

0 0 X + m - -Jn 

(0 0 0 \ 

“ i’fX + 2/1 -f- m)k^Ez — (/< "f* ® ® ® 1 

\o 0 1/ 

Thus in order to produce the simple shear a hydrostatic tension of 
magnitude -^(X -H 2/i + m)k^ must be applied to the medium and, m 
addition, a simple linear compression ( «= negati\ e tension) of magw 
tude (/I -f -Jn)A.* parallel to the z-axis, these forces being added to the 
shearing force furnished bj the linear theory If these surface forces 
nhich are necessary, in addition to the shear forces furnished by the 
linear theory, to mam tarn the simple shear are not supplied, themediiuD 
will contract (because of the absence of the hydrostatic tension) , fut 
thermore, it will expand m the direction of the z axis and contract la 
anj direction perpendicular to this, because of the absence of thesimple 
linear compression m the direction of the z axis (see the followng sec- 
tion, m which simple tension is treated) Since both the hydrostatic 
pressure ( = absence of hydrostatic tension) and the simple linear tension 
(= absence of simple linear compression) are second-order effects, the 
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mnKiiiludp of oacli having k' ns a factor, we may use the linear theory 
in appraising the effects duo to them. Thus the dcfonnation due to the 
hydrostatic pres-mre is furni'^hcd 1 ) 3 ' the formulas x' = (1 — ^)i, y' 

- {1 - a'ni. « {1 - ^)r; ^ 

11(0 next sect ion that the deformation due to the simple linear tension is 


y + , ■ 


J' , 


V 


1 


/r j y\ 






where c 


. „ . , , r, p(3X + 2\l) . 

IS Poisson > ratio and h — — ^ is Young s 


2{X -p X + 

modulus. On romhining tiiesc deformations with the simple .shearx = 
a -f Lb, y h, z r-. r, we obtain (on neglecting poiver.s of k higher than 
tiie s''cond) 

jr" = (1 - a)r)a + kh, y" = (1 - ak-)b, z" = (1 + yk-)c 

+ m 
2(3X +~^’ 

In ortier to vcrif.v the correctness of this solution we have morel 3 ' to 
start with the deformation 


X 4- 2g + TTi r / 1 \ + -J-n X + 2p 

where c* i — “lu-j — nl. v = — — 

2(3X4- 2a) ^ EV ' -1 / ^ E 2(3X4 


X (1 — ak’)a 4- kb, 


V 


a - ak")b, c = (I 4- yk^‘)c 

whore « and y are undetermined constants. \Vc obtain (alwa 3 ’’s 
neglecting all powers of k higher than the second) 

/l - nfr k 0 

./ - ( 0 1 - air 0 

\ 0 0 14- ylrj 

h ~ 2ak- k 0 

.M - I k 1 4- (1 - 2a)lr 0 

\ 0 0 14- 2-,k-j 

~iJr U 0 \ 

ik (I ~ a)lr 0 J. /, (i- - 2o 4- y)lr, h - -\k-, 

ylrj 


prs 


(1 


0 


0“ 


(i " + y)’r~ m 


k'Ei 


■2o;.-' /; 0 

(1 - 2a)lr 0 
0 0 2yk- 


<0 0 (A 


--n/.- 0 0 0 
\0 0 L 
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QA /{^d ~ 2 a + 7 ) + -Jm + ;j(l — 2 a)}fc^ pfc 0 \ 

J— =1 pX |X(i - 2a + 7 ) 4- im + p(l - 2a)}fc® 0 ) 

\ 0 0 {X(i-2a+7)+im + 2M7-in]JL7 

/O pk 0\ 

In order that J should reduce to ( 0 0 I, the undetermined 

\o 0 0/ 

constants a and 7 must satisfy the two equations 


2(X + p)a — X 7 = -^X + + p 

2Xa — (X + 2p)y = -^X + “ "iti, 


and a simple calculation (perform it) shous that a and 7 have the 
values given above 

The values of a and 7 given above correspond to an applied force per 
unit area on the planes, into which planes parallel to the o-plane are 
deformed, of amount pk parallel to the y axis and an applied force per 
unit area on the planes, into which planes parallel to the b'plane are 
deformed, of amount pk parallel to the x^axis The second of these 
applied forces is a shearing force since the orientation of any plane 
parallel to the h-plane is not changed by the deformation (any such 
plane being deformed into a plane parallel to it), but the first is not 
Any plane element of area of the a plane is deformed into an element of 
area i\hich is perpendicular to the vector r(l + (7 — — fc, 0), 

for, on setting a *= constant m the equations defining the deformation, 
N\e obtain dx « kdb, dy ^ (1 — ak^)db, rfs = (1 4* and so 

dS* dS® dS* =* 1 + (7 — a))l* — fc 0 (powers of k above the second 
being methodically neglected) If we wish, then, the forces on those 
planes, into which planes parallel to the o-plane are deformed, to be 
shearing forces, the first column vector of Ta must be perpendicular to 
the vector t)(l 4- (7 — ce)k^, —k, 0) We start out with the Inal 
deformation 

X = (1 — afc*)a + A.f>, y = {I — pk^)b, z = (1 4* 'Yk^)c 


where a, 0, and 7 arc undetermined constants We find, as before, 


M 


( 1 - al* k 0 \ 

0 I - 0k^ 0 ), 

0 0 1 + 7^7 

( 1 - 2afc® k 0 \ 

k 1 + (I - 0 

0 0 1 + 2yky 
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/-air 

u 

o\ 


V 

Ik 

('i - 0r 

0 Ui 

= (i — a — /3 4- y)k‘, h - 


\ 0 

0 ylr/ 


Or, 

. [x(|. 

^ 5 1 

-a~() + yi + - 777 1 

IrE, 4- 


■ilr. 



/~2a}r k 0 \ j /O 0 0\ 

M /,- (1-2^)/" 0 --7i/r 0 0 0 . 

\ 0 0 2ylr/ \0 0 1/ 

/}X(i — a ~ /3 + ■>■) "f + #7(1 ~ 2a)]k~ fih 0 \ 

I fik |X('i- — a — + 'i') ■}■ ^771 4* p(l — 2^)}^" 0 j. 
\ 0 0 wo, - a- 0 + y) +h» + 2^y\}r -in/ 


Tlio, thn‘e oqiijilions that .‘•crvo to determine «, (S, and y are 

(X + 2/1)0 + X/5 - X7 = l\ 4- ini. 

X« 4* (X 4" 2a)/? — X'/ = -i-X 4" 4* Mj 

X« 4- X^ - (X 4- 2ti)y ^ 4/X 4- ^771 - -J-77. 

I'rom the finst t^vo of thc-sc \vc obtain /3 = a 4- {>•, and, on .snl)stituting 
thi•^ Value in tlie ficeond and third, we ohtjiiji 


2(X 4" a)« “ 

2X« — (X 4- 2a)T ~ — j-a, 

and ao 

^ m ^ 7iX __ 711 ii{\ 4- a) 

2{rjX 4- 2a) '' si^x"T'Ta)’ 2(3X + 2a) 

With thase values of «, /3, y, T„ — J — is 

Oij 

fulr a/.' 0\ 

T„ - a/‘ 0 0 

\ 0 0 0/ 


and the Fh(‘nring force per unit area on phincs into wliieh jdanes 
!>:tr,n!lel {tj the fi-planc and to the h-p!ane are deformed is nk. 'I'ho 
d'hnaiation is 


n^-')f7 4- kl>, II - II - (o 4- i)k-]k, 


(1 -r yfr)c. 


and the eotnprosion ratio -- is funushed bv tlu! formula 

Po 


idei ./)■ 




(>77! — n I 

+ 2i7)| ' 



106 


SIMPLE SHEAR AND TENSION 


The stress matrix T is 




nk 0 \ 
0 0 
0 0 / 


If the initial hydrostatic pressure p 0, the values of the elastic con 
stants (X, tx, I, m, n) are different since these are functions of p The 
equations determining a, 0, and y are also slightly modified since 

— and, hence, Ta — J~ involve p We have (show this) 

dy] ojj 


-P + (X(i — a — /3 + y) + (fi — p)A: 

4- /i(l ~ 2a) 4* pa\k^ 

nk ~p + |X(-^ — a — ^ 4- v) 4- 4- ft(\ — 2/3) 0 

0 0 — p 4- {X(^ — a — 0 4- 7) 4- 4- 2 m7 ^ 

-in — py]k^ 


An element of area of magnitude dSa of the a plane is transformed into 
the matrix element of area 



whose magnitude is dS^ (1 4* (i- 4- 7 — Hence the resolved 

part of the force per unit area on dS*, normal to it, is (show this) the 
product of 


~P 4* {X(i - a - /3 4- 7 ) 4- - 2aii 4- p{« 4* )3 - 7 )jfc* 

hy I — Ci 4- 7 ~ le ,—p+ {X(^ — a ~ |3 4- 7 ) 4- iw - 2af( 4 - 
p[^ 4- and if this is to be — p (so that the additional applied 

force on the matrix element of area dS* is a shearing force) the unde- 
termined constants a, 0, y must satisfy the relation 


(X 4- 2 ^ — p)a 4- XP — X 7 “ -J-X 4- ini 4 - ip 

An element of area of magnitude dS^ of the b plane is transformed into 
the matrix element of area 

/ ® 

dS* “ dSa I I 4" (7 ” a) 

\ 0 

and so the resolved part of the force per unit area on dS*, normal to it, 
13 (show this) 

— P 4 - {X(^ — a — ^ 4- 7 ) 4* im 4- m(1 — 2^) 4* p{0 4- 7 — 
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Ht-nro, in order that tlic ndditiona! applied force on (IS‘ be a shearing 
forre, the undetermined constants «, 0, y must satisfy the relation 

(X + p)rt + {X -p 2 a — p)0 — (X -f p)y === -j^X + + a- 

Similarly (show this), on considering an element- of area of the c-plane, 
we obtain the relation 

(X -r 7>)o + (X + 7>)0 — (X + 2 a — p)y ~ -UX + jrm — -^n. 

When rr, 0, y arc a-ssigned the valuers that satisfy these equations, we 
have 

f~P + (p - (a - pV: 0 \ 

7 ‘„ “ I ft): - p + 7)(o - y)):- 0 j, 

\ 0 0 -p + pUt + 0)):-/ 


and, since -- r- i — (7 — n — 0 )):~. 

Pn 

/ - 7) + 1 2a -f p{y - ^ 1 k- {ft - p)k 0\ 

= (a - P)k -p 0 . 

\ 0 0 -pj 

Thus the initial pres>un* a(Tc(:t,s even the lirsUordcr approximation, the 
efff'Cl being an ajiparent reduction of the modulus of rigidity a by p. 
.Since, however, a in itself a function of p all we caji say is that p — p 
jdays, when p 9^ 0, the role played by ao 'vlion p = 0, ao being the 
v.'ihie of a when ;> ~ 0. 

On adding together the second and third of the three equations that 
K'r\'e to determine n, 0, and y, avo obtain 



(X + p)ft d- (X -r f‘)(0 ~ y) ~ ■^•X +■ ^7?! -f- dra ~ -J-m. 

and. on combining this with the first oi the throe equations, eve obtain 

.7 „ d- ?»Ha — p) 4 - {^p — l n) i\ d- 2a — p)_ — ^ ^(X -f p) 

3X^ *5* 2/1* p(2X 4 * m) 

jup d- A-p(X -f a) 

.“iXa d- 2a' - pC^TTaT’ 

f)n subtracting the third of our thn'e equations from the .second, wc 
obtain 

, a?' d- 

0^-7 

p - r 
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and these formulas determine a, P, and y, when X, n, m, and n are 
knoun as functions of p Conversely, a knowledge (obtained hj 
experiment) of how the compression ratio, for example, vanes with 
p furnishes information as to the nature of the dependence on p of the 
function ()3 — 7) + a of X, n, m, n, and p 

The theory is more complicated when the medium that is being 
subjected to the simple shear is, like wood, non-isotropic We con 
sider first a medium that is elastically insensitive to any rotation 
around the a-axis (the direction of the shear) There are, then, three 
additional second-order elastic constants and seven additional third 
order elastic constants (the medium having fifteen elastic constants m 
all, five of the second order and ten of the third order) However, 
since Ii and many of the dements of ij are infinitesimals of the second 
order and not, as is generally the situation, of the first order, not all 
these ten additional elastic constants will affect our second order 
approximation Adopting the notation of Exercise 12 , p 95 , where 
the three additional second order elastic constants are denoted by 
(Si, $2, 63) and the seven additional third order clastic constants by 

(^41 » 8io)» the addition to the matrix — due to the non-isotropy 

drj 

of the medium is 


/hvi + ^2(^1 + ’ll) 0 
I 0 isvi 4 * hvj 

\ 0 0 


/~i\a + — 2o — ^ + 7) 

fcM 0 

\ 0 


: ). 

0 0 

“^20 + {57 0 

0 “ 52 “ + — ^), 


(Remember that, when evaluating — » ^ must be \vritten symmetrically 
dri 

and so 71273 — 74*, for example must be replaced by ybbVce — VbeVeb) 
Thus the seven additional third-order elastic constants do not affect 
our second order approMmation (only the three additional second- 
order constants ( 5 i, S2, 5 s) being involved) The addition to Ta is 

the same as the addition to ~ (since the addition to — is of the second 
ay dy 

order in k)j and the addition to tiie stress matrix T is also, for the same 

reason, the same as the addition to ~ The equations that serve to 

dy 

determine the constants a, P, and 7 are (the initial stress being taken 
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to be zero) 

fx -f 2;j + 5i + 2o:)« -f (X -f- ojlfl — (X -f 5 ;)y = iX -f i-m -f 

{X *f 4- (X -t- 2;i)/3 — (X -i- 03)7 = -^X -j- hn -p n. 

(X h)n + fx 4- h)fi - fX -P 2p)-y = LX + Im - Jn + ios- 

Wlicn tlio values of n. and 7 wln»'h satisfy these equations are used, 

/nh- 0 \ /2pJr fik 0 \ 

T„ - UL 0 0 . T = ;d: 0 0 , 

\ 0 00/ \ 0 0 0/ 

and 

~ = 1 -f (a -p j8 — t)!*"' 

Pa 

XVlion the. inediurn is insensitive (0 any rotation around the 6-axis, 
the additional terms in arc of the form 


4* Pih^i 4' ” ’ll.') 4- -h 

4 - ~ VtT) 4- PihiViVl ~ Vh‘) 4* 5s7?;(7;c' 4* V^) 

-p l>)h(.ve.‘ 'P nU) 4* Oloi’t-j’jcf’;.'! •” vi) 4- »7r>(nc' — ’?•(■))> 


and so the additional terms in — arc 

Oij 


1 • 2 


(h^riz + 05153 -f* 05 TJg“ 0 

0 ^ 1^2 "**■ ^ 2 ^^ 1 "f* ^ 2 ) 'f* OgT/r," 0 

» 2 
•r 05>;f. 

0 ^ 2^2 “ 1 “ 03:71 4" 05 : 7 ^“ 

^ . 


/P:(A — ^) 4* 5j7 4- 0 *^'5 10 


14 1 


0 


0 


oi(-^ — d) -p 62(1 — « — 2 ,S -p 7) 

4" ’i'O'i 4" 4 05 

•^SlO 0 02(4 ~ 0 ) — 5 ;rt -p -^Oil 


The additional third-order elastic constants o^, 65 , and Sm are now 
inVflved in our second-order apitroximation. Furthermore, there is 
Ui'-.v a shearin? 2 : force on plane.s into which planes parallel to the c-planc 
fire dtfoniKsl (this shcarint: force beinp; J-SicL' per unit area parallel 

Irt the JT-axis). 

If the nutlsum F elnstically insensitive to any rotation around the 
•• » 

r-rAp, the additional terms in — (over and abo%'e the terms that would 

{«. } 5 :x-‘ nt if the mtHliitm were icotropic), are furnished by the matrix 

ftdif'-a thi-c't 
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( 52 ^ 3+53772 —53776 

— 53776 53773 -f 53771 

0 0 


/527 + 53(-|- — (5) 


0 \ /o 53 0\ 

0 = -1.(53 0 0 

5rt3 + *2(^1 + 773)/ \o 0 0/ 
0 0 

52 Y — 53a 0 

0 ill + SiCi - a - S + 27)y 


Thus the presence of the additional second-order elastic constant 53 
affects the first-order or linear approximation, the effect being an appar- 
ent reduction of the modvlus of rtgiditt/ m by S$ We leave to the reader 
the task of writing out the equations that serve to determine the 
constants a, 0 , and 7 in this, and the preceding, case and also the task 
of developing the appropriate theory when the medium is under an 
initial hydrostatic pressure 


3. Simple tension; the linear theory 

The problem of simple tension is that of the deformation of a homo- 
geneous cylinder under the application of equal and opposite forces 
applied unifomilj over the two ends of the cylinder We take the 
c-a\is parallel to the generators of the cylinder and consider Bret a 
cylinder that is elastically isotropic and jmtially unstressed We 
start off w ith the trial deformation 

z = (1 — ak)a, y «» (1 — ak)b, 2 = (1 + k)c 

where k is an undetermined constant which is so small that wo neglect 
any power of it higher than the first <7 is a second undetermined 
constant that is, as we shall see a function of the elastic constants of 
the medium (being independent of the applied force) Our trial 
deformation is such that J, and hence 17, are diagonal, the three diag- 
onal elements of 77 are 

771 = —ffk, 172 = 773 = i 

and 80 7 i = (1 - 2 <r)k ^ = X/1E3 + 2 nv is diagonal its first and 

017 

second diagonal elements having the common value {X(l — 2cr) — 
Since we are neglecting powers of k higher than the first, the 

matn\ Ta = J — is the same as the matrix — > and, since there is no 
077 dr) 

applied force on. the sides of the cjlmder, v must be such that the first 
(and second) diagonal elements of Ta are zero (for the matrix element 
of area has its third coordinate zero on the sides of the cylinder and 
thus at least one of its first two coordinates is different from zero) 
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Hcncc 


X 

2(X + p)' 


c i« known ns Poh^on’f ratio; for many media its value is about 0.3. 
If I in the (initial) Icnplh of the cylinder the two ends of the cylinder 
ntny h<‘ taken in ho f^iven, before the deformation takes place, by the 
equations c ~ 0, r ~ 1. If the magnitude of the applied force per 
unit (initial) area on the deformed position of the end c = I i.s denoted 
bv /% we have (since the third diagonal element of 7'a is X/j 2 pr ;3 == 
ix(l - 2^) + 2,.)/: = 2^(1 + a)k), 2^i(l + a)k = P and .so k = 
P 

■ — We denote 2pt(l c) bv E and wo term E the Young's 

2ji(l 4- tr) 

modulus of the medium; 


E 


+cr) ^ 


n(3X 4- 2/j) 
X 4- p 


Tlien k - —• Thus our as.sumption that k is .so small that its square 

may be neglectcsl rnean.s that P is .so .small in comparison with E that 
p 

the squ.'irc of ~ may be neglected. Since dx dy = (1 — ah)’ da db ~ 

! 1 “ 2ak) da dh, any element of area of a cross section of the cylinder is 
decre.a.scd in magnitude by the factor 1 — 2ak', if, then, we denote by 
P' the magnitude of the stress on the end of the cylinder (i.c., the 
nirignitude of the applied force per \init final area), we have P'(l — 
2«/.') *-• P, and so 

P /^U - 2<Tk) 

■' " E 


On multiplying both .sides of this equation by 1 4- 2ffk and neglecting 
k\ we obtiiin 


k 


r 

e' 


In luirds, (o the degree of approximation furnished by the linear theory 
it m.nkcs no ditlcrencc whether P denotes the applied force per unit 
initnd area or the applied force per unit final area. In either event 
p 

k ^ y,' .and the defonnation L furnished by the formula.« 


_2' -js 





r. 



112 


SIMPLE SHEAR AND TENSION 


The relatn e extension ^ ^ 13 equal to ^ The graph of P, as a 

function of the relative extension, 13 a straight line through the origin 
having slope E 

If the initial stress, instead of being xero, is a hj drostatic pressure p, 

the first diagonal element of — is —p -f- |A(1 — 2<t) — 2/itr}i and so 
3ij 

the first diagonal element of To = J — is — p + {A(I — 2c) — 2 imt -f 
3ij 

pc\K An element of area which was onginall> normal to the a axis 
remains normal to the a-axis after the deformation but its magnitude is 
increased by the factor 1 + A.(l — a) The applied force on such an 
clement being a pressure p, we must have — p + {X(l — 2a) — 2iia -4- 
pa]k = — p{l + k{l — ff-)} andsoXfl — 2£r) — 2fic + p “ 0 Hence 


a 


X-hp 
2(X + ;.) 


Thus the elTcct of the initial pressure is to change the formula furnish- 
ing Poisson’s ratio to the extent of replacing the numerator of the 
fraction by X 4- p Of course, X and n are themselves functions of p 
80 that c depends on p not only explicitly (through the term p in the 
numerator X + p of the fraction) but also implicitl> through X and m 

The third diagonal element of ^ is -p-h |X(1 — 2(y) 4- 
dn 

and so the third diagonal element of T® is — p -h {X(l — 2a) + 2ji — 
pjk An element of area which was originally normal to the c-axis 
remains normal to the r-a\is but its magnitude is decreased bj the 
factor 1 — 2ak If, then, the applied force is of magnitude P per unit 
initial area (in addition to a pressure p), we must have 


“P + {X(l - 2<r) -h 2p - pjit = P - p(l - 2 <tA.), 

and 80 

\\{\ - 2<r) 4 - 2^» - p - 2p<r)t = P 
or, equivalcntlj , 
p 

I *= — where P = X(1 — 2a) 4* 2^ — p — 2p«y 
L 


Aswhenp = 0, it docs not matter whetherPdenotes the magnitude of 
the applied force per unit initial area or the magnitude of the applied 
force per unit final area (since we arc m^Iecling /.’) Since X(1 — 
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&■) 2n<? ~ p wc have 


E - 2{p - p)(l + <r) 


(ft ~ P)(3X -f -f p) 

X + ft 


Tinn the cEcct of the initial pressure, in the foiTnnlas furnishing <r 
an'l E, is to increase X by p and to decrease by p (so that the denom- 
inators of the fractions that furnish o- and E arc unaffected forinalh'; 
they are, of course, affected by the fact that X ft is a function of 7)). , 
If the initial stress is a non-scalar stress, we know that the medium 
cannot be elasticalb’ isotropic (why?). Let us consider an initial 
stre-'s matrix all of whose elements arc zero save the one in the third 
row and third column, and let us denote this non-zero clement by P. 
We assume, further, that the medium is elastically insensitive to any 
rotation around the c-axis. Then <i(jj) is of the form 


^iv) ~ E'Js d — ~ h~ ~ 2ft/c -j- Oil i??.-! -f 63(771772 — 77 r,*) 


and so, .«ince 77 is, bv hvpothe.sL, a diagonal matrix, is a diagonal 

di? 


matrix whose diagonal elements are 


Xf 1 -f 2a»7i + 62773 -b 63772, Xl 1 -b 274772 + 60773 -b 63771, P XI i 

-b 2/1773 -b 61773 -b di{Ii -b 773), 


i.e., 


|X(l — 2 <t) — 2 /i <7 -b 62 — 53crj/;, {X(l — 2<r) — 2 /j(t -b 62 — 65<rj/,‘, 

P -b 1X(1 - 27 t) -{- 2/1 -b 61 -b 262(1 - 


lienee, .since the fir,^t and second diagonal elements of T„ 
be zero (why?), a is determinoil by the formula 


- J 


t)(^> 

“ must 
dv 


X "b 60 

2X -b 2/1 63 

The (sjuution determining k is 


IXil - 2c-) d- 2;i -b 6, -b 2o;(l ~ c) d- P\h - AP 

r4ic-re P -b AP is the applied force per unit iuitijil urea. The cocflicient 
of k i-i YoungV mcKiulus E and, since X( 1 — 2<rj -- O/icr — 62 -b 6:0', 
wclu-o*' 

/: - 2.11(1 -r <7) -b 6s -b 6;(1 - 2 <t) -b 6;« + P. 
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Since k IS the relati\e extension ^ we have 


z — c _ ^ 

c “ “e 

where E is a function of P both explicitl> and through the elastic con- 
stants fi, a, 5i, Sj, and 83 On denoting z — ch) Ac, we obtain, on 

letting AP —* 0, from the approximate equation ~ = — the exact 

c E 

equation 

tfc _ dP 

7 ““e' 


and, on integrating this, 
loi 


(-f)=rf 


where I is the length of the cj Under when there is no applied force and 
Al is the extension when the applied force per unit final area is P' 
Note P' is the applied force per unit final area, and not the applied 
force per unit initial area, since P is the applied force per unit area, 
the area being evaluated w hen the c> Under is under the tension P 


In order to perform the integration 


Jo B’ 


we must know how the 


Young's modulus E of the medium vanes with P If wo make the 
drastic assumption that n and a are independent of P and neglect the 
lack of isotropy introduced bj the application of the tension P, E 
is the linear function 2p(l -h (r) + P of P and we obtain 


log(l+y)=loK(l+0 

or, cqui\ aicntly, 

I Eo 


where Eo = 2/i(l + e) is the Young's modulus of the medium when 
unstressed Thus, under these unrealistic hj-pothcscs, the cxtcndwl 
linear theory furnishes the same result as the simple linear theoiy, 
namely 

The relatuc extension is proportional to the npplie<l stress (or force 
per unit final are^' 

This result ike’s law in its simplest form The UirimiK 1 

( 
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j* mav be regarded as the general form of statement of 

\pj <)n 
Hooke’s law. 

If we suppose that E may be sufficiently approximated, over the 
iiiteri'.'U (0, 1^’), by the linear function Eq -f ctP, where a is a constant 
that is def<Tmincd by experiment, we obtain 

f) ' ^ i) 

or, equivalently, 



If n happens to be negative this formula indicates that as 

E 

P' ---• In fact, on setting a - — /3, we obtain 
— a 


I 


0 


1 


- 0 



- L 


.and, as P' — * —• ( 1 — /3 — * 0. Thus the extended linear thcorj' 

is capable of explaining the phenomenon of flow of a cylinder under 
Icu'-ion; but. in order for it to do so one must take into account the 
de|?en(lcncc of the elastic constants on the applied tension. It scem.s 
hardly likely that the mere dependence of /i and e on P would be 
fufficient to account for this phenomenon of flow; it i.s to be expected 
that the elastic eonstant.s oj, So. oj (which measure, by their differenc&s 
from rcro, the lack of isotropy of the medium) must be considered. 


4. Simple tension; the .“icconrl-ordcr njijiroxinmlion 

Wc consider ftr.-:t a medium, .supjroscd to be initially unstressed, that i.s 
i^!)tropic. Wc assnmc the trial deformation x — (1 ~ ah tih')a, 

P 

y (I — eh -f i^h~)h. r — (I -{- !• 4- 5/;')c where /.* = ” and where a. 

Jj 

E .are the ehcitic ron.stants (Poisson’s ratio and Young’s modulus) of 
t!ie •Imph* liitear theoty*: & and o are two new undetermined constants. 
V, r raglect metln'Klicnlly powc^^ of h above the second, rj i.s a diagonal 
matrix wh.O'-e di.ngonal elements are 

— eh -t- (5 d" Ee"}},"" = t;i 
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and 80 

/i = (1 - 2a)k + (20 + 5 + ^ 

Is = (ff® - 

The constant a is such that the coefficient of k m the first diagonal 
element of X/iEj + 2;xij la zero and so the first diagonal element of 

^ 13 the product of A,* by X(20 + S + i + 2;i(0 + V®) + 
O’? 

1(1 — 2ff)* + 2mff(l <r) — na Since we are neglecting poAvers of 
k above the second, this expression also furnishes the first diagonal 
element of Ta (nhy?) and bo the two undetermined constants 0 and 5 
must satisfy the following relation (why?) 

2CX + /t)0 + = -X(i + ff*) - - 1(1 - 

— 2ffiff(l 'bo)*}' no 

The constant E is such that the linear part of the third diagonal ele- 
ment of To IS P, the applied force per unit initial area Since the third 
aa 

diagonal element of ^ is {X(l — 2cr) -j- 2n]k + {X(20 4* J + 
ff“) + it(2B + 1) + 1(1 - 2.r)“ + 2m(l - <t’^) + rur^k^ the third 
diagonal element of T® “ J is P + iX(20 — 2# + ^®)-}- 

/i(25 + 3) + f(l — 2<r)® + 2m(l — »*) + On setting the 

coefficient of in this expression equal to zero, we see that the unde- 
termined constants 0 and 5 roust satisfy the relation 
2X0 + (X -b 2^)5 * -Xd - 2<7 -b «r*) - 3 m - 1(1 - 2(t)2 

— 2m(l — (T*) — TUr®, 

and, on solving for 0 and 5 the two linear equations that they satisfy, 
we find 

X(3X -b 4 m) m‘^1 X(3X= + 6 Xm + 4m=) 

8(X + fr (3X + 2u)(x + rt* 8^(3X + 2rt(X + ri’ " 

j ^ _3 (3X + 2ri 

2 (3X + 2ri(X + ii)‘ 2(X + m)' “ 



4,.(3X + 2rt(X + " 

(a “ c) P 

The graph of the relative extension ^ plotted against ^ ^ 

IS a parabola passing through the origin with slope 1 and second 
den\ ative 25 
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The d(‘pr;n{lcncc of Younp’p modulus E on the applied tension, 
areordinp to the second-order approximation, is obtained as follows: 

(s — c) P 

On writing the relative extension in the form —> we have 

c I'j 


!L 

E' 


- k -f 5}r 



and so 



or, ecjuivalently, 



Thus the constant a of the extended linear theory is —6 if the second- 
order approximation is adopted. In order that the .second-order 
approximation may furni.sh (through the method of the extended linear 
theory) an explanation of the phenomenon of flow, it is necc.ssary that 
f he po*-ilive; this being the ease, the medium will flow before the 


applie<l force per unit area attains the value 


Eo 


To examine the ctTcct of an initi.al Itydrostatic pressure p (following 
the theory of the second-order approximation), we ob.serx’c that the 


tir.-t diagonal element of — is 

dr? 

"P -r |X(l — 2<r) — -f- !X(2/J 5 -b -i- -f < t ‘) -j- /i(2/3 

d- <r") -i- /(I - 2<r)" -f 2m«r(l -j- .r) - 

•» jt 

and s:o the finst diagonal element of 7 « = J — is 

t)’? 
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the coefficients of we obtain the equation 

2(X + p),3 + (X + p)a = -X(i -a-i- 3tr®) - - 2<rf 

— 2ffl<T(l + er) 4" ««■ + pff 

connecting the undetennined constants 0 and 6 On using the relation 
p = 2(X + n)a — X we may write the right hand side of this equation 
m the form + <r)* ~ fia* — 1(1 — 2(r)® — 277iff(l + (t) + rw 

To obtain a second equation connecting 0 and 5 we observe that the 

third diagonal element of — is 
Bit 

-p + |X(1 - 2.7) + 2;.1J: + (X(2fl + S + i + O + p(2S + 1) 

+ ((1 - 2.7)’ + 277.(1 - 

The third diagonal element of T* is the product of this expression 
by 1 -h k + 5k^ i e 

~p + fX(l - 2cr) + 2m- Pik + tH20 + S + f - 2cr + O 

+ m(2S + 3) + f(l - 2ff)* + 2m(l - <t*) + - p^lfc* 

and this must be equal to the sum of P (the applied force per unit 
initial area) and the product of —p by the ratio of dx dy to da db 
namely 1 — 2<rfc + (2/3 + The Youngs modulus E is deter 

mined (m accordance with the linear theory where the initial pressure 
IS different from zero) so that the linear parts of the t^^ o sides of the 
equation obtained m this way are equal Equating the coefficients of 
we obtain the equation 

2(X + p)0 + (X + 2m - p)« - -X(| - 2.T + V*) - 3 m 
- 1(1 - 2ff)* - 2ni(l - 9^) 

connecting the undetermined constants /3 and 5 The determinant 
of the two linear equations that 0 and 5 must satisfy is 

2(m - p)(3X + 2 m + p) 

and so 0 and 5 are determined unambiguously by the two equations 
unless p — M The theory breaks down when p = fi because the 
effective rigidity of the medium is m “ P and so the smaller is m ~ P 
the more liquid like (i e less able to support tension) is the medium 
It IS important to remember that fi is a. function of p it ma} n ell be 
that M “ P 15 for positive values of p, an increasing function of p 
and so M “ P IS never zero 
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1. The comprc'ision of a spherion! shell 


We consider nn ela'^tically isotropic spherical shell whose inlcrnal and 
external radii, when the shell is free from stress, are 7?, and Rr, respec- 
tively. The problem we wi^h to .solve is the following: What is the 
deformation of the .shell under an applied internal pressure p, and an 
applied external pressure p,? From reasons of symmetry each particle 
of the shell undcrgoe.s a radial displacement whose magnitude is a 
function of the initial distance r of the particle from the center of the 
shell. If wo denote this displacement by ku, the space polar coordi- 
nates (relative to a reference frame whose origin is at the center of 
the shell) of the final position of a particle whose initial coordinate.s 
had the space polar coordinates (r. 0, <f>) are (r ku, 0, <i>). Here k is 
a constant which is such that in the linear theory we neglect powers of 
k above the first whereas in the .second-order approximation we neglect 
powers of k above the .«ccond. I'he .space ])olar coordinates of the 
vector da arc (dr, rdO, r sin Od<i>), and the .space polar coordinates of 
the vector dx are (dr -f /; du, (r -{- ku) dO, (r -f ku) .sin 0 Thus 

da is the 3X1 matrix 


and dx is the 3 X 1 mat rix 



^dr -f k du 
(r -f ku) do 
^(r -r ku) sin 0 dd>) 


Hence the 3 X 3 matrix ./ that traii'-fonns, by means of the formula 
dz « J da, da into dx is the diagonal matrix 


J 


-f kid 
0 


itO 



0 


0 
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where the pnme indicates differentiation with respect to r It follows 
that the space polar coordinate of the strain matrix ij are furnished by 
the elements of the 3X3 diagonal matrix 



Neglecting powers of k above the second we have 

/, = k(u’+ 2fj + + 2^}' 



Since the medium is by hypothesis isotropic and free from initial 
stress 


= MiEa + 2/^7} •{- (f/,* - 2mh)Bs + 2mhy, -f* » co j? 

Stj 

(In this equation, <(> is, of course the energy density and not the third 
dA d<^ 

space polar coordinate ) Thus — is diagonal and so To = •/ T" ^ 

07} 

diagonal We assume that the shell is free from applied mass forces, 
and we find the conditions that v must satisfy by equating to zero the 
3X1 matnx diVa Ta* or, equivalently, since To* = To the 3X1 
matrix divo To Note The axes of the space polar coordinate refer- 
ence frame that is attached to the final position (r -h ku, 6, 4>) of any 
particle of the medium have the same directions as the axes of the 
space polar coordinate reference frame that is attached to the initial 
position (r, 6, <t>) of this particle In other words the rotation matrix 
Ri that transforms the latter of these two reference frames into the 
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former tlic identify mntrix Es. Thus the equation div« iTa*Ri*) = 
1), wlnrii express's that the Tna.«s, or body, force? arc zero, reduces to 
fiiv„ (7'„*} =- 0. In evaluating dhv. Tn we must use the expression for 
the divergence of a 3 X 3 matrix in space polar coordinates. This 
<'\prt'-~inn, and its derivation, arc given in Introduction to Applied 
Mothniuitici-. p. 112.’ For a diagonal 3 X 3 matrix who.?e diagonal 
dement? fr, 00, arc functions of r alone, the divergence reduces to 
the 3 X 1 matrix 

it' + “ f? — - <#><?', (Oi — o\ 

r r r r / 

and (irice the second and Ihiid diagonal elements of 'fa arc equal it 
follow.- that u must satisfy the single condition furni.shcd by the equa- 

tion rr' " rr (00 -f = 0. 

r r 

We first obtain the solution furnished by the linear theory by 
nt’plccling power? of k above the finst. 7’o reduces to the diagonal 
matrix whose diagonal elements arc 

-f- 2 2fiu' 




and so i< must sali-^fy the ditierential equation 


k(^ + 2 ;.) 




u" + 2- 



- 0 . 


Ihc general .solution of this equation is m = Ar + — where A and B are 

r~ 

eoustftnts which must be determined by the facts that 7' a dS''‘ hn.s the 
v: \w' — p, dS' when r — R; and the value — pc dS^ when r = Re- The 

3 X 1 matrix dX' is the product of the 3X1 matrix dS'^ by + k 

.Hs’onJ and third elements of each of these matrices being zero). 

A 

Ihnee .-1 and B are dt'tcrmined by setting ri- — —p{ ^1 -f k when 
R, and rf — p, -f k when r = /?,. Thu? p, and pe are 
’ -b t r.a,:- ]). ^Sl!nlARb.•i^!. John WiLy A i^n?, tOtS. 


r 
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of the same order of magnitude as the products of X and nhy k and, 
to the degree of approximation afforded by the linear theory, we may 

replace — p, and — + by — p, and — p„ 

respectively When we substitute « = Ar + ^ m the expression for 
fr we obtain 

S = il!j(3X + 2riX-^j. 

and so the two equations that serve to determine A and B are 
(3X + 2„)A-^-0=- 

Hence 

, ^ - V,K.’ (p, - v,)RW 

ki,a\ + 2;.)(R.’ - R.“)’ ” 4i^(R,’ - B,‘) 

We denote by Ui this linear theory solution, thus the radial displace- 
ment feui that is furnished by the linear theory is 

p,R‘ - p.B.» (p, - v,)R,V 1 

' ” (3X + 2p)(R.* - R’) MR,‘ - R,‘) r* 

The solution furnished by the second order approximation is of the 
form u = Ui -j- kw where, m determining the equation that w must 
satisfy, we may replace u by u\ in the coefficient of fc* in Ta (for ve 
neglect, in the second order approximation, powers of k above the 
second) Ta is the 3X3 diagonal matnx whose diagonal elements are 
furnished by the formulas 

rr = i: |x -f 2^^ + + ft® jx Q w'® + 

+ 3;iu'® + Z + 2 + 27rt + n j;j|’ 

^ = ft jx + 3 4- 2 fi + ft® jx Q 3 
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wiicrc it is understood that, in the coefficients of h', v is to be replaced 
by whereas, in the oocflicients of it is to be replaced by «i -f he. 
On caicnlntinf!; diva T„ we find that the coefficient of /; is zero (pre- 
ri‘-cly because i/j was determined by means of the linear theory'), and 
tm eqnafinp the coefficient of f:~ to zero we find that te must satisf}’ the 
equation 


(X + 2,0 + 2 y - 2 - 2(X + 3p + 2m) 




lienee le CV -f — + 
r 


J) (X + 3ft + 2m)B‘ 


1S(X + 3m + 2m)B- 


where C anrl D are two 


(X + 2m) J"* 

utidelermined constants. The radial di.splncement furnished by the 
-ecoiid-order approximation is 


hn -- ht^ + h-w = Q:A + k‘C)r + 


{l-B + IrD) 


+ 


(X + 3m + 2m)lrB’ 

(X + 2fi)r^ 

'I'lie constants .-I and B have been determined, by the linear theorj’, 
i-o that the coefficient of }: in fr ~ when r = /f,-, and — 

when r We now determine C and 1) so that the coefficient of k~ 

in r? ~-2pi~j} when r — and — —2p,~! when r ~ R^. For 
?h)s to he the case we must huv(> 

.3X d- 2 m)0 - - — - (a 

, fax -i' 1 I»m)(X + 'b* -r 2m) R" 


X + 2m 


OR- , /O GAR . 3n\ 


atal !! similar eqviatjoti obtained bv replaeinc R, bv R, and p, bv p,. 
- 2 ; 
k 


On repl’srinp — by (f»X fpjA — '"yj and collect ini; term.s wc 
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(3X 


+ 2y.)C - = ^|x + M A2 


+ (6X 4- 4;i + Cm - 


■(' 


>(7X + 10m - 8m) 
X + 2;i 




On replacing R, by we obtain a second equation connecting C and 
D, and on solving these two equations (or C and Dwe obtain 

A} + U(7X + 10/1 - 8m) n 1 
i(X + 2 m)(3X + 2^) 3X + 2;t) 


(1 _ 91 + n\ 

' V2 3X + 2^) 

■ - i ^3X + 2;i + 3m 




AB 4- 


(7X 4- lOfx - 8m) 
4(X 4- 2n) 


4- 




jSpcctal Cases 
1 The solid sphere 


pressure simply by p) A 


Here R, = 0 and so (if we denote the applied 

914-n> 


. B - 0, C 


fc»(3X 4- 2 m)* 


A(3X + 2m) 

The radial displacement is 

n 9f4-n \ p^ j 
■ V2 3X + 2m/ (3X 4- 2m)*| ’ 


/I 9l + n \ 
A2 3X 4- 2m/ 


l3X 4- 2 m 

It IS easy to verify that this expression agrees with the result furnished 
by the second order approximation m the treatment of the hydrostatic 
compression of an isotropic medium (carry out this verification) 
{Hint If the strain matrix is —eEs the coefficient of r above =* — c — 
4^* since, if we denote this coefficient by (1 4- 5)* = (1 — 2e)^) 

2 The spherical caitly in an infinite solid Here pj = 0, Rj = * 
and so (if w e denote the pressure on the surface of the cavity simply 
by p and the initial radius of the cavity simply by R) A = 0, B = 

The radial 


4A:m ^ '' 
displacement is 


0, B = 

[£ + 


7X 4~ 10m - Bot 
4(X4-2m) ■' 

px 4- IQm - Sm 
1 4(X4-2m) 




3 The thin spherical shell Denoting the thickness of the shell bj 
e, we have B* = Ki 4- « On keeping m the expressions for A and B 

only the lowest powers of ■—> we obtam 



THi: (’fJMJ^HKSSIOX OF A CIRCULAR TUBE 


12.3 


I t 0>i - }h) jtR __ O’ij: 7>:* II2 - 

''\?(.1X -f 2,0 T /?,-*“ ”' 'i 2 m’*' t ' ' 


<nv; flicrt'forc, 7.vl =- It foU(mi> from the exprc-'-sions for 

3X -f 2,1 

(' fituf O (hii( 

{/} 9l + n\ ir,,«= pi(7X + 10,1 - Sm) 

" |\2 “ 3X (3X + (X -f 2„)(3X + 2,;) 


-f U-. 

^ 3X + 2m) 


-f37 

3X 4- 2m V 


:1X -I- 2m + ^tn 


27 ‘ "ifx+vr 


r,\ ^■■ 

Tin'- rndtn! tlKplaccmont is found by .substituting tlic.s(' vniucs in the 
«'\pro} -lion 

, . j.^r u i. + 2n!)A--7F-j 
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whose magnitude is a function of the initial distance r of the particle 
from the axis of the tube If ViC denote this displacement by ku, 
the plane polar coordinates (relative to a reference frame ■whose origin 
IS on the axis of the tube) of the final position of a particle ivhose initial 
position had the plane polar coordinates (r 0) are (r -}* ku 6) Here 
A, IS a constant which is such that in the linear theory we neglect powers 
of k above the first whereas in the second order approximation we 
neglect powers of k above the second The plane polar coordinates 
of the vector da are (dr, r d6), and the plane polar coordinates of the 
vector dx are (dr + k du, (r + ku) dd) Thus da is the 2X1 matrix 

( whereas da: is the 3X1 matrix 2X2 

matrix J that transforms by means of the formula dx ~ J da, da into 
dx IS the diagonal matrix 

/I + Au 0 \ 

where the pnme denotes differentiation with respect to r It follows 
that the plane polar coordinates of the (plane) strain matrix i; are 
furnished by the elements of the 2X2 diagonal matrix 



Neglecting powers of k above the second we have 

r 

Since the medium is by hjpothesis isotropic and free from initial 
stress, 

— \I1E3 -j- 2ftri + (f/i* — 2niI'>)E3 + 2/n/iJj 
drj 

Note Since we are dealing with a two-dimensional problem so that 
the strain matrix »j is two-dimensional rather than three dimensional, 
the invariant 1 3 does not appear and we have only the two third order 
elastic constants I and m rather than the three third order elastic 
constants I, m, and n 
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~ nnd Ta ^ — sre diagonal. We assume that our tube is free 

bn 3r} 

5 rr>m applied mass forces and so we find the conditions that w must 
^atis^y by setting divn To = 0 (why?). For a diagonal 2X2 matrix 
iviiO’C diagonal elements fr, Off are functions of r alone, the divergence 
(in plane polar coordinates) reduces to the 2 X 1 matrix 

^5 t ' + ^ (tt — fO), 0^ 

(tee Introdnrh'ou to Applied Mathematics, p. 98)’ and therefore v must 
F!iti«fy the single condition furnished 1 ) 3 ' the equation fr' + " (fr — 00) 
f- n. 

We first obtain the linear theor 3 ' solution b 3 ‘ neglecting powers of I: 
alutvc the first. Tn reduces, then, to the 2X2 diagonal matrix whoso 
diagonal elements are 


rr 




(X -j- "b X ~ 


oh jxn' + (X + 2m) ’-'j 
and so m must satisf}* the differentia! equation 


ru" + tt' - - - 0. 
r 


/? 


T)i(> general solution of this equation is u = Ar H where A and B 

r 

arc constants which must bo determined b.v the facts that T„ dS° has 
tl'.f’ v.nhie —pidS' wlion r ~ /?,• nnd the value —Pr dS‘ when r — 7?,. 
The 2 X 1 m.atrix dS' is the product of the 2X1 matrix dS''- bi* 1 -f 

t thus .1 and /i are determined b 3 * setting fr ~ —Pi(^) 

v.'nen r - /?,■ nnd ff + k when r — /?,. It follow.s that 

.nnd p. are of the stime order of magnitude as the products of X 
ntid fs in- k. and so (to the degree of approximation afforded by’ the 
litirsr theory) we may' equate rf to — p,- when r - /?,- and to — p. 


v.y. 


B 


on r w \\licn we substitute v Ar -|- •" in the cxprc.s.«ion for 


Fifterb 1). Mura.»it;5i:in, John ttltcy A: Sterna, lOtS. 
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— T— 5 J D = 0 The radial displacement is 

4fe2(X + u)* 

-P^ 11 + ^2 + ! 

2(X + |x)l ^ X+/i/2(X + M)j 


2 The thin-walled circular tube Denoting the thickness of the walls 
of the tube by e we have R* = R» + « On keeping m the expres- 
sions for A and B only the lowest powers of > we obtain 
ic, 


lA 


Pi - Vt R. 
4(X-|-m) «’ 


^ ^ (p» - Vt) ^ 
R* An ( 


say, and therefore kA ‘ 
formulas 


X + p 


0 C and D are furnished by the 




k^D - - 


j3X + + Am 2(n + 2/n)) 


\ + ti \+2n} 

The radial displacement of any point on the inner wall of the thin tube 


(iu)r-*, 


'-M- 

-X + M 


3X + Gm 4- Am ^ 
X + M 


nin + 2m) 

(X + ^i)(X 4* 2/i) 


4- Cm — X] 

2 (X + 2 u) I 


B, 


3. The torsion of a circular cylinder 

We choose the axis of our c^ linder as the axi'* of a system of cylindrical 
coordinates (r, d, z), and we denote by r 4- fc«, ^ 4- kz, z + kw the 
final coordinates of the particle of the cylinder whose initial coordinates 
were (r, $, i) Here k is e. coiatant parameter which is such that its 
square and higher pow ers may be neglected in the linear theory whereas 
its cube and higher powers may be neglected in the second-order 
approximation We assume that u is a function of r alone, that ta is a 
function of z alone, and that the medium is isotropic and free from 
initial stress, and we try to so determine the unknown functions u “ 
n(r), w = v(z) that 



TOHSKLX or A CIRCULAR CVLIKDER 


i;u 


(1) Thf* (Irformrd cylinder is in equilibrium under the action of 
^t!r^nfe forces alone, the mass* or body forces being zero. 

(2? 7’here are no applied forces on the sides (or curved portion of 
the bounding surface) of the cylinder, and the forces applied to eatdi 
cml of the cylinder redtice to a couple about the axis of the cylinder. 

The cylindrical coordinate.*: of the vector da arc (dr, r dO, dz) whcrca.s 
the cylindrical coordinates of the vector dx arc {(I + hd) dr, (1 -{- 

k~) {rdO -t- krdz), (1 }:ir') dz\ where a' denotes the derivative of 

r 

H with respect to r and ic' denotes the derivative of ic with respect to 
z. Kftli'. The cylindrical coordinate*: of dx have reference to a 
fr.atne oldained hy rotating the reference frame in which the cylindrical 
courdin.ates of da are (dr, r dO, dz) through an angle kz around the z-axis. 
The S X 3 matrix that transforms, by means of the formula dx ~ J da, 
ih<- 3 X 1 matrix 

/ dr \ 


into the 3 X I matrix 


da -- I r do 
\ dz 


(1 -r kid) dr ' 

d- k~jrdO-\- k(r + ku)dz 
(1 -f kid) dz 


/ I d- kid 0 ^ \ 

d- 0 Id-/.-‘ /.-(rdWm) • 

'0 0 id- kid / 

Hence th.e cylindrie.'d coordinates of the strain matrix rj arc given by 
llic formulii 


id 0 0 
.. » 1 


“r n 


f«')- 0 


(le'j* d- 


tpor.en- o: above the .-econd b-ins; nej^ioctcdL It follows that 

/, . k(‘d 4- 'i 4. .A 4. r- f,P*‘ • ‘ 4 ^ , -I 

\ r ' ' J ' \ 2 ' *^57 + 5 r d- s- . 
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1 , 
--ru 


We first determine the linear theorj solution Here 



«' H h w' ) + 


\ 4, H 


+ 2/i - 




+ w')+2Mtx>V 


and To = J T~ 13 the same as this Since wo ere using orthogonal 
01 ) 

curvilinear coordinates the matrix ivhose divergence must bo equated 
to zero (to express the fact that the mass or body, forces are zero) 
is Ta*J2i* nhere I2i is the rotation matrix that transforms the cylin- 
drical coordinate reference frame at (r, $, z) into the cylindrical coor- 
dinate reference frame at {r + ku, 9 + Arz, 2 + kw) 

( cos kz — sin kz 0\ /O 0\ 

sin A 2 cos Az 0 j = £?} + A; I 2 0 0 J + 

0 0 1 / \0 0 0 / 

Since we are neglecting (when obtaining the linear theory solution) 

powers of k above the first the matrix •whose divergence must be set 
equal to zero 13 





+ 2p- 




Using the formula for the divergence of a 3 X 3 matrix in cylindrical 
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coordinates (fgc Inirodmlion io Applied MaOicvmlm, p. 112),^ we 
obtain tlic following equations; 

^ = 0 , w" = 0 . 

r r“ 


lienee n is a linear combination of r and - and, since n must be bounded 

r 

at r ” 0, it follows that v is a multiple of r. Also w' is a constant and 
so we maj' take, without loss of generality (why?), tu to be a multiple 
of 2 ; 

u — At, w — liz. 


TIk* fact that there is no apjdicd force on the sides of tlic cylinder 
yields 2(X + n)A + X/i — 0. The coordinate in the direction of 
the axis of the cylinder of the force, per unit initial area, on the end 
z I of the cylinder is constant and equals 2X/1 + (X + 2p.)B. In 
order that the force on the end 2 = / of the cylinder should reduce to a 
couple we mvist have, (hen, 2XA + (X -f 2;i)/i = 0. Tlic two cqua- 
tiniis 2(X n)A XR = 0, 2X/1 (X -}* 2^)R = 0 yield A — 0, 
li =" 0, Thus the solution furnished by the linear thcor>' is vi = 0, 
it'i 0; thus any cross section z = constant of the cylinder is merely 
rotrderl through an angle kz. The force, per unit initial area, applied 
to the end z ~ I of the cylinder is perpendicular to the radius to the 
jioint of application and is of magnitude kiir. Hence the system of 
forces applied to the end s ~ / of the cylinder is equivalent to a couple 


or torque around the axis of the cxdinder of amount 



If o is the angle through which the end s = 



I of the cylinder 


is rotated we have « — hi and thus we obtain the formula 


Tortjue — 


a^r^tR^ 

~"2r 


v-hirh Ff-rves to determine the 77ioduliis of ripiditij n. 

In order to obtain the solution furni.shed by tlie second-order 
approximation, we set, since ui and trj are zero, it — h\ and ir k<-' 
''herec and j art' undetermined functions of randr. respectively; 

c r(rU >.(z}. 

Ne;’!cc{ing pov-ers of k above the s<‘cond, wc obtain 
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/O 0 0\ 

jj = fc I 0 0 ) 

\0 ir 0 / 


'V 0 0 

-H I 0 0 

^0 0 fi + 


where v denotes the derivative of v with respect to r and s denotes the 
derivative of s with respect to z Hence 


and s 


h - -ifeV 

/r’ 0 0\ 
co:j = -iJkMO 0 0 
\0 0 0 / 


'0 0 
0 0 
^0 nr 


9 


+ k^ 


+ 2nv 


X + 3 ^ 

+ ( 5 ’' + 2 '"-^) 




+ 2,1 


+ 2,iS 


X + ^ + s ^ 

+ (2^ + "+2’“) 


J 

The matrix whose divergence must be equated to zero (to express the 

fact that the mass or body forces are zero) is — J*Ri* ^\here 

di? 

/O -z 

Ri = B 3 + fcl2 0 0j+ f J ~ Ez + k 




/O 0 0\ 

0 0 r - 

\0 0 0/ 
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Hence 

/O -2 0\ 

HiJ ^Et + h[z 0 r + • 

\0 0 0 / 

/ 0 2 0 \ 


J*U\* - Fa + A- 1 —2 0 0 I + • • • and so 
\ 0 r 0/ 


0 0 0 


0 r 

Vo r 0 


X + |; + s'^ + Uv' 

(\ 1 IN 

+ (^-X + -m--7.j 


+ G ''+'‘+ r ") 


X 4- " + s'^ + 2/is' 

+ Qx + „+i„,)r= 


On equatinR the divergence of this matrix to zero, ve obtain the 
foiknving equations: 


(X + 2ft) ~ = ^2/1 — X - f« + ^ 71^ r, s" ~ 0. 

Since (' is bounded at r = 0, we have 

^ (2m - X - 7?t 4- f77)r* 

Cr + — 

whore C is an undetermined constant, .s is a linear function of z, 
and we may set (why?) 

s Dz 


■vhere D h an undetermined constant. Since - ./ ~ is given bv 

3t) 


the furmula 
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/O 0 0\ 

. « i 0 0 + 

\0 fir 0/ 

X ^ 4- s'^ 4- 2nv' 


+ (5’' + '‘ + 5“)’- 


X ^t' 4- ~ 4" + 2/1^ 

+ G ''+''+ i”')'' 


the fact that there is no applied force on the sides of the cylinder yields 

2(\ + rfC + XD + |20X>. + 24/ + 4^m + (2X + «)»)= 0 

^here R is the initial radius of the cjlinder The axial coordinate of 
the force per unit initial area applied to the end a « f of the cylinder is 

the product of 1' by 2XC + (X + 2^)D + ± 8t»" + 3Xn 

8(X 4" 2 ft) 

and -we choose the constants C and D bo that the integral of this over 
the end of the cylinder is zero Thus \\c obtain the relation 

2XC + (X + 2^)D + g jX,.+ + 

1G(X + 2;*) 

U'’hen C and D satisfy the two linear equations just determined, the 
force applied to the end z = f is a couple about the axis of the c> linder 

of amount (pro\e this) Under this moment (or torqiie) not 

onlj docs the c>lindcr rotate (as predicted by the linear theory) but 
also its radius increases by the amount 

K^R jc 4- /^sj 

1 8(X4-2;i) I 
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nnd its Icnplh increases by tlie amount JrDl. Since C and D are 
multiples of I?" it follows that the increase of the radius of the cylinder 
is n multiple (possibly negative) of JrR^ AYhereas the increase of the 
length of the cylinder is a multiple (possibly negative) of IrR-. Accur- 
ate measurements of the changes of the radii and lengths of cylinders 
(of reasonably large radius) would serve to detennine the third-order 
elastic constants m and n (it being supposed that the second-order 
elastic constants X and n arc already known). 
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